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INTRODUCTION 


A new method has been developed for calculating the 
pressure distribution and aerodynamic characteristics (lift, 
drag, and pitching moment) of wing-body-tail combinations in 
subsonic and supersonic potential flow. This report describes 
the computer program developed to perform the numerical calcu- 
lations . 


METHOD OF SOLUTION 


The configuration surface is subdivided into a large 
number of panels, each of which contains an aerodynamic singu- 
larity distribution. A constant source distribution is used 
on the body panels , and a vortex distribution having a linear 
variation in the streamwise direction is used on the wing and 
tail panels. The normal components of velocity induced at 
specified control points by each singularity distribution are 
calculated and make up the coefficients of a system of linear 
equations relating the strengths of the singularities to the 
magnitude of the normal velocities. 

The singularity strengths which satisfy the boundary con- 
dition of tangential flow at the control points for a given 
Mach number and angle of attack are determined by solving this 
system of equations using an iterative procedure. Once the 
singularity strengths are known, the pressure coefficients are 
calculated, and the forces and moments acting on the configu- 
ration determined by numerical integration. A detailed descrip- 
tion of the method is given in Part I of this report. 


PROGRAM DESCRIPTION 


The computer program is written in CDC FORTRAN IV, version 
2.3 for a SCOPE 3.0 operating system and library file. It is 
designed for the CDC 6000 series of computers, occupies 70,000 
(octal) words, and operates in OVERLAY mode. The program re- 
quires five peripheral disc files in addition to the input and 
output files. 



PROGRAM INPUT DATA 


The input to this program consists of two basic parts, 
namely, the numerical description of the configuration geometry 
as adapted from reference 1, and an auxiliary data set speci- 
fying the singularity paneling scheme, program options, Mach 
number, and angle of attack. The program input is illustrated 
by the sample case presented in Appendix III. 


Description of Input Geometry Cards 

The configuration is defined to be symmetrical about the 
xz plane, therefore only one side of the configuration need be 
described. The convention used in this program is to present 
that half of the configuration located on the positive y side 
of the xz plane . The number of input cards depends on the 
number of components used to describe the configuration, and 
the amount of detail used to describe each component. 

Card 1 - Identification .- Card 1 contains any desired iden- 
tifying information in columns 1-80. 

Card 2 - Control integers .- Card 2 contains 24 integers, 
each punched right justified in a 3-column field. Columns 


73-80 may 

be used in 

any desired manner. Card 2 contains the 

following 




Columns 

Variable 

Value 

Description 

1-3 

JO 

0 

No reference area 



1 

Reference area to be read 

4-6 

J1 

0 

No wing data 



1 

Cambered wing data to be read 



-1 

Uncambered wing data to be read 

7-9 

J2 

0 

No fuselage data 



1 

Data for arbitrarily shaped fuse- 
lage to be read 



-1 

Data for circular fuselage to be 
read (With J6=0, fuselage will be 
cambered. With J6=-l, fuselage 
will be symmetrical with xy-plane. 




With J6=l, entire configuration 
will be symmetrical with xy-plane) 

10-12 

J3 

0 

No pod (nacelle) data 



1 

Pod (nacelle) data to be read 
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Columns Variable Value 


Description 


13-15 J4 0 No fin (vertical tail) data 

1 Fin (vertical tail) data to be 

read 

16-18 J5 0 No canard (horizontal tail) data 

1 Canard (horizontal tail) data to 

be read 

19-21 J6 0 A cambered circular or arbitrary 

fuselage if J2 is nonzero 
1 Complete configuration is sym- 

metrical with respect to xy-plane, 
which implies an uncambered circu- 
lar fuselage if there is a fuse- 
lage 

-1 Uncambered circular fuselage with 

J2 nonzero 

22-24 NWAF 2-20 Number of airfoil sections used 

to describe the wing 

25-27 NWAFOR 3-30 Number of ordinates used to define 

each wing airfoil section. If the 
value of NWAFOR is input with a 
negative sign, the program will 
expect to read lower surface 
ordinates also 

28-30 NFUS 1-4 Number of fuselage segments 

31-33 NRADX(l) 3-30 Number of points used to represent 

half-section of first fuselage 
segment. If fuselage is circular, 
the program computes the indicated 
number of y- and z-ordinates 


34-36 

NFORX (1) 

2-30 

Number of stations 
lage segment 

for 

first fuse 

37-39 

NRADX(2) 

3-30 

Same as NRADX (1) , 
fuselage segment 

but 

for 

second 

40-42 

NFORX (2) 

2-30 

Same as NFORX (1) , 
fuselage segment 

but 

for 

second 

43-45 

NRADX ( 3 ) 

3-30 

Same as NRADX (1) , 
fuselage segment 

but 

for 

third 
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Columns 


Variable 


Value Description 


46-48 

NFORX(3) 

2-30 

Same as NFORX (1) , but for third 
fuselage segment 

49-51 

NRADX ( 4 ) 

3-30 

Same as NRADX(l), but for fourth 
fuselage segment 

52-54 

NFORX ( 4 ) 

2-30 

Same as NFORX (1) , but for fourth 
fuselage segment 

55-57 

NP 

0-9 

Number of pods described 

58-60 

NPODOR 

4-30 

Number of stations at which pod 
radii are to be specified 

61-63 

NF 

0-6 

Number of fins (vertical tails) 
to be described 

64-66 

NFINOR 

3-10 

Number of ordinates used to 
describe each fin (vertical tail) 
airfoil section 

67-69 

NCAN 

0-2 

Number of canards (horizontal 
tails) to be described 

70-72 

NCANOR 

3-10 

Number or ordinates used to define 


each canard (horizontal tail) air- 
foil section. If the value of 
NCANOR is input with a negative 
sign, the program will expect to 
read lower surface ordinates also, 
otherwise the airfoil is assumed 
to be symmetrical 

Cards 3, 4, ... - remaining input data cards . - The remain- 
ing input data cards contain a detailed description of each com- 
ponent of the configuration. Each card contains up to 10 values, 
each value punched in a 7-column field with a decimal point and 
may be identified in columns 73-80. The cards are arranged in 
the following order: reference area, wing data cards, fuselage 

data cards, pod data cards, fin (vertical tail) data cards, and 
canard (horizontal tail) data cards. 

Reference area card: The reference area value is punched 

in columns 1-7 and may be identified as REFA in columns 73-80. 

Wing data cards: The first wing data card (or cards) con- 

tains the locations in percent chord at which the ordinates of 
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all the wing airfoils are to be specified. There will be exact- 
ly NWAFOR locations in percent chord given. Each card may be 
identified in columns 73-80 by the symbol XAFJ where J denotes 
the last location in percent chord given on that card. 

The next wing data cards (there will be NWAF cards) each 
contain four numbers which give the origin and chord length 
of each of the wing airfoils that is to be specified.^ The 
card representing the most inboard airfoil is given first, 
followed by the cards for successive airfoils. These cards 
contain the following: 


Columns 

Contents 


1-7 

x-ordinate of airfoil 

leading edge 

8-14 

y-ordinate of airfoil 

leading edge 

15-21 

z-ordinate of airfoil 

leading edge 

22-28 

airfoil streamwise chord length 

73-80 

card identification, WAFORGJ where 
J denotes the particular airfoil. 


thus WAFORGl denotes the most 
inboard airfoil 

If a cambered wing has been specified, the next set of wing 
data cards is the mean camber line cards. There will be NWAFOR 
values of delta z referenced to the z-ordinate of the airfoil 
leading edge, each value corresponding to a specified percent 
chord location on the airfoil. These cards are arranged in 
the order which begins with the most inboard airfoil and pro- 
ceeds outboard. Each card may be identified in columns 73-80 
as TZORDJ where J denotes the particular airfoil. Note that 
the z-ordinates are dimensional. 

Next are the wing ordinate cards. There will be NWAFOR 
values of half-thickness specified for each airfoil expressed 
as percent chord. These cards are arranged in the order which 
begins with the most inboard airfoil and proceeds outboard. 

Each card may be identified in coliamns 73-80 as WAFORDJ where 
J denotes the particular airfoil. 

Fuselage data cards: The first card (or cards) specifies 

the X values of the fuselage stations of the first segment. 

There will be NFORX(l) values and the cards may be identified 
in columns 73—80 by the symbol XFUSJ where J denotes the number 
of the last fuselage station given on that card. 
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the fuselage is circular/ the next card (or cards) gives 
the fuselage cross sectional areas, and may be identified in 
columns 73~80 by the symbol FUSARDJ where J denotes the number 
of the last fuselage station given on that card. If the fuse- 
lage is of arbitrary shape, NRADX(l) values of the y-ordinates 
for a half-section are given and identified in columns 73-80 
as YJ where J is the station nvimber. Following the y-ordinates 
are the NRADX(l) values of the corresponding z-ordinates for 
the half-section identified in columns 73-80 as ZJ where J is 
the station number. Each station will have a set of y and z, 
and the convention of ordering the ordinates from bottom to 
top is observed. 

For each fuselage segment a new set of cards as described 
must be provided. The segment descriptions should be given in 
order of increasing values of x. 


Pod data cards: The first pod (nacelle) data card speci- 

fies the location of the origin of the first pod. The card 
contains the following: 


Columns 

Contents 




1-7 

x-ordinate of origin 

Of 

first 

pod 

8-14 

y-ordinate of origin 

of 

first 

pod 

15-21 

z-ordinate of origin 

of 

first 

pod 

73-80 

card identification, PODORGJ 
where J denotes the pod number 


. ^he next pod input data card (or cards) contains the x— 
ordinates, referenced to the pod origin, at which NPODOR values 
of the pod radii are to be specified. The first x value must 
be zero and the last x value is the length of the pod. These 
cards may be identified in columns 73-80 by the symbol XPODJ 
where J denotes the pod number. 

For each additional pod, new PODORG, XPOD, and POOR cards 
must be provided. Only single pods are described but the pro- 
gram assumes that if the y-ordinate is not zero an exact dupli- 
cate is located symmetrically with respect to the xz-plane, a 
y-ordinate of zero implies a single pod. 

Fin data cards : Exactly three data input cards are used 

to describe a fin (vertical tail) . The first fin data card 
contains the following: 
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Coliamns 


Contents 


x-ordinate on inboard airfoil 
leading edge 

y-ordinate of inboard airfoil 
leading edge 

z-ordinate of inboard airfoil 
leading edge 

chord length of inboard airfoil 

X— ordinate of outboard airfoil 
leading edge 

y-ordinate of outboard airfoil 
leading edge 

z— ordinate of outboard airfoil 
leading edge 

chord length of outboard airfoil 

card identification, FINORGJ where 
J denotes the fin number 

The second fin input data card contains NFINOR values of 
X expressed in percent chord at which the fin airfoil ordinates 
are to be specified. The card may be identified in columns 
73-80 as XFINJ where J denotes the fin number. 

The third fin input data card contains NFINOR values of 
the fin airfoil half-thickness expressed in percent chord. 

Since the fin airfoil must be symmetrical, only the ordinates 
on the positive y side of the fin chord plane are specified. 

The card identification FINORDJ may be given in columns 73-80 
where J denotes the fin number. 

For each fin, new FINORG, XFIN, and FINORD cards must be 
provided. Only single fins are described but the program 
assumes that if the y-ordinate is not zero an exact duplicate 
is located symmetrically with respect to the xz-plane, a y 
ordinate of zero implies a single fin. 

Canard data cards: If the canard (or horizontal tail) 

airfoil is symmetrical, exactly three cards are used to describe 
a canard, and the input is given in the same manner as for a 
fin. If, however, the canard airfoil is not symmetrical 


1-7 

8-14 

15-21 

22-28 

29-35 

36-42 

43-49 

50-56 

73-80 
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(indicated by a 
input data card 
The information 
is as follows: 

negative value of NCANOR) , a fourth canard 
will be required to give the lower ordinates, 
presented on the first canard input data card 

Columns 


Contents 

1-7 


x-ordinate of inboard airfoil 
leading edge 

8-14 


y-ordinate of inboard airfoil 
leading edge 

15-21 


z-ordinate of inboard airfoil 
leading edge 

22-28 


chord length of inboard airfoil 

29-35 


x-ordinate of outboard airfoil 
leading edge 

36-42 


y-ordinate of outboard airfoil 
leading edge 

43-49 


z-ordinate of outboard airfoil 
leading edge 

50-56 


chord length of outboard airfoil 

73-80 


card identification, CANORGJ 
where J denotes the canard number 


The second canard input data card contains NCANOR values of 
X expressed in percent chord at which the canard airfoil ordi- 
nates are to be specified. The card may be identified in 
columns 73—80 as XCANJ where J denotes the canard number. 

The third canard input data card contains NCANOR values of 
the canard airfoil half-thickness expressed in percent chord. 
This card may be identified in columns 73-80 as CANORDJ where 
J denotes the canard number. If the canard airfoil is not sym- 
metrical, the lower ordinates are presented on a second CANORD 
card. The program expects both upper and lower ordinates to 
be punched as positive values in percent chord. 

For another canard, new CANORG, XCAN, and CANORD cards 
must be provided. 
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Description of Auxiliary Input Cards 


Card 1.1 - Identification .- Card 1.1 contains any desired 
identifying information in columns 1-80 . 

Card 1.2 - Boundary condition and control p oint definition. - 
Non planar boundary conditions are always applied on a body, 
hoLve? SLd 1.2 permits the selection of boundary condrtrons to 
appiron a wing, fin (vertical tail), or canard (horizontal tail) 
This card also^selects the output print options. This card con 
tains the following: 


Columns Variable Value 
1-3 LINBC 0 


1 


4-6 THICK 0 

1 

7-9 PRINT 0 

1 

2 

3 

4 


Description 

Control points on surface of wing, 
fin (vertical tail) , and canard 
(horizontal tail) . This is re- 
ferred to as the nonplaner bound- 
ary condition option. 

Control points in plane of wing , 
fin (vertical tail) , and canard 
(horizontal tail) . This is 
referred to as the planar bound- 
ary condition option. 


Do not calculate wing thickness 
matrix 

Calculate wing thickness matrix 
if LINBC = 1 

Print out the pressures and the 
forces and moments 

Print out option 0 and the span- 
wise loads on the wing , fins , and 
canards 

Print out option 1 and the veloc- 
ity components and source and 
vortex strengths 

Print out option 2 and the steps in 
the iterative solution 
Print out option 3 and the axial 
and normal velocity matrices 


A negative value of print adds the panel geometry print out 
to the output indicated for options 1 through 4 • 

LINBC, THICK, and PRINT are punched as right justified 
integers. THICK is not used if LINBC = 0. 
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Card 2.1 - Revised configura tion 
control integersT - The contents of car 
rxght justified integers as follows: 


Paneling description 
u ^ . 1 are punched as 


Colxomns Variable 
1-3 KO 

4-6 K1 

7-9 K2 

10-12 K3 

13-15 K4 


16-18 K5 


19-21 K6 


Value 

0 

1 

0 

1 

3 


0 

1 


0 

1 

3 


0 

1 

3 


Description 

No reference lengths 
Reference length data to be read 

No wing data 

Wing data to be read, wing has a 
sharp leading edge 
Wing data to be read, wing has a 
round leading edge 

No body data 
Body data follows 

Not used 

No fin (vertical tail) data 
Fin (vertical tail) data to be 
read , fin has a sharp leading 
edge 

Fin (vertical tail) data to be 
read, fin has a round leading 
edge 

No canard (horizontal tail) data 
Canard (horizontal tail) data to 
be read, canard has a sharp 
leading edge 

Canard (horizontal tail) data to 
be read, canard has a round 
leading edge 

Not used 


22-24 KWAF 


25-27 KWAFOR 


0, Number of wing sections used to 

2- 20 define the inboard and outboard 

panel edges. If KWAF = 0, the 
panel edges are defined by NWAF in 
the geometry input 

0 ! Number of ordinates used to 

3- 30 define the leading and trailing 

edges of the wing panels . if 
KWAFOR =0, the panel edges are 
defined by NWAFOR in the geometry 
input 
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Columns Variable Value 


Description 


28-30 KFUS 


31-33 KRADX(l) 0, 

3-20 


34-36 

KFORX(l) 

0, 

2-30 

37-39 

KRADX(2) 

0, 

3-20 

40-42 

KF0RX(2) 

0, 

2-30 

43-45 

KRADX ( 3 ) 

0, 

3-20 

46-48 

KFORX ( 3 ) 

0, 

2-30 

49-51 

KRADX ( 4 ) 

0, 

3-20 

52-54 

KFORX ( 4 ) 

0, 

2-30 


The number of fuselage segments. 
The program sets KFUS = NFUS 

Number of meridian lines used to 
define panel edges on first body 
segment. There are three options 
for defining the panel edges. 

If KRADX(l) = 0, the meridian 
lines are defined by NRADX(l) in 
the geometry input. If KRADX(l) 
is positive, the meridian lines 
are calculated at KRADX(l) equally 
spaced PHIKs . If KRADX(l) is 
negative, the meridian lines are 
calculated at specified values of 
PHIK 

Number of axial stations used to 
define leading and trailing edges 
of panels on first body segment. 

If KFORX(l) = 0, the panel edges 
are defined by NFORX(l) in the 
geometry input 

Same as KRADX(l) , but for second 
body segment 

Same as KFORX(l), but for second 
body segment 

Same as KRADX(l), but for third 
body segment 

Same as KFORX(l), but for third 
body segment 

Same as KRADX(l), but for fourth 
body segment 

Same as KFORX(l), but for fourth 
body segment 


The program is restricted to 600 body singularity panels. 
For this program there is an additional restriction that t e 
total number of singularity panels in the axial direction 
body (fuselage) cannot exceed 30. The arbitrary body (fuselage) 
capability of this program is limited to those shapes for whic 
the radius is a single-valued function of PHIK for each cross 
section of the body. 
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^ Card 2.2 - Addition al revised configuration panelinq 
cont r ol integers .- The conhf^n<-fi of o 7 

punched as right justified integers as follows: 


Columns 

Variable 

Value 

Description 

1-3 

KF(1) 

0 , 

2-20 

Number of fin sections used to 
define the inboard and outboard 
panel edges on the first fin. 

If KF(1) =0, the root and tip 
chords define the panel edges 

4-6 

KFINOR(l) 

0 , 

3-30 

Number of ordinates used to 
define the leading and trailing 
edges of the fin panels on the 
first fin. If KFINOR(l) = 0, the 
panel edges are defined by NFINOR 

7-9 

KF(2) 

0 , 

2-20 

Same as for KF(1), but for second 
fin 

10-12 

KFINOR(2) 

0 , 

3-30 

Same as for KFINOR(l) , but for 
second fin 

13-15 

KF(3) 

0 , 

2-20 

Same as for KF(1), but for third 
fin 

16-18 

KFINOR(3) 

0 , 

3-30 

Same as for KFINOR(I) , but for 
third fin 

19-21 

KF(4) 

0 , 

2-20 

Same as for KF(1), but for fourth 
fin 

22-24 

KFINOR(4) 

0 , 

3-30 

Same as for KFINOR(l) , but for 
fourth fin 

25-27 

KF(5) 

0 , 

2-20 

Same as for KF(1), but for fifth 
fin 

28-30 

KFINOR(5) 

0 , 

3-30 

Same as for KFINOR(I) , but for 
fifth fin 

31-33 

KF( 6 ) 

0 , 

2-20 

Same as for KF(1) , but for sixth 
fin 

34-36 

KFINOR( 6 ) 

0 , 

3-30 

Same as for KFINOR(l) , but for 
sixth fin 
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Columns 


Variable 


Value Description 


37-39 

KCAN ( 1 ) 

0, 

2-20 

Number of canard sections used to 
define the inboard and outboard 
panel edges on the first canard. 
If KCAN(l) = 0, the root tip 
chords define the panel edges. 

If KCAN(N) negative, no vortex 
sheets carry through the body and 
concentrated vortices are shed 
from the inboard edge of the 
canard or tail surface 

40-42 

KCANOR ( 1 ) 

0, 

3-30 

Number of ordinates used to define 
the leading and trailing edges of 
the first canard. If KCANOR (1)=0, 
the panel edges are defined by 
NCANOR 

43-45 

KCAN ( 2 ) 

0, 

2-20 

Same as for KCAN(l) , but for 
second canard 

46-48 

KCANOR (2) 

0, 

3-30 

Same as for KCANOR (1) , but for 
second canard 

49-51 

KCAN ( 3 ) 

0, 

2-20 

Same as for KCAN(l), but for 
third canard 

52-54 

KCANOR ( 3 ) 

0, 

3-30 

Same as for KCANOR (1) , but for 
third canard 

55-57 

KCAN ( 4 ) 

0, 

2-20 

Same as for KCAN(l), but for 
fourth canard 

58-60 

KCANOR (4) 

0, 

3-30 

Same as for KCANOR ( 1) , but for 
fourth canard 

61-63 

KCAN ( 5 ) 

0, 

2-20 

Same as for KCAN(l), but for 
fifth canard 

64-66 

KCANOR (5) 

0, 

3-30 

Same as for KCANOR (1) , but for 
fifth canard 

67-69 

KCAN (6) 

0, 

2-20 

Same as for KCAN(l), but for 
sixth canard 

70-72 

KCANOR (6) 

0, 

3-30 

Same as for KCANOR(l), but for 
sixth canard 
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The program is restricted to a total of 600 singularity 
panels on the wing-fin-canard combination. 

For this program there is an additional restriction that 
the total number of singularity panels in the spanwise direction 
on the wing-fin-canard combination cannot exceed 20. 


Cards 3,4,...- remaining input data cards . - The 
remaining input data cards contain a detailed description of 
the singularity paneling of each component of the configuration. 
Each card contains up to 10 values, each value punched in a 
7-column field with a decimal point and may be identified in 
columns 73-80. The cards are arranged in the following order: 
reference lengths, wing data cards, fin (vertical tail) data 
cards, canard (horizontal tail) data cards, fuselage (body) data 
cards, and finally Mach number and angle of attack case cards. 
Note that the present program will not handle a pod and there- 
fore there are no pod panel inputs. However, if the geometry 
input contains a pod description it will be read and ignored. 

Reference length card: This card may be identified as REFL 

in columns 73-80 and contains the following: 


Columns 

1-7 

8-14 

15-21 

22-28 

29-35 

36-42 

43-49 


Variable 
RE FA 


REFB 

REFC 

REFD 

REFL 

REFX 

REFZ 


Description 

Wing reference area. If REFA = 0, 
the reference area is defined by 
the value of REFA in the geometry 
input 

Wing semispan. If REFB = 0, a 
value of 1.0 is used for the 
reference semispan 

Wing reference chord. If REFC = 0 
a value of 1.0 is used for the 
reference chord 

Body (fuselage) reference diameter 
If REFD = 0, a value of 1.0 is 
used for the reference diameter 

Body (fuselage) reference length. 
If REFL = 0, a value of 1.0 is 
used for the reference length 

X coordinate of moment center 

Z coordinate of moment center 
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Wing data cards: The first wing data card is the wing 

leading edge radius card and is required only when K1 = 3. 

This card contains NWAF values of leading edge radius expressed 
in percent chord. It may be identified in columns 73-80 as 
RHOJ where J denotes the number of the last radius given on 
that card. 

Next is the wing panel leading edge card. This card con- 
tains KWAFOR values of wing panel leading edge locations ex- 
pressed in percent chord. This card may be identified in columns 
73-80 as XAFKJ where J denotes the last location in percent chord 
given on that card. Omit if KWAFOR = 0. 

The last wing data card gives the wing panel side edge data. 
This card contains KWAF values of the y ordinate of the panel 
inboard edges. This card may be identified in columns 73-80 as 
YKJ where J denotes the last y ordinate on that card. These 
values are arranged in the order which begins with the most 
inboard panel edge and proceeds outboard. Omit if KWAF = 0. 

Fin (vertical tail) data cards: The first fin data card 

is the fin leading edge radius card and is required only when 
K4 = 3. This card contains NF values of leading edge radius 
expressed in percent chord, one value for each fin. It may be 
identified in columns 73-80 as RHOFIN. 

Next is the fin panel leading edge card for the first fin. 
This card contains KFINOR(l) values of fin panel leading edge 
locations expressed in percent chord. This card may be identi- 
fied in columns 73-80 as XFINKJ where J denotes the fin number. 
Repeat this card for each fin. 

The last fin data card gives the fin panel side edge data 
for the first fin. This card contains KF(1) values of the 
z ordinate of the panel inboard edges. This card may be identi- 
fied in columns 73-80 as ZFINKJ where J denotes the fin number. 
These values are arranged in the order that begins with the most 
inboard panel edge and proceeds outboard. Repeat this card for 
each fin. 

Canard (horizontal tail) data cards: The first canard data 

card is the canard leading edge radius card and is required only 
when K5 = 3. This card contains NCAN values of leading edge 
radius expressed in percent chord, one value for each canard. 

It may be identified in columns 73-80 as RHOCAN. 

Next is the canard panel leading edge card for the first 
canard. This card contains KCANOR(l) values of canard pane i. 
leading edge locations expressed in percent chord. This card 
may be identified in columns 73-80 as XCANKJ where J denotes the 
canard number. Repeat this card for each canard. 
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The last canard data card gives the canard panel side edge 
data for the first canard. This card contains KCAN(l) values 
of the y ordinate of the panel inboard edges. This card may be 
identified in columns 73-80 as YCANKJ where J denotes the canard 
number. These values are arranged in the order that begins with 
the most inboard panel edge and proceeds outboard. Repeat this 
card for each canard. 


Fuselage (body) data cards: The first body card is the 

body meridian angle card. This card contains KRADX(l) values 
of body meridian angle expressed in degrees and may be identi- 
fied in columns 73-80 as PHIKJ where J denotes the body segment 
number. The convention is observed that PHIK = 0 . at the bot- 
tom of the body and PHIK = 180. at the top of the body. Omit 
unless KRADX(l) is negative. Repeat this card for each fuse- 
lage segment. 

The second body card is the body axial station card. This 
card contains KFORX(l) values of the x ordinate of the body 
axial stations and may be identified in columns 73-80 as XFUSKJ 
where J denotes the body segment number. Omit if KFORX(l) = 0. 
Repeat this card for each fuselage segment. 


Mach number and angle of attack card: This card may be 

identified in columns 73-80 as MALPHA and contains the following: 


Columns 


Variable 


Description 


1-7 MACH The subsonic Mach number (includ- 

ing the value MACH = 0 . ) or the 
supersonic Mach niimber at which it 
is desired to calculate the aero- 
dynamic data 

8-14 ALPHA The angle of attack expressed in 

degrees at which it is desired to 
calculate the aerodynamic data 


A series of Mach number and angle of attack combinations 
for the same geometry may be calculated by repeating this card 
with the desired values. 


A value of MACH = -1. on this card signifies the termination 
of the present case. Geometry cards for a new case can follow 
such a terminal card. 
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PROGRAM OUTPUT DATA 


All output is processed by a standard 132 characters-per- 
line printer. The output from each run is always preceded by 
a complete list of the input data cards. The amount and type 
of the remaining output depend on the PRINT option selected, 
the number of panels used, and whether the configuration being 
analyzed is an isolated wing, an isolated body, or a complete 
wing-body-tail combination. The program output options are 
described below: 

PRINT = 0 The program prints the case description, Mach 

number and angle of attack, followed by a table 
listing the panel number, control point coordi- 
nates (both dimensional and non-dimensional) , 
pressure coefficient, normal force, axial force, 
and pitching moment. Separate tables are printed 
for the body and wing panels, noting that any 
tail, fin or canard panels are included with the 
wing output. If the planar boundary condition 
option has been selected, the results for the 
wing upper surface are given in one table, fol- 
lowed by a separate table giving the results for 
the wing lower surface. Additional tables giving 
the total coefficients on the body, the wing and 
the complete configuration follow the pressure 
coefficient tables. These include the reference 
area, reference span and reference chord, the 
normal force, axial force, pitching moment, lift, 
and drag coefficients, and the center of pressure 
of the component. 

PRINT =1 In addition to the output described for PRINT = 0 , 

the program prints out additional tables giving 
the normal force, axial force, pitching moment, 
lift and drag coefficients, and the center of 
pressure of each column of panels on the wing and 
tail surfaces. In addition, the indices of the 
first and last panel in the column are listed, to- 
gether with the span, chord and origin of the 
column . 

PRINT =2 In addition to the output described for PRINT = 1, 

the program prints out tables listing the panel 
number, the source or vortex strength of that 
panel, and the axial velocity u, lateral velocity 
V, and vertical velocity w at the panel control 
point. The normal velocity is also calculated for 
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body panels. Separate tables are printed for 
the body and wing panels, noting again that any 
tail, fin, or canard panels are included with 
the wing output. If the planar boundary con- 
dition option has been selected, separate tables 
are given for the wing upper and lower surfaces. 

PRINT =3 In addition to the output described for PRINT = 2 , 

the program prints out the iteration number, and 
the source and vortex strength arrays obtained 
at each step of the iterative solution procedure. 

PRINT =4 In addition to the output described for PRINT = 3 , 

the program prints out tables of the axial and 
normal velocity components which make up the ele- 
ments of the aerodynamic matrices . The program 
prints out the matrix row number, and gives the 
number of elements in that row. A maximum of 
four matrix partitions will be printed if this 
option is selected, each of which is identified 
by number and its influence description prior to 
printing the velocity component tables. 

If a negative value of PRINT is selected, the program 
prints all the information described above for the positive 
values , together with the complete panel geometry description 
of the configuration following the list of input cards. This 
consists of tables giving the wing panel corner points, control 
points, inclination angles, areas, and chords. If the config— 
u^^tion has a horizontal tail, fin or canard, additional tables 
are printed giving the same information as listed above for the 
wing. Finally, if the configuration includes a body, the body 
panel corner points, control points, areas, and inclination 
angles are listed. 

The program output is illustrated by the sample case pre- 
sented in Appendix III. 
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PROGRAM STRUCTURE 


The program is designed to operate in OVERLAY mode. The 
main overlay progreun is designated USSAERO, and calls the 
three primary overlay programs GEOM, VELCMP, and SOLVE. In 
turn, GEOM calls seven secondary overlay programs CONFIG, 
NEWORD, WNGPAN, NEWRAD, BODPAN, NUTORD, and TALPAN; while 
VELCMP calls three secondary overlay programs BODVEL, LINVEL, 
and WNGVEL. The overlay structure is illustrated on Figure 1. 

The complete program consists of 14 overlay programs and 
19 subroutines. Detailed descriptions of each program and sub- 
routine are given in Appendix I. These descriptions give the 
purpose of the program or subroutine, outline the method used, 
and list the names of the principal variables and constants. 


AUXILIARY FILES 


The program designates TAPE 5 as the input file and TAPE 6 
as the output file. In addition, five auxiliary files are uti- 
lized for temporary storage and data transfer within the pro- 
gram. These files are designated TAPE 7, TAPE 8, TAPE 9, TAPE 
10, and TAPE 11. 

TAPE 7 is used primarily for the storage of the panel geo- 
metry data. The first three records are written by program 
WNGPAN, and contain the wing panel geometry data. If the con- 
figuration has a fin, canard, or horizontal tail, the first 
three records are rewritten by program TALPAN , and subsequently 
contain all wing and tail panel geometry data. The fourth rec- 
ord is written by program BODPAN, and contains the body panel 
geometry data. Additional records are written on this file in 
program VELCMP if the aerodynamic matrix partitions are further 
subdivided into blocks. The elements of the diagonal block 
matrices are stored in individual records on this file, behind 
the panel geometry data. 

TAPE 8 is used exclusively to store the u, v, w velocity 
component arrays, and each record in this file contains one row 
of the velocity component arrays from a given matrix partition. 
In the first partition, NBODY records are written on this file 
by program BODVEL. In the second partition, another NBODY rec- 
ords are written by either program LINVEL or WNGVEL. However, 
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Figure 1 - Program Overlay Structure 




























if the planar boundary condition option with thickness is 
selected, program LINVEL writes an additional NBODY records 
on the file. In the third partition, NWING records are written 
by program BODVEL, and in the fourth partition, another NWING 
records are written by programs LINVEL or WNGVEL. As before, 
LINVEL writes an additional NWING records on this file if the 
planar boundary condition option with thickness is selected. 

Thus, a total of 2 (NBODY + NWING) records are always written 
on TAPE 8, and a maximum of 3 (NBODY + NWING) records if the 
planar boundary condition option with wing thickness is selected. 

TAPE 9 is first used in program CONFIG to store the con- 
figuration geometry data. Five records are written on this 

containing the reference area, wing geometry data, body 
geometry data, fin geometry data, and canard or horizontal tail 
geometry data. Dummy records are written for missing compo- 
nents. TAPE 9 is re-initialized in program VELCMP, and used to 
store the normal velocity arrays. Each record contains one row 
of normal velocities from a given matrix partition. In the 
first partition, NBODY records are written on this file by pro- 
gram BODVEL. In the second partition, another NBODY records 
are written by program LINVEL or WNGVEL. In the third parti- 
tion, NWING records are written by program BODVEL, and in the 
fourth partition, an additional NWING records are written by 
program LINVEL or WNGVEL. Thus, a total of 2 (NBODY + NWING) 
records are written on TAPE 9 . 

TAPE 10 is first used in program NEWRAD as temporary stor— 
age for the body panel corner point coordinates. It is re- 
initialized by program VELCMP, and then used to store the ele- 
ments of the diagonal block matrices, if the matrix partitions 
are further subdivided into blocks. Each record contains one 
row of normal velocities from a given diagonal block matrix 
partition. The records are written at the same time the normal 
velocity arrays for the remainder of the row are written on 
TAPE 9. Thus, a total of 2 (NBODY + NWING) records are also 
written on TAPE 10 . These records are subsequently read by pro- 
gram VELCMP, transferred to TAPE 7, and the file re-initialized 
a second time. TAPE 10 is finally used to store the elements 
of the inverse diagonal block matrices, or the inverse diagonal 
partition matrices, if the matrix is not subdivided into blocks. 
In the former case, the elements of each inverse diagonal block 
matrix are written as a single record on TAPE 10 by subroutine 
DIAGIN , or in the latter case, the elements of each inverse di- 
agonal partition matrix are written on the file by subroutine 
PARTIN. 

TAPE 11 is used to transfer two geometry parameter arrays, 
^^itten as a single record in program VELCMP, to program SOLVE. 

No further use is made of this file. 
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OPERATING INSTRUCTIONS 


The program deck and data deck are loaded in the following 
sequence: job card, system control cards, end-of-record card, 

program deck, end-of-record card, data deck, end-of-file card. 
The data deck is described above in the Program Input Data 
section. 
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APPENDIX I 


PROGRAM AND SUBROUTINE DESCRIPTIONS 
(Arranged in Alphabetical Order) 


This appendix contains a brief outline of the purpose, 
method, and use of each program or subroutine. The 
principal variables and constants in each are listed 
in the order of their first appearance, and identified 
as input or output data. 



PROGRAM BODPAN 


PURPOSE ; 
METHOD : 


USE: 


To revise the axial spacing on the body and 
compute the body panel geometry. 

For each body segment, the x, y, and z co- 
ordinates of the cross sections are read from 
TAPE 10. If the value of KFORX of the segment 
is positive, an array of new axial stations 
for the segment is read in; otherwise the orig- 
inal axial stations are retained. 

The body panel geometry is established by a 
linear interpolation along body meridian lines 
of the y and z coordinates at the new axial 
stations. The interpolation is started with 
the first ring of panels at the nose and con- 
tinued until the last ring of panels on the 
last segment is reached. The corner point 
coordinates, the control point coordinates, 
the inclination angles, and area are calcu- 
lated for each panel in sequence. 

The panel control point coordinates, the panel 
dihedral angle 0 , the panel inclination angle 6 , 
the corner point coordinates and the panel areas 
are stored in the COMMON block POINT, and the 
entire sequence of arrays written as a single 
record on TAPE 7 following the wing and tail 
panel geometry arrays. The remaining body 
geometry parameters are stored in COMMON blocks 
PARAM and BTHET . Finally, if the print option 
is negative, the corner point coordinates, 
control point coordinates, inclination angles, 
and areas are written on the output file. 

CALL OVERLAY (LWB, 1, 5) 

Input : 

PRINT Print option 

NFUS Number of body segments 

KFORX, Number of axial stations on segment 
NFORX, 

JMAX 
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SUBROUTINES 
CALLED : 

ERROR 
RETURNS : 


KRADX, Number of meridian lines on segment 
NRADX, 

KRAD 

XB, Array of original axial stations on 

XFUS segment 

XJ Array of revised axial stations on 

segment 

YB, Arrays of y and z coordinates on 

ZB segment 

Output : 

NBODY Total number of body panels 

NFU Body segment number 

NP Panel number 

IP, Panel identification constants 

IQ 

XC, Arrays of panel corner points 

YC, 

ZC 

XPT, Arrays of panel control points 
YPT, 

ZPT 

THET, Array of panel dihedral angles 
THETA 

DELTA Array of panel inclination angles 

AP, Array of panel areas 

AREA 


PANEL 

The program calls EXIT if NBODY > 600. 



PROGRAM BODVEL 


PURPOSE : 


METHOD : 


To compute the three components of velocity 
induced at specified control points by the 
body panels. 

The X, y, and z coordinates of the control 
point, and the corresponding panel inclination 
angles 6 and 6 are read from COMMON block POINT. 

Starting with the first body segment, the body 
panel corner point coordinates and inclination 
angles are also read from COMMON block POINT 
for each row and column of panels. Considering 
a single body panel, the corner point and con- 
trol point coordinates are transformed to a 
new coordinate system with origin at the first 
corner of the panel and inclined at an angle 0 
'^ith respect to the horizontal. The velocity 
components induced by this inclined constant 
source panel at the given control point are 
computed in the panel coordinate system by sub- 
routine SORPAN , which is called twice to obtain 
the influence of panels located on both right 
and left sides of the body. These velocity 
components are combined and transformed back 
to the reference coordinate system to obtain 
the final u, v, and w components of velocity, 
and the velocity normal to the panel at the 
control point. This process is repeated for 
each panel on the body, following which the u, 

V, and w velocity component arrays are written 
on TAPE 8 , and the array of normal velocities on 
TAPE 9. 

If the control point is in the same ring of 
panels on the body as the influencing panel and 
the body has more than 60 panels , the normal 
velocity at the control point is written on 
TAPE 10 , and its value set to zero in the array 
written on TAPE 9 . This procedure sets up the 
diagonal blocks of the aerodynamic matrix for 
later use in the iterative solution procedure . 

If the print option is selected, the axial and 
normal velocity arrays are written on the output 
tape . 

The process is repeated for each control point. 
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USE: 


CALL OVERLAY (LWB, 2, 1) 


EM, Mach number 

MACH 

NBODY Number of body panels 

LBC Planar boundary condition option 

parameter (logical) 

PRINT Print option parameter 

NPART Matrix partition number 

NMAX Maximum order of diagonal block matrices 

KFUS Number of body segments 

KRADX Number of body panel meridian lines in 
segment 

KFORX Number of body axial stations in segment 

JMAX, Total number of axial stations on body 
MAX 

NPOINT Number of control points 

jrji Array of wing supersonic trailing edge 

indicators 

THET Array of panel inclination angles at 
control point 

THETA Array of body panel inclination angles 

DELTA Array of body panel incidence angles 

DELTI Array of wing panel incidence angles 

XPT, Arrays of wing control point coordinates 

YPT, 

ZPT 

XBT , Arrays of body panel control point 

YBT, coordinates 

ZBT 
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Arrays of body panel corner points 


XC, 

YC, 

ZC 

Output ; 


I 

Control point index 


J 

Body panel index 


ISKIP 

Wing supersonic trailing edge 
indicator 


KF 

Body segment index 


L, 

NC 

Column index 


N 

Row index 


NROW, 

NS 

Number of rows of panels on body 

NCOL 

Number of columns of panels on 

body 

Jl, 
J2 , 
JSl, 
JS2 

Body panel numbers in diagonal 
matrices 

block 

K 

Panel corner index 


SINTI 

sin 0(1) 


COSTI 

cos 9(1) 


XPTI, 

YPTI, 

ZPTI 

Coordinates of control point I 


DI 

tan (S (I) 


DA 

tan 6 (J) 


SINT 

sin 9 (J) 


COST 

cos 0 ( J) 


SINTR 

sin (0(J) - 9(1)) 



SUBROUTINES 
CALLED : 

ERROR 
RETURNS : 


SINTL 

sin (9(J) + 9(1)) 

COSTR 

cos (9 ( J) - 9(1)) 

COSTL 

cos (9 (J) + 9(D) 

XCOR, 

YCOR, 

ZCOR 

Coordinates of panel corner points 
in panel coordinate system 

CX 

Panel chord length 

XI, 

YI, 

ZI 

Coordinates of control point I in 
panel coordinate system 

XJ, 

YJ, 

ZJ 

Coordinates of body panel J control 
point in panel coordinate system 

UR, 

VR, 

WR 

Velocity components at control point I 
induced by body panel J on right side 
of body, in body panel coordinate system 

UL, 

VL, 

WL 

Velocity components at control point I 
induced by body panel J on left side of 
body, in body panel coordinate system 

UB, 

VB, 

WB 

Arrays of velocity components at control 
point I in reference coordinate system 

VI, 

WI 

Arrays of velocity components at control 
point I in control point panel co- 
ordinate system 

AN 

Array of velocities normal to control 
point panel I 

DN 

Array of normal velocities in diagonal 
block matrices 


SORPAN 

None 
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SUBROUTINE COMCU 


PURPOSE : 


METHOD : 


To fit a composite cubic through n points 
(^i' Yi) i.e., a separate cubic between each 
pair of adjacent points, such that the n-1 
cubics are so determined that each matches its 
neighbors in function value and in the first 
two derivatives . 

Rather than solve simultaneously for the 4 (n-1) 
cubic coefficients, the approach here is to 
solve simultaneously for the slopes of the 
composite cubic at the given n points. Thus a 
linear system of order n, rather than 4n-4 is 
involved. It can be shown that a necessary 
and sufficient condition for continuity of the 
second derivative is that 




(Xi-Xi.^) (x^^.,-x^) 


(x. -X. ) ^ (y . ^ -y • ) 

1 1-1 “'i+i -‘i 




for i = 2, 3, . . ., n-1 

This yields n-2 equations in the n unknowns# 
yl,i = l, 2, . . .,n. For the 1®^ and n"^" 

equations of the linear system, the boundary 
conditions on yj and yn are used. This has 
been generalized to permit any combination of 
a given y' or y' ' at the end points, e.g., 
yJ and y^' can be given as the boundary 
conditions. The second derivative of a cubic 
through two points can be expressed as a 
function of the first derivatives and of the 
given point coordinates as follows : 


X 



y" 


X2-X1 



and 


X -X 

n n - 1 


y.. = .3 ^ ^ 

n Xj^““Xj^_ 2 ri"" 1 n 
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Whether the boundary conditions involve first 
or second derivatives (or both) and no matter 
what the spacing of the so long as the 
form a strictly monotone sequence, the coeffi- 
cient matrix of the linear system is tridiagonal 
(all elements are zero except on the principal 
diagonal, the first subdiagonal, and the first 
superdiagonal). When n is large, a considerable 
time saving and an enormous storage saving can 
result if the special structure of this matrix 
is taken advantage of. Hence, this subroutine 
stores the matrix elements in 4n locations (as 
opposed to n^) and then solves the system. 

The actual coefficients of the n-1 cubics of the 
composite cubic are not found by COMCU. Since 
on any subinterval x , x^+j, a cubic is uniquely 
determined by the known two points and two 
slopes, the calling program can find the four 
coefficients of each cubic independently and 
may often need to do so for only one of the n-1 
cubics. In any case, the subroutine CUBIC2 
specifically finds a cubic, given two points 
and the slope at each point. 

USE: CALL COMCU (DA, DB, S, X, Y, L, M, N, NDA,NDB) 

Input : 

X Array of x-abscissae of input points 

Y Array of y-ordinates of input points 

N Number of input points 

NDA Order (1 or 2) of derivative at X(l) 

NDB Order (1 or 2) of derivative at X(N) 

DA Value of derivative at X(l) 

DB Value of derivative at X(N) 

L Code 


= 1, 

if 

single 

precision is 

to 

be 

used 

= 2, 

if 

double 

precision is 

to 

be 

used 
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Output ; 


SUBROUTINES 
CALLED : 

ERROR 
RETURNS : 

RESTRICTIONS : 


S Array of first derivatives 

M Error return 

= 0 - success 

^ 0 - error detected 

None 

If overflow occurred, M = 1. Otherwise, M = 0. 

The x-abscissae must form a strictly monotone 
sequence. N < 400. 
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PROGRAM CONFIG 


PURPOSE : 


METHOD : 


USE: 


To input the geometrical description of the con- 
figuration using the same input data as program 
START of reference 1. 

The configuration reference area is read from the 
input file if JO ^ 0, otherwise the reference area 
is set equal to unity. The reference area is then 
written on TAPE 9. If Jl ^ 0, the wing geometry 
data is read from the input file in the order spec- 
ified in reference 1. The program computes the 
upper and lower surface coordinates of the wing 
airfoils, and writes the entire wing geometry array 
as one record on TAPE 9. 

If J2 0, the body geometry data is also read from 
the input file in the order specified in reference 
1 for each body segment. For arbitrary cross- 
sections, the y and z ordinates of the body seg- 
ment are read in, but for circular cross-sections, 
the body cross-sectional area is read in and the 
corresponding radius calculated by the program. 

The entire body geometry array is then written as 
one record on TAPE 9 . 

If J3 ^ 0, the pod geometry is read in, but no 
further use is made of this data. 

If J4 0 , the fin geometry data is read in. The 
program computes the coordinates of the fin airfoils 
and writes the fin geometry array as one record on 
TAPE 9. Similarly, if J5 ^ 0, the tail or canard 
geometry data is read in, the tail airfoil coordi- 
nates calculated, and the tail geometry array writ- 
ten on TAPE 9 . 

If one or more of the above components is missing, 
the program writes a dummy record on TAPE 9 and 
continues . 

CALL OVERLAY (LWB, 1, 1) 

Input : 

JO Reference area parameter 

Jl Wing definition parameter 
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J2 

J3 

J4 

J5 

J6 

REFA 

ABCD 

NWAF 

NWAFOR 

WAFORG 

XAF 

WAFORD 

TZORD 

NFUS 

NRADX, 

NRAD 

NFORX, 

NFUSOR 

XFUS 

ZFUS 

SFUS 

FUSARD 

NP 

NPODOR 


Body definition parameter 

Pod definition parameter 

Fin definition parameter 

Canard or tail definition parameter 

Body camber parameter 

Reference area 

Dummy array 

Number of wing airfoil sections 

Number of ordinates used to define each 
wing airfoil section 

Origin and chord length of each wing 
airfoil (x, y, z, c) 

Array of percent chords for wing 
airfoil ordinates 

Array of half-thickness ordinates in 
percent chord 

Array of mean camber line ordinates 

Number of body segments 

Number of points used to define half- 
section of body segment 

Number of axial stations on body 
segment 

Array of axial stations on body segment 

Array of body camber ordinates 

Array of y and z ordinates used to 
define half-section of arbitrary body 
segment 

Array of body cross-sectional areas 
Number of pods 

Number of axial stations on pod 


NF 

NFINOR 

FINORG 

XFIN 

FINORD 

NCAN 

NCANOR 

CANORG 

XCAN 

CANORD 
Output ; 
REFA 
WAFOR 

TZORD 

WAFORD 

J2TEST 

FUSRAD 

FINOR 


Number of fins 

Number of ordinates used to define fin 
airfoil sections 

Origin and chord lengths of fin 
airfoils (x, y, z, c) 

Array of percent chords for fin air- 
foils 

Array of fin airfoil half-thickness 
ordinates in percent chord 

Number of tails or canards 

Number of ordinates used to define tail 
or canard airfoil sections 

Origin and chord lengths of tail or 
canard airfoils (x, y, z, c) 

Array of percent chords for tail or 
canard airfoils 

Array of airfoil half-thickness ordinates 


Reference area 

Array of wing half-thickness ordinates 
(percent chord) 

Array of wing camber line ordinates 

Array of x, y, z, coordinates defining 
upper and lower surfaces of wing (not 
used) 

Parameter to specify body camber and 
cross-section definition 

Array of body radii (circular cross- 
sections only) 

Array of fin half-thickness ordinates 
(percent chord) 
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FINCR 


Array of fin camber line ordinates (set 
zero) 


FINORD, Arrays of x, y, z coordinates defining 
FINX2 , left and right surfaces of fin (not 
FINX3 used) 

CANOR Array of tail or canard half-thickness 

ordinates 

CANCR Array of tail or canard camber line 

ordinates 

CANORX, Arrays of x, y, z coordinates defining 
CANORD, upper and lower surfaces of tail or 
CANORl canard (not used) 

BLOCK Dummy array used for storing geometry 

data on TAPE 9 

SUBROUTINES 
CALLED: None 

ERROR 

RETURNS : None 
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SUBROUTINE CUBIC2 


PURPOSE : 


METHOD : 


USE: 


SUBROUTINES 
CALLED : 

ERROR 
RETURNS : 


To fit a ciibic to two points, being given the 
slope at each. 

The subroutine sets up the system of four simul- 
taneous equations expressing the four given 
conditions and solves it for the coefficients 
of the cubic. 

CALL CUBIC2 (X, Y, D, C, M) 

Input : 


X 

Array 

of 

x-coordinates 

Y 

Array 

of 

y-coordinates 

D 

Array 

of 

first derivatives 

Output : 



C 

Array 

of 

cubic coefficients 


Error 

return 

M ' 

= 1 - 

success 


^ 1 - 

error detected 


None 


If M = 2, overflow occurred. If M = 3, 

X(l) = X(2). Otherwise, call is successful, 
and M = 1 


RESTRICTIONS 


X(l) X(2) 



SUBROUTINE DERIV 


PURPOSE : 

METHOD : 
USE: 


To fit a chain of cubic curves through a set of 
N points (x^, y.) having continuous first and 
second derivatives at the intermediate points 
and specified first or second derivative at the 
end points . 

The method outlined in subroutine SCAMP4 is 
applied. 

CALL DERIV (X, Y, N, NDA, DA, FD) 


Input 


Array of x values 
Array of y values 
Number of points 

The order of the derivative to be 
specified at the first point 

The value of the derivative to be 
specified at the first point 


Output : 


FD Array of first derivatives at the points 

SUBROUTINES 
CALLED: SCAMP 4 


ERROR 

RETURNS : None 
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FUNCTION DERIVl 


PURPOSE : 


METHOD : 


USE: 


SUBROUTINES 
CALLED : 

ERROR 

RETURNS: 


To find the first derivative of the quadratic 
through three given points at a specified one 
of these points . This provides a good approx- 
imation to the slope of a function at a point, 
particularly if the other two points used are 
nearby . 

The subroutine simply finds the unique poly- 
nomial of degree two through the given points 
then evaluates its first derivative at the 
specified point. 

D = DERIVl (X, Y, N) 


Input : 


X Array of x-coordinates 


Y 

N 


Array of y-coordinates 
Code 

= 1, 2 , or 3 indicating point at which 
derivative is desired 


None 

None 


RESTRICTIONS The x-coordinates must be distinct, but need 
not be in any order or evenly spaced 



FUNCTION DERIV2 


PURPOSE : 
METHOD : 

USE; 


SUBROUTINES 

CALLED; 

ERROR 
RETURNS ; 

RESTRICTIONS ; 


To find the second derivative of the cvibic 
through four given points (x- , y-) at an 
arbitrary point whose x coordinates if given. 

The subprogram simply finds the unique poly- 
nomial of degree three through the given points , 
then evaluates its second derivative at the 
desired x, which need not be one of the four 
given x^^ . 

D = DERIV2 (X, Y, XX) 

Input ; 

X Array of x coordinates 

Y Array of y coordinates 

XX X coordinate of point at which second 

derivative is desired 


None 


None 

The X coordinates must be distinct but can be 
in any order and unevenly spaced. 


40 



SUBROUTINE DIAGIN 


PURPOSE : 
METHOD ; 


USE: 


To invert the diagonal blocks of the matrix 
and store the results on TAPE 10 . 

If the order of the body matrix partition 
exceeds 60, the diagonal blocks of the body 
matrix are read from TAPE 7, the block matri- 
ces inverted, and the inverse matrices stored 
on TAPE 10. Otherwise, the complete body 
matrix partition is read from TAPE 9 , the 
matrix inverted, and the inverse stored on 
TAPE 10. 

A similar procedure is followed for the wing 
matrix partition. 

CALL DIAGIN 


Input : 


NWING 

NBODY 

NMAX 

NDIM 

NBBLOK 

NWBLOK 

NBROW 

NWROW 

D 


Number of wing panels 

Number of body panels 

Maximum order of diagonal block 
matrices (60) 

Matrix dimension statement size 

Number of diagonal blocks in body 
matrix partition 

Number of diagonal blocks in wing 
matrix partition 

Order of diagonal blocks in body 
matrix partition 

Order of diagonal blocks in wing 
matrix partition 

Array of matrix elements 


Output : 

NB Body diagonal block number 



NW Wing diagonal block mimber 

NROW Number of rows in diagonal block 
NCOL Number of columns in diagonal block 
D Array of inverse matrix elements 


SUBROUTINES 


CALLED : 

None 

ERROR 


RETURNS : 

None 
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SUBROUTINE FORMOM 


PURPOSE : 
METHOD : 


To calculate the force and moment coefficients 
on body, wing and tail components. 

Depending on the component being analyzed and 
the boundary condition option selected, exe- 
cution of this subroutine follows one of three 
paths. In all three paths, the panel inclina- 
tion angles, control point coordinates, areas 
and chords are obtained from COMMON block POINT . 
The pressure coefficients are obtained from 
COMMON block SCRAT. 

If the component is a body, the normal force, 
axial force, and moment about the origin of 
coordinates is computed for each panel and the 
results summed. The total axial force, normal 
force and pitching moment of the body are stored 
in COMMON block FORM. The body panel force and 
moment arrays are written on the output file, 
together with the panel number, control point 
coordinates, and pressure coefficient. The con- 
trol point coordinates are non-dimensionalized 
by dividing by the body reference length or 
diameter, and both dimensional and non-dimensional 
coordinates are presented in this array . Finally , 
the normal and axial force coefficients, the 
pitching moment coefficient about the origin of 
coordinates, the lift and drag coefficients, 
and the center of pressure of the body expressed 
as a fraction of the reference chord are computed 
and written on the output file. 

If the component is a wing or tail, and the non- 
planar boundary condition option has been select- 
ed, the subroutine first calculates the chord 
length of each column of panels, and assigns the 
chord length and the axial coordinate of the 
leading edge to each panel in the column. The 
normal force , axial force and moment about the 
origin of coordinates is computed for each panel, 
and the results summed. The wing panel force 
and moment arrays are then written on the output 
file, together with the panel nxamber , control 
point coordinates, and pressure coefficient. The 
control point coordinates are non-dimensionalized 
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by dividing them by the chord length or 
reference semi-span, and both dimensional and 
non-dimensional coordinates are presented in 
the array . Next , the normal and axial force 
coefficients, the pitching moment coefficient 
about the origin of coordinates, the lift and 
drag coefficients and the center of pressure 
of the wing, expressed as a fraction of the 
reference chord, are computed and written on 
the output file. If a non-zero print option 
has been selected, the siobroutine proceeds to 
calculate the spanwise load distribution on 
the wing and tail. The forces and moment acting 
on each panel in a given column are summed, and 
the force and moment coefficients and center of 
pressure of the column computed and written on 
the output file. The summation includes panels 
on both the upper and lower surfaces of the 
column. 

If the component is a wing or tail and the 
planar boundary condition option has been 
selected, the subroutine is called twice. On 
the first pass, the forces and moment acting 
on the upper surface is calculated, followed 
by the forces and moment on the lower surface 
in the second pass. In this case, the sub- 
routine performs an interpolation to determine 
the wing or tail panel slope, pressure coeffi- 
cient, and control point at the panel centroid 
prior to the panel force and moment calculations. 
The remainder of the calculation procedure is 
similar to that described above for the wing 
with non-planar boundary condition option. The 
axial force, normal force, and pitching moment 
of the wing upper surface are stored in COMMON 
block FORM on the first pass, and added to 
those calculated on the lower surface in the 
second pass to obtain the total forces and 
moments acting on the wing. This information 
is then used to calculate the total force and 
moment coefficients, lift, drag, and center of 
pressure of the wing and tail, and the results 
stored on the output file. The spanwise load 
distribution on the wing and tail is also calcu- 
lated if the print option is other than zero, 
using the procedure described above. 


USE: 


CALL FORMOM (NPAN, NPASS , ALFA, COMPT) 


Input : 
NPAN 

NPASS 

COMPT 

ALFA 

ALPHA 

PRINT 

LBC 

MACH 

NBODY 

NWING 

ARRAY 

CHORD 

DZTDX 

XC 

TITLEl, 

TITLE2 

DELTA 

THET 

NSEG 


Number of panels on the component 
being analyzed 

Pass number 

Component identification integer 

COMPT = 1 Body component 

COMPT = 2 Wing or tail component 

Angle of attack (radians) 

Angle of attack (degrees) 

Print option (integer) 

Planar boundary condition option 
parameter (logical) 

Mach number 

Number of body panels 

Number of wing and tail panels 

Array of panel geometrical parameters 

Array of wing panel chord lengths 

Array of wing panel half-thickness 
slopes 

Array of panel corner point coordinates 
Case identification arrays 

Array of panel incidence angles , or 
panel camber slopes on wing with planar 
boundary condition option 

Array of panel inclination angles 

Number of wing and tail segments 
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NCOL 

NROW 

SINS, 
COSS 

XCPT 

LOCPT 

CP 

XPT, 

YPT, 

ZPT 

XLEW 

SPNW 

AREA 

RE FA 

REFB 

REFC 

REFD 

REFL 

REFX, 

REFZ 

Output ; 

XON 

SGN 

SIND, 

COSD 

SINT, 
COST 


Number of columns of panels in segment 

Number of rows of panels in segment 

Array of sines and cosines of segment 
inclination angle 

Array of chord fractions for control 
point location 

Array of control point location indicators 
Array of panel pressure coefficients 
Arrays of control point coordinates 

Array of wing colximn leading edge origins 

Array of wing column spanwise widths 

Array of panel areas 

Reference wing area 

Reference wing span 

Reference wing chord 

Reference body diameter 

Reference body length 

Coordinates of moment reference point 


X coordinate of body nose 

Wing upper and lower surface sign 
parameter array 

Trigonometric function arrays of panel 
incidence angle 6 

Trigonometric function arrays of panel 
inclination angle 9 



SIAL, 

COAL 

NC 

NR 

NRl 

TAND 

XS, 

PT 

RL 

RT 

CHD 

XLE 

NP 

IP 

CN 

CT 

CM 

XP, 

yp, 

ZP 

FI, 

F2, 

F3 

FAK 

DCN 


Trigonometric functions of angle of 
attack a 

Number of columns of panels in segment 
Number of rows of panels in segment 
NR + 1 

Slope of wing upper or lower surface 
at centroid, planar boundary condition 
option 

Control point location in chord 
fractions 

Chord fraction of control point from 
leading edge 

Chord fraction of control point from 
trailing edge 

Array of wing column chord lengths 

Array of wing column leading edge 
origins 

Number of panels on component 
Panel number 

Normal force coefficient 
Axial force coefficient 
Pitching moment coefficient 
Panel control point coordinates 

Direction cosines of panel normal vector 

Panel area factor 
Panel normal force 
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DCT 

Panel axial force 


DCM 

Panel pitching moment 


XQ, 

Non-dimensional panel control point 


YQ, 

coordinates 


ZQ 



CL 

Lift coefficient 


CD 

Drag coefficient 


DXN 

Distance of center of pressure from 
origin, in reference chord lengths 


11 

Number of first panel in wing column 


12 

Number of last panel in wing colxomn 


IZ 

Index counter 


DELY 

Wing coliimn spanwise width 


XL 

Wing column leading edge origin 


CNS 

Normal force on wing upper surface 


CTS 

Axial force on wing upper surface 


CMS 

Pitching moment on wing upper surface 

SUBROUTINES 
CALLED : 

None 


ERROR 
RETURNS : 

None 
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PROGRAM GEOM 


PURPOSE: To input configuration geometry and specify panel 

subdivision of the components. A complete descrip- 
tion of the input geometry cards is given in the 
program listing. 

METHOD: The case identification and initial configuration 

parameters are read from the input file. The 
secondary overlay program CONFIG is then called 
to complete the input of the configuration descrip- 
tion, The auxiliary case identification is then 
read, followed by the boundary condition and print 
option. Finally, the revised configuration para- 
meters used for specifying the panel subdivisions 
are read. Depending on the values of the revised 
configuration parameters, the program calls the 
secondary overlay programs NEWORD, WNGPAN , NEWRAD, 
BODPAN, NUTORD or TALPAN, which interpolate the 
input geometry and calculate the corner points , 
control points and inclination angles of the panels 
on the wing, body, or tail. 

USE: CALL OVERLAY (LWB, 1, 0) 



JO , KO 
Jl, Kl 
J2, K2 
J3, K3 
J4, K4 
J5, K5 
J6, K6 
NWAF 


NWAFOR 


Case description. 

Auxiliary case description. 

Reference area and length parameters. 

Wing definition parameters. 

Body definition parameters. 

Pod definition parameters. 

Fin definition parameters. 

Canard or tail definition parameters. 

Body camber parameters . 

Number of wing airfoil sections. 

Number of ordinates used to define each 
wing airfoil section. 
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NFUS, 

KFUS 

NRADX 

NFORX 

NP 

NPODOR 

NF 

NFINOR 

NCAN 

NCANOR 

KWAF 

KWAFOR 

KRADX 

KFORX 

KF 

KFINOR 

KAN 


Number of body segments 

Number of points used to define half- 
section of body segment 

Number of axial stations on body segment 
Number of pods 

Number of axial stations on pod 
Number of fins 

Number of ordinates used to define each 
fin airfoil section 

Number of tails (canards) 

Number of ordinates used to define each 
tail airfoil section 

Number of streamwise panel edges on wing 

Number of ordinates used to define the 
leading and trailing edges of the wing 
panels 

Number of meridian lines used to define 
panel edges on body segment 

Number of axial stations used to define 
leading and trailing edges of panels on 
body segment 

Number of sections used to define the 
streamwise panel edges on fin 

Number of ordinates used to define the 
leading and trailing edges of the panels 
on fin 

Number of ordinates used to define the 
leading and trailing edges of the panels 
on tail or canard 


REFA, 


REFAR 

Wing 

REFB 

Wing 

REFC 

Wing 


reference area 
reference span 
reference chord 



PROGRAMS 
CALLED : 


REFD Body reference diameter 

REEL Body reference length 

REFX, Coordinates of moment reference point 

REFZ 

LINBC Boundary condition selection parameter 
(integer) 

THICK Wing thickness selection parameter 

(integer) 

PRINT Output print selection parameter 

(integer) 


Number of body panels 

Number of wing and tail panels 

Number of tail panels (not used) 

Number of control points on wing and tail 
Boundary condition parameter (logical) 
Wing thickness parameter (logical) 

Tail parameter (logical) 


KOL 

Number 

tail 

Of 

columns 

of panels on wing and 

NSEG 

Number 

of 

wing and tail segments 

BLOCK 

Dummy array used 
data on TAPE 9 

for storing geometry 

OVERLAY 

(LWB, 

1, 

1) 

(CONFIG) 

OVERLAY 

(LWB, 

1/ 

2) 

(NEWORD) 

OVERLAY 

UWB, 

1, 

3) 

(WNGPAN) 

OVERLAY 

(LWB, 

1, 

4) 

(NEWRAD) 

OVERLAY 

(LWB, 

1, 

5) 

(BODPAN) 

OVERLAY 

(LWB, 

1, 

6) 

(NUTORD) 

OVERLAY 

(LWB, 

1, 

7) 

(TALPAN) 
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ERROR 
RETURNS : 


The program is terminated if: 

(a) An end of file is read on TAPE 5 

(B) KOL > 20 

(c) KRADX(l) > 21 
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SUBROUTINE ITRATE 


PURPOSE : 


METHOD : 


To solve the boundary condition equations by an 
iterative procedure and determine the strengths 
of the body sources, and wing and tail vortices. 

The first approximation to the body panel source 
strengths is obtained by post-multiplying the 
inverted body diagonal block matrices written 
on TAPE 10 by the body normal velocity array. 

The first approximation to the wing and tail 
panel vortex strengths is obtained in a similar 
manner. If the absolute value of the print 
option is greater than two, the approximate 
source and vortex strengths are written on the 
output file. 

The body normal velocity array is then revised 
by subtracting an incremental normal velocity 
array from the original normal velocity array. 

The incremental values are obtained in two steps. 
In the first step, the matrix giving the influ- 
ence of the body sources on the body panel con- 
trol points is read from TAPE 9 and multiplied 
by the approximate body source strengths . In 
the second step, the matrix giving the influence 
of the wing and tail vortices on the body panel 
control points is read from TAPE 9 and post- 
multiplied by the approximate wing and tail 
vortex strengths. The incremental normal veloc- 
ity array on the body control points is the sum 
of these two contributions. 

The wing and tail normal velocity array is 
revised in a similar manner. The revised normal 
velocity arrays are then used to obtain a second 
approximation to the source and vortex strengths 
by repeating the above procedure. This iteration 
procedure is repeated until the maximum number 
of iterations has been completed. 

If the order of the body partition does not exceed 
60 , the same procedure is followed except that 
the first step in the determination of the incre- 
mental normal velocities on the body is omitted. 

If the order of the wing partition does not exceed 
60, the same procedure is followed except that 
the second step in the determination of the 
incremental normal velocities on the wing is 
omitted. 
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USE : 


CALL ITRATE 


Input : 

IMAX 

NBODY 

NWING 

NBBLOK 

NWBLOK 

NBROW 

NWROW 

NB 

NW 

D 

A 

PRINT 
Output ; 
RB 

RW 

NBLOK 

NROW 

NCOL 

GB 

GW 


Maximum number of iterations 
Number of body panels 
Number of wing panels 

Number of diagonal blocks in body matrix 
partition 

Number of diagonal blocks in wing matrix 
partition 

Number of rows in diagonal blocks in 
body matrix partition 

Number of rows in diagonal blocks in 
wing matrix partition 

Array of body normal velocities 

Array of wing normal velocities 

Array of diagonal block matrix elements 

Array of normal velocity matrix elements 

Print option parameter 

Array of revised normal velocities on 
body 

Array of revised normal velocities on 
wing 

Number of diagonal block matrices 

Number of rows in diagonal block matrix 

Number of columns in diagonal block 
matrix 

Array of body source strengths 
Array of wing vortex strengths 
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DNB 

Incremental normal 

velocities 

on body 


DNW 

Incremental normal 

velocities 

on wing 


TIME 

Elapsed time 



SUBROUTINES 
CALLED : 

SECOND 




ERROR 
RETURNS : 

None 






SUBROUTINE INVERT 


PURPOSE : 
METHOD : 


USE: 


Matrix inversion subroutine 

Subroutine INVERT is a simple matrix inversion 
procedure based on Gauss- Jordan elimination 
without pivoting. 

CALL INVERT (A, lA, NMAX) 

Input I 

A Name of matrix to be inverted 

lA Number of rows and columns in matrix A 

NMAX Maximum dimensions specified for A in 
calling subroutine 

Output : 

A Inverse of A 


SUBROUTINES 
CALLED : None 

RESTRICTIONS: NMAX not greater than 115 

ERROR 

RETURNS: Subroutine calls EXIT if matrix is singular 
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PROGRAM LINVEL 


PURPOSE : 


METHOD : 


To calculate the three components of velocity 
induced at specified control points by source 
and vortex distributions on panels located in 
the plane of the wing or tail surfaces. 

The X, y, and z coordinates of the control 
point, and the corresponding panel inclination 
angles 9 and 6 are read from COMMON block POINT . 

Starting with the first wing segment, the panel 
leading and trailing edge slopes are calculated 
and stored. The program then computes the veloc- 
ity components induced by the panel corner points 
along the inboard edge of the first column of 
panels. These calculations are performed by sub- 
routines VORVEL and SORVEL, which return the 
three components of velocity induced by constant 
and linearly varying vortex and source distribu- 
tions. These subroutines are called twice to 
obtain the contributions of both left and right 
wing panels. In addition, a second call to sub- 
routine VORVEL is required at panel trailing 
edge corner points if the panel spacing is not 
uniform. 

To compute the velocity components induced by 
the panel corner points along the outboard edge 
of this and the remaining colvimns of panels, 
the procedure is repeated. However, for the 
remaining columns of panels , advantage is taken 
of the fact that the velocity components along 
the inboard edges of a given column of panels 
are the same as those computed at the outboard 
edges of the previous column of panels. There- 
fore, the inboard velocity components are not 
recomputed, but stored in temporary arrays prior 
to the calculation of the outboard velocity 
component arrays . 

Once the velocity components induced by the panel 
corner points along the outboard edge of a given 
column of panels are computed, the inboard and 
outboard influences of each panel in the column 
are combined to obtain the resultant panel veloc- 
ity components. First the velocity components 
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USE: 


induced by the right and left wing panels are 
calculated, using appropriate combination rules 
for the source and vortex panels, and applying 
special rules for leading and trailing edge 
panels. Finally the contributions of the left 
and right wing panels are combined, the velocity 
components transformed back to the reference 
coordinate system, and the velocity normal to 
the panel at the control point computed. 

The procedure is repeated for each column of 
panels in each wing segment, until all wing 
panels are accounted for. At this point the u, 
V, and w components of velocity induced by the 
source panels are written as a single record on 
TAPE 8, followed by the u, v, and w components 
of velocity induced by the vortex panels. If 
the thickness option is not requested, only the 
vortex panel arrays are written on this tape. 

The normal velocities are then written as a 
single record on TAPE 9. If the control point 
is in the same column of panels on the wing as 
the influencing panel, and the wing has more 
than 60 panels, the normal velocity at the con- 
trol point is written on TAPE 10 and its value 
set to zero in the array written on TAPE 9 . 

This procedure sets up the diagonal blocks of 
the aerodynamic matrix for later use in the 
iterative solution procedure. Finally, if the 
print option is selected, the axial and normal 
velocity component arrays induced by the vortex 
panels and source panels are written on the 
output tape. 

The process is repeated for each control point. 
CALL OVERLAY (LWB, 2, 2) 


Input : 


Note : 

The word wing includes any tail, fin 
canard in the following descriptions 

MACH 

Mach number 

PRINT 

Print option parameter 

THK 

Wing thickness option parameter 
(logical) 
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NPART 

Matrix partition number 

NMAX 

Maximum order of diagonal block 
matrices 

NWING 

Number of wing panels 

NSEG 

Number of wing segments 

NROW 

Number of rows of panels in segment 

NCOL 

Number of columns of panels in 
segment 

NWT 

Tail segment identification parameter 

NPOINT 

Number of control points 

IT 

Array of wing supersonic trailing 
edge indicators 

XPT, 

YPT, 

ZPT 

Arrays of control point coordinates 

THET 

Array of panel inclination angles at 
control points 

DELTA 

Array of panel incidence angles at 
control points 

XC, 

YC, 

ZC 

Arrays of wing panel corner point 
coordinates 

GOSS 

cos 9 ( J) 

SINS 

sin 9 (J) 

CHORD 

Array of wing panel chords 

Output : 


I, II 

Control point index 

J, 

JSAVE 

Wing panel index (vortex panels) 

K, 

KSAVE 

Wing panel index (source panels) 


L 


Panel row index 



M 

N 

NP, 

NPSAVE 

NWTHK 

BETA 

SUB 

SUPTE 

CON, 

BCON 

CONT, 

BCONT 

ISKIP 

NR, NS 

NC 

NT 

Jl, J2 , 

JSl, 

JS2 

SINTI 

COSTI 

XI, YI, 
ZI 

DI 

BLE ,BL 

BTE 

FLAG 


Panel column index 
Wing segment index 
Panel number 

Number of wing panel source distributions 

Mach number parameter 

Subsonic flow parameter (logical) 

Supersonic trailing edge parameter 
(logical) 

Constants for vortex panel velocity 
components 

Constants for source panel velocity 
components 

Wing supersonic trailing edge indicator 

Number of rows of panels in segment 

Number of columns of panels in segment 

Tail segment identification parameter 

Wing panel numbers of diagonal block 
matrices 

sin 0 (I) 
cos 0(1) 

Coordinates of control point I 
tan 6(1) 

Array of panel leading edge slopes 

Panel trailing edge slope 

Logical variable denoting presence of 
additional column of vortex panels 
extending from the center line to the 
inboard edge of the wing 


EPOS 


Panel leading edge slope sign 
parameter (logical) 

COST cos 0(J) 

SINT sin e(J) 

SINTR sin (9 (J) - 0 (I) ) 

COSTR cos (9(J) - 0(D) 

SINTL sin (0 (J) + 0(1) ) 

COSTL cos (0(J) + 0(D) 

XC, Arrays of panel corner point 

YC, coordinates 

ZC 

DX, Control point coordinates in panel 

DY, reference system 
DZ 

AL, Difference between panel leading 

AB, and trailing edge slopes 

AT 

CL, Panel chord length along edge 

CT, 

CC 

ABA Absolute value of (AL - AT) 

ACL Absolute value of (CL - CT) 

ML Panel edge indicator 

AMP Reciprocal of panel chord 

X Dummy variable 

UCOR, Velocity components induced by outboard 
VCOR, corners of right wing panels containing 
WCOR constant vortex distributions 

ULOR, Velocity components induced by outboard 
VLOR, leading edge corner of right wing panels 
WLOR containing linearly varying vortex 

distributions 
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UTOR, Velocity components induced by outboard 
VTOR, trailing edge corner of right wing 
WTOR panels containing linearly varying 

vortex distributions 

RCOR, Velocity components induced by outboard 

SCOR, corners of right wing panels containing 

TCOR constant source distributions 

RLOR, Velocity components induced by outboard 

SLOR, corners of right wing panels containing 

TLOR linearly varying source distributions 

UCOL, Same as UCOR, VCOR, WCOR for outboard 

VCOL, corners of left wing panels 

WCOL 

ULOL, Same as ULOR, VLOR, WLOR for outboard 

VLOL, corners of left wing panels 

WLOL 

UTOL, Same as UTOR, VTOR, WTOR for outboard 

VTOL, corners of left wing panels 

WTOL 

RCOL, Same as RCOR, SCOR, TCOR for outboard 

SCOL, corners of left wing panels 

TOOL 

RLOL, Same as RLOR, SLOR, TLOR for outboard 

SLOL, corners of left wing panels 

TLOL 

UCIR, Same as UCOR, VCOR, WCOR for inboard 

VCIR, corners of right wing panels 

WCIR 

ULIR, Same as ULOR, VLOR, WLOR for inboard 

VLIR, corners of right wing panels 

WLIR 

UTIR, Same as UTOR, VTOR, WTOR for inboard 

VTIR, corners of right wing panels 

WTIR 

RCIR, Same as RCOR, SCOR, TCOR for inboard 

SCIR, corners of right wing panels 

TCIR 
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RLIR, Same as RLOR, SLOR, TLOR for inboard 

SLIR, corners of right wing panels 

TLIR 

UCIL, Same as UCOR, VCOR, WCOR for inboard 

VCIL, corners of left wing panels 

WCIL 

ULIL, Same as ULOR, VLOR, WLOR for inboard 

VLIL, corners of left wing panels 

WLIL 

UTIL, Same as UTOR, VTOR, WTOR for inboard 

VTIL, corners of left wing panels 

WTIL 

RCIL, Same as RCOR, SCOR, TCOR for inboard 

SCIL, corners of left wing panels 

TOIL 

RLIL, Same as RLOR, SLOR, TLOR for inboard 

SLIL, corners of left wing panels 

TLIL 

ULR, Velocity components induced by right 
RCR, wing panels containing linearly varying 

VLR, vortex distributions with zero strength 

SCR, along leading edge 

WLR, 

TCR 

ULL, Same as above for left wing panels 
RCL, 

VLL, 

SCL, 

WLL, 

TCL 

UCR, Velocity components induced by right 
VCR, wing panels containing linearly varying 

WCR vortex distributions with zero strength 

along trailing edge 

UCL, Same as above for left wing panels 
VCL, 

WCL 
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RLR, Velocity components induced by right 
RTR, wing panels containing linearly varying 

SLR, source distributions with zero strength 

STR, along leading edge 

TLR, 

TTR 

RLL, Same as above for left wing panels 
RTL, 

SLL, 

STL, 

TLL, 

TTL 

UTR, Velocity components induced by right 
VTR, wing panels containing linearly varying 

WTR source distributions with zero strength 

along trailing edge 

UTL, Same as above for left wing panels 
VTL, 

WTL 

UC , Arrays of velocity components induced 

VC, by vortex panels at control point I 

WC 

USAVE , Velocity component storage arrays 

VSAVE , 

WSAVE 

AC Array of normal velocities induced by 

vortex panels at control point I 

ASAVE Normal velocity storage array 

BC, Velocity tangential to control point 

BT panel I 

UT, Arrays of velocity components induced 

VT , by source panels at control point I 

WT 

AT Array of normal velocities at control 

point I induced by source panels 

DC Array of normal velocities induced by 

vortex panels in diagonal block matrices 
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SUBROUTINES 
CALLED : 


ERROR 
RETURNS : 


VORVEL , 
None 


SORVEL 



PROGRAM NEWORD 


PURPOSE : 


METHOD : 


To revise chordwise panel spacing on the wing 
and to compute new airfoil ordinates by inter- 
polation . 

The program first checks the input data to 
determine if the wing has a round leading edge. 

If so, an array of wing leading edge radii are 
read in. The program then checks if the percent 
chord locations of the panel edges are to be 
redefined. If so, an array of revised chordwise 
locations are read in, otherwise the panel edges 
are used as originally defined. 

For each wing section, the original camber and 
thickness distributions are rewritten as one 
dimensional arrays. A chain of cubic curves having 
continuous first derivatives is fitted between 
each pair of points on these two curves , and the 
four coefficients of the cubic curve calculated 
within each interval. For wing sections having 
round leading edges with infinite leading edge 
slope, the slope at the second percent chord loca- 
tion is calculated by fitting the curve 
z = / 2 px + ax + bx^ through the first three points. 
The program then calculates the coefficients of 
the cubic curves through the remaining points in 
the usual way, starting with the slope determined 
from the first derivative of the above formula. 

The revised camber and thickness ordinates and 
slopes are then calculated at the new chordwise 
locations by the formulas 

Z = Cl + C 2 X + CjX^ + C4X^ 

dz/dx = C 2 + 2 C 3 X + 3ci,x^ 

where the coefficients correspond to the interval 
of the curve under consideration. For wings 
having round leading edges, the formula given in 
the previous paragraph is used in the first interval. 
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USE: 


CALL OVERLAY (LWB, 1, 2) 


Input : 
K1 

NWAF 

NWAFOR 

KWAFOR 

XAF 

XAFK 

TZORD 

WAFORD 

RHO 

Output : 

NWAR 

ZORD 

TORD 

NBA, 

DA 

DZC, 

DZCDX 

DZT, 

DZTDX 

A, B 

C, CC 


Wing leading edge definition parameter 

Number of wing airfoil sections 

Number of ordinates used to define wing 
airfoil section 

Number of ordinates used to define wing 
panel leading and trailing edges. If 
KWAFOR = 0, NWAFOR ordinates are used. 

Array of percent chords for airfoil 
ordinates (NWAFOR values) 

Array of percent chords for panel leading 
and trailing edges (KWAFOR values) 

Array of camber line ordinates 

Array of half-thickness ordinates 

Array of leading edge radii 


Number of intervals in curve 

Array of camber line ordinates 

Array of half-thickness ordinates 

Number and value of derivative at 
initial point on curve 

Array of camber line slopes 
Array of half-thickness slopes 

Coefficients of leading edge curve 

Coefficients of cubic curves in each 
interval 
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TZORK 


Array of revised camber line ordinates 
DZCDXK Array of revised camber slopes ordinates 

WAFORK Array of revised half-thickness ordinates 

DZTDXK Array of revised half-thickness slopes 

SUBROUTINES 
CALLED ; DERIV 

ERROR 

RETURNS : None 
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PROGRAM NEWRAD 


PURPOSE : 
METHOD : 


USE: 


To revise the body meridian line spacing. 

For each body segment, there are three options 
for redefining the meridian lines. Considering 
the first segment, if KRADX(l) = 0, the merid- 
ian lines are not changed. If KRADX(l) is 
positive, the meridian lines are relocated at 
KRADX(l) equally spaced values of the meridian 
angle <)) . If KRADX(l) is negative, an array of 
arbitrary meridian angles is read in. 

If the body has a circular cross section, the 
y and z coordinates are calculated at each 
axial station as follows: 

y = r cos({) 

z = + r sin(|) 

where the body radius r and camber Zq have 
been previously calculated in program CONFIG. 

If the body has an arbitrary cross section, 
the y and z coordinates are obtained by linear 
interpolation at the new values of (j> of the 
original y and z coordinates read in program 
CONFIG. 

The X, y, and z coordinates are written on 
TAPE 10 , and the procedure repeated for the 
remaining body segments . 

CALL OVERLAY (LWB, 1, 4) 

Input ; 

J2TEST Parameter to specify body cross 
section and camber definition 

NFUS Number of body segments 

NRADX, Number of meridian lines on segment 
KRADX 

PHIK Array of meridian angles on segment 



NFORX, Number of axial stations on segment 
NFUSOR 


XFUS 

Array of axial stations on segment 

FUSRAD 

Array of body radii on segment 

ZFUS 

Array of body camber ordinates on 
segment 

SFUS 

Array of y and z coordinates used in 
9-^bitrary cross section definition 

Output : 


KFUS 

Number of body segments 

NF, 

NFU 

Body segment number 

NRAD, 

KRAD 

Number of meridian lines in segment 

KTEST 

Arbitrary body indicator 

NEWPHI 

Logical variable controlling input of 
new meridian angles 

PHIR 

Meridian angle (radians) 

DELE 

Incremental meridian angle 

XB 

Array of axial stations on segment 

YB, 

ZB 

Arrays of y and z coordinates on 
segment 

YF, 

ZF 

Temporary arrays of y and z coordinates 

RAD 

Body radius 

CAM, 

ZC 

Body camber ordinate 

PHIN 

Array of original meridian angles for 
arbitrary cross section body 

MAX 

Maximum number of body axial stations 


SUBROUTINES 
CALLED : 


None 


ERROR 
RETURNS : 


The program will call EXIT if KRAD > 60. 



PROGRAM NUTORD 


PURPOSE : 
METHOD : 


USE: 


To revise chordwise panel spacing on fin, canard 
or tail and compute new airfoil ordinates. 

The program first tests to determine if the 
component under consideration is a fin (vertical 
tail) , a canard, or a horizontal tail. The 
program then initializes the appropriate con- 
stants , and reads in an array of leading edge 
^^dii if the component has a round leading edge. 

Each horizontal or vertical tail component is 
then treated as an additional wing segment , and 
the procedure follows the steps described under 
program NEWORD. 


CALL OVERLAY (LWB, 1, 6) 


Input : 




J4 

Fin definition parameter 


J5 

Tail or canard definition parameter 

K4 

Fin leading edge definition 

parameter 

K5 

Tail or canard leading edge 
parameter 

definition 

NF 

Number 

of fins 


NC 

Number 

of tails and canards 


NFINOR 

Number 

of ordinates defining 

fin airfoil 

NCANOR 

Number 

canard 

of ordinates defining 
airfoils 

tail and 

KFINOR 

Number of ordinates defining 
leading and trailing edges 

fin panel 

KANOR 

Number 

canard 

of ordinates defining tail or 
panel leading and trailing edges 

XAF, 

XT 

Array of percent chords for airfoil 
ordinates 
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XAFK 


TALCR 

TALOR 

RHO 

Output : 

FIN 

NT 

J1 

JL 

KL 

NWAFOR 

KWAFOR 

NWAR 

NDA, 

DA 

ZORD 

TORD 

DZC, 

DZCDX 

DZT, 

DZTDX 

A, B 

C, CC 

TZORK 


Array of percent chords for panel 
leading and trailing edges 

Array of airfoil camber line ordinates 

Array of airfoil half-thickness 
ordinates 

Array of airfoil leading edge radii 

Fin identification variable (logical) 
Number of fins , tails or canards 
Tail or canard camber identifier 
Tail definition integer 

Airfoil leading edge definition integer 

Number of ordinates defining airfoil 

Number of ordinates defining fin, tail 
or canard panel leading and trailing 
edges 

Number of intervals in curve 

Number and value of derivative at 
initial point on curve 

Array of camber line ordinates 

Array of half —thickness ordinates 

Array of camber line slopes 

Array of half —thickness slopes 

Coefficients of leading edge curve 

Coefficients of cubic curves in each 
interval 

Array of revised camber line ordinates 




DZCDXK 

Array 

of 


WAFORK 

Array 

of 


DZTDXK 

Array 

of 

SUBROUTINES 
CALLED : 

DERIV 



ERROR 
RETURNS : 

None 




revised camber line slopes 
revised half— thickness ordinates 
revised half— thickness slopes 
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SUBROUTINE PANEL 


PURPOSE : 


METHOD : 


USE: 


To calculate direction cosines of the normal 
vector, the centroid, area, and inclination 
angles of an arbitrary quadrilateral panel. 

The four corners of the panel are numbered in a 
clockwise direction. A diagonal vector T^ connects 
points 1 and 3, and a diagonal vector T 2 connects 
points 2 and 4. The normal vector N is obtained 
by taking the cross product of these diagonal 
vectors, and the direction cosines determined by 
calculating the projections of this vector in the 
reference coordinate system. The plane of the 
panel is defined to be perpendicular to the normal 
vector and to pass through a point whose co- 
ordinates are the averages of the coordinates of 
the four input points. The input points are 
then projected into the plane of the panel, and 
transformed to the reference coordinate system. 

A new panel coordinate system Cf U is introduced 
with the average point of the panel as origin. 

The coordinates of the centroid and the panel 
area are calculated in this new system, and the 
centroid transformed to the reference system. Two 
angles are used to define the panel orientation. 
The incidence 6 is the angle between the x axis 
and the line of intersection with the panel of a 
plane passing through the x axis and perpendicular 
to the panel. The inclination 9 is the angle 
between the y axis and the line of intersection 
of the panel with the yz plane. These two angles 
are calculated in terms of the direction cosines 
of the normal vector. 

CALL PANEL (IP, IQ, J, K, L, NP , AP) 

Input : 

IP, IQ Panel identification code 

j Panel row number 

K Panel column number 

L Surface identification code 

NP Panel number 


75 




XC 

X coordinate of panel corner point 


YC 

y coordinate of panel corner point 


ZC 

z coordinate of panel corner point 


Output : 



NX, 

NY, 

Direction cosines of the normal vector 


NZ 



AVX, 

AVY, 
AVZ 

Coordinates of average point 


D 

Distance from corner point to plane of 
panel 


XI, 

ETA 

Coordinates of corner point in panel 
coordinate system 


XIO, 

ETAO 

Coordinates of centroid in panel co- 
ordinate system 


XPT, 

YPT, 

ZPT 

Coordinates of centroid in reference 
coordinate system 


DELTA 

Panel incidence angle 


THET 

Panel inclination angle 


AP 

Panel area 

SUBROUTINES 
CALLED ; 

None 


ERROR 
RETURNS ; 

None 
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SUBROUTINE PARTIN 


PURPOSE : 


METHOD : 


USE; 


For wing-body combinations, to invert the 
matrix partitions and store the results on 
TAPE 10. For isolated wings or bodies, to 
solve the boundary condition equations and 
determine the body source strengths or wing 
vortex strengths. 

This subroutine is called only if the order of 
the matrix partition (or the full matrix in 
the case of isolated wings or bodies) does not 
exceed 60. 

For wing-body combinations, the partitions are 
read from TAPE 9, inverted, and the inverse 
matrix written on TAPE 10. 

For isolated wings or bodies , the matrix is 
read from TAPE 9, inverted, and the inverse 
post-multiplied by the normal velocity array 
to obtain the body source strengths or wing 
vortex strengths. 

CALL PARTIN 

Input ; 

NWING 

NBODY 

NDIM 

A, 

D 

NW 
NB 

Output : 

NPANEL Total number of panels 


Number of wing and tail panels 
Number of body panels 
Matrix dimension statement size 
Array of matrix elements 

Array of wing and tail normal velocities 
Array of body normal velocities 


A, 

D 


Array of inverse matrix elements 



TIME 


Elapsed time in seconds 



GW 

Array 

of 

wing and tail vortex strengths 


GB 

Array 

of 

body source strengths 

SUBROUTINES 
CALLED : 

SECOND 




ERROR 
RETURNS : 

None 
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SUBROUTINE PRESS 


PURPOSE : 
METHOD : 


To compute the pressure coefficient at a 
panel control point. 

The u, V, and w components of velocity are 
given in terms of the reference coordinate 
system. They are transformed into a new co- 
ordinate system aligned with the free stream 
velocity vector, and used to determine the 
pressure coefficient by means of the exact 
isentropic formula: 

-1} 

where = (1 + u' ) ^ + + w’ ^ 

u' = u cosa + w sina 


w' = w cosa - u sina 


If the Mach number M is zero, the pressure 
coefficient is determined by the incompressible 
formula : 


CP = 1 - 


The subroutine also calculates the stagnation 
pressure coefficient, the critical pressure 
coefficient, and the vacuum pressure coefficient 
by the following formulas : 


CPSTAG 

CPCRIT 

CPVAC 


YM^ 

2 


{t 

h 


1 + 


Y-1 


Y-1 


YM^ ^ Y+1 Y+1 


m2] - 1 I 

m2] - 1 I 


2 

YM2 
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USE: 


CALL PRESS (NP, XMACH , ARA, U, V, W, CPP , 
CPSTAG, CPCRIT, CPVAC) 



Input : 



NP 

Panel number 


XMACH 

Mach number 


ARA 

Angle of attack in radians 


u. 

Velocity components at panel control 


V, 

point 


W 



Output : 



CPP 

Pressure coefficient at panel control 
point 


CPSTAG 

Stagnation pressure coefficient 


CPCRIT 

Critical pressure coefficient 


CPVAC 

Vacuum pressure coefficient 

SUBROUTINES 

CALLED: 

None 


ERROR 
RETURNS : 

None 
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SUBROUTINE SCAMP4 


PURPOSE : 


METHOD : 


Given a set of n points (x^^, whose abscissae 

form a strictly monotone sequence, a first or 
second derivative at Xj , and a first or second 
derivative at x^ , to find the smoothest possible 
curve passing rigorously through the given points, 
satisfying the specified boundary conditions, and 
possessing continuous first and second derivatives. 
The criterion for smoothness is the minimization 
of the integral of the square of the second deriv- 
ative, from Xj to Xj^ , over all functions having the 
stated properties. Accordingly, the curve found is 
a chain of cubics, i.e., a separate cubic defined 
on each interval (xj^, Xj^+j). The coefficients of 
each such cubic are explicitly found in the form 

y = Cg + CjX + C 2 X^ + CjX^ 

The most economical (in time and space) and most 
accurate method of finding such a chain of cubics 
is to solve first for the n slopes y^^ of the curve. 
This is done by the composite cubic subroutine 
COMCU, which solves an nth order linear system, the 
coefficient matrix of which is tridiagonal. Having 
found the slopes at each of the n given Xj^ , one can 
determine the coefficients of each cubic separately 
by using CUBIC2, which finds the cubic through two 
points, being given the slope at each. The coeffi- 
cients of all the n -1 cubics can be obtained by 
using the subject routine (SCAMP4) which serves as 
a vehicle for calling COMCU (once) and CUBIC2 (n-1 
times) . SCAMP4 has an option to compute the re- 
quired boundary conditions (first or second deriv- 
atives at the end points) if these are not known by 
the calling program; in this case, the computation 
of first derivatives at x^ and x^ is recommended. 

The cubic coefficients found by SCAMP4 are either 
stored in a 4 by n-1 array or are arranged in the 
composite curve format, i.e., in a single linear 
array where each segment is specified by a block of 
seven consecutive words: x^ , Xj^^_j , 3., Cg , Cj, C 2 , 

C 3 . The calling program should dimension the coef- 
ficient array as a doubly subscripted variable in 
the former case and singly subscripted in the latter 
case . 
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USE: 


CALL SCAMP4 (X, Y, N, NDA, NDB, DA, DB , C, S, M) 
Input ; 


SUBROUTINES 
CALLED : 


X 

Y 

N 

NDA 


NDB 

DA 


DB 


M 


Array of x values 
Array of y values 
Number of points 

The order (1 or 2) of the derivative to 
be given at X(l) . If derivative is to 
be computed by SCAMP4 , NDA < 0 . 

The order of the derivative to be given 
at X(N). Similar to NDA. 

The value of the derivative at X(l) . 

If derivative is to be computed by 
SCAMP4, leave blank. 

The value of the derivative at X(N). 
Similar to NDA. 

Code . 

7 ^ 12, if the cubic chain coefficients 
are to be stored in a doubly dimen- 
sioned 4 X (N-1) array. 

= 12, if the cubic chain coefficients 
are to be stored in a singly dimen- 
sioned array 


Output t 

C 

S 

M i 


Array of cubic chain coefficients 
Array of first derivatives 
Error return 
= 0 - success 
^ 0 - error detected 


COMCU, CUBIC2, DERIVl, DERIV2 
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ERROR 
RETURNS : 


M = -1 indicates N < 2. 1 ^ M 
an error return from COMCU. M 
error return k on the i^h call 
(M = 100 . j + k) . 


^ 7 indicates 
large indicates 
to CUBIC2 



SUBROUTINE SECOND 


PURPOSE ; 
METHOD : 

USE; 

SUBROUTINES 
GALLED ; 

ERROR 
RETURNS : 

RESTRICTIONS ; 


To return elapsed CPU time in seconds. 

Control Data Corporation SCOPE Library 
subroutine . 

CALL SECOND (TIME) 

Output : 

TIME Elapsed CPU time in seconds 
None 


None 

Limited to CDC computers using SCOPE 3.0 
operating system and library tape. 
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PROGRAM SOLVE 


PURPOSE : 


METHOD : 


To solve for the strengths of the body sources 
and wing vortices which satisfy the boundary 
condition of tangential flow at the panel con- 
trol points , and to determine the corresponding 
pressure distribution, forces and moments on 
the configuration. 

The program first calculates the array of normal 
velocities required to satisfy the boundary 
conditions at the wing and body panel control 
points. The panel inclination angles 0 and 6 
are obtained from the geometry arrays on TAPE 7 , 
and the angle of attack a from COMMON block 
PARAM. 

If the planar boundary condition and wing thick- 
ness options have been selected, the program 
next computes the normal velocities induced on 
the body and non-coplanar wing or tail segments 
by the wing source distribution. These normal 
velocities are subtracted from those previously 
calculated to obtain the resultant normal veloc- 
ities at the control points. 

The coefficients of the equations to be solved 
have previously been stored in row order on 
TAPE 9. Three different procedures are followed 
to solve the equations depending on the order 
of the matrix of coefficients. If the configu- 
ration to be analyzed consists of an isolated 
wing or body, and the order of the matrix does 
not exceed 60, the equations are solved in sub- 
routine PARTIN by direct inversion of the matrix. 
If the configuration consists of a wing-body 
combination , and the order of the wing and body 
partition does not exceed 60, subroutine PARTIN 
inverts the diagonal partitions of the matrix 
and writes the inverse matrices on TAPE 10 . An 
iterative procedure described in subroutine 
ITRATE is then applied to solve the equations. 

For any configuration for which the order of 
the wing or body partition exceeds 60 , the di- 
agonal blocks of the matrix are read from TAPE 7 , 
inverted, and written on TAPE 10 by subroutine 
DIAGIN. Subroutine ITRATE is then called to 
solve the resulting equations by an iterative 
procedure . 
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Once the strengths of the source and vortex 
distributions have been determined, the program 
calculates the three components of velocity and 
pressure coefficient at each panel control point 
starting with the body panels. The velocity 
components corresponding to unit strength source 
and vortex distribution are obtained from TAPE 8 
The first three records on this file contains 
the velocity components at body control points 
induced by the body source panels , the wing 
source panels (if present) , and the wing vortex 
panels . The program multiplies these by the 
corresponding source and vortex strength, and 
sums the products to obtain the resultant veloc- 
ity component arrays at each body control point. 
The magnitude of the normal velocity at the body 
•^o^^trol points is also calculated at this point. 
If the absolute value of the print option is 
9^®3.ter than one, the three components of veloc— 
ity and the normals are written on the output 
file. The program then calls sxibroutine PRESS 
to obtain the pressure coefficients at the body 
panels, and subroutine FORMOM to integrate the 
pressures and calculate the force and moment 
acting on the body. 

The velocity components at the wing and tail 
panel control points are computed next. The 
remaining three records containing the velocity 
components at wing and tail control points in- 
duced by the body source panels, the wing source 
panels (if present) and the wing vortex panels 
are read from TAPE 8 . The program multiplies 
these by the corresponding source and vortex 
strengths and sums the products to obtain the 
resultant velocity component arrays at the wing 
and tail panel control points. If the absolute 
value of the print option is greater than one, 
the velocity component arrays are written on 
the output file. The program then calls sub- 
routine PRESS to obtain the pressure coeffi- 
cients, and subroutine FORMOM to calculate the 
force and moment acting on the wing. 


If the planar boundary condition option has been 
selected, two passes through this section are 
required to obtain the velocity components, 
pressure and forces on both upper and lower 
surfaces . 


USE: 


The program writes the values of the stagnation 
pressure coefficient, the critical pressure 
coefficient, the vacuum pressure coefficient, 
and the elapsed time on the output file prior 
to returning. 

CALL OVERLAY (LWB, 3, 0) 

Input : 

NBODY Number of body panels 

NWING Number of wing and tail panels 

NWTHK Number of wing and tail panel source 

distributions 

NMAX Maximum order of diagonal block matrices 
PRINT Print option parameter (integer) 

MACH Mach number 

MATIN Matrix inversion parameter 

LBC Planar boundary condition option 

parameter (logical) 

THK Wing thickness option parameter (logical) 

ALPHA Angle of attack in degrees 

CHORD Array of wing and tail panel chords 

DZTDX Array of wing and tail panel half- 

thickness slopes 

ARRAY Wing or body panel geometry arrays 

DELTA Array of panel incidence angles 

THET Array of panel inclination angles 

UA, Arrays of velocity components induced 

VA, by unit strength source and vortex 

WA distributions 

Output : 

EM Mach number 
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NPASS 


Number of passes through program 


COMPT 

ALP 

SINAL, 

COSAL 

SINT, 

COST 

TANDEL 

NW 

NB 

U, 

V, 

W 

NS 

GB 

GW 

CP 

CP STAG 
CPVAC 
CPCRIT 
SGN 


Wing or body component indicator 
(integer) 

Angle of attack in radians 

Trigonometric functions of angle of 
attack 

Trigonometric functions of inclination 
angle 6 

tan 6 

Array of normal velocities required 
to satisfy boundary conditions at 
wing control points 

Array of normal velocities required 
to satisfy boundary conditions at 
body control points 

Arrays of velocity components at 
control points 


Array of normal velocities at control 
points 

Array of body panel source strengths 

Array of wing and tail panel vortex 
Strengths 

Array of pressure coefficients at control 
points 

Stagnation pressure coefficient 

Vacuum pressure coefficient 

Critical pressure coefficient 

Wing upper and lower surface sign 
parameter 


TIME 


Elapsed time in seconds 



SUBROUTINES 
CALLED ; 


DIAGIN, PARTIN, ITRATE , PRESS, FORMOM, 


ERROR 

RETURNS : None 


SECOND 



SUBROUTINE SORPAN 


PURPOSE : 


METHOD ; 


USE: 


To calculate the three components of velocity 
induced at a given control point by a constant 
source distribution on a quadrilateral panel 
having longitudinal taper and inclined at an 
angle 6 to the free stream direction. 

Formulas for the three components of velocity 
u, V, and w are given in Part I of this report. 
The subroutine obtains the information necessary 
to evaluate these formulas through COMMON block 
BODCOM, and returns the velocity components 
through the subroutine pareuneter list. 

The first step in the calculation procedure is 
the adjustment of the z coordinate of each panel 
corner to ensure that it lies in the plane con- 
taining the panel leading edge and its control 
point. The influence functions F, G, and H are 
then calculated at the specified control point 
for each corner point in sequence. The final 
result is obtained by combining the four values 
of these functions. 


CALL SORPAN (UPM, VPM, WPM) 


Input : 



EM 

Mach number 


SA 

tan 6 (6 is 

panel inclination angle) 

CX 

Panel chord 

length 

XC, 

YC, 

Arrays of corner point coordinates 

ZC 



XI, 

YI, 

Coordinates 

of control point 

ZI 



XJ, 

ZJ 

Coordinates 

of panel control point 
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Output ; 



BT2 , 
BTA 

Mach number parameters 


SM 

Panel side edge slopes 


DX, 

DY, 
DZ 

Coordinates of control point referred 
to corner point 


D 

Compressed distance from corner point 
to control point 


DPM 

Scaled compressed distance from corner 
point to control point 


XPM, 

YMX, 

ZAX, 

AYM 

Transformed control point coordinates 


RPM 

Compressed distance from side edge to 
control point 


TA 

Panel edge slope parameter ,( 1 . + BT2.SA^) 


TAM 

Panel edge slope parameter, 
(1. + BT2 (SA^ + SM^) ) 


E, 

F, 

G , 
H 

Corner point influence functions 


E14, 
F14 , 
G14 

Combined corner point influence functions 


R4PI 

Reciprocal of 4 tt (or 2 tt if EM > 1) 


UPM, 

VPM, 

WPM 

Velocity components at control point in 
panel coordinate system 

SUBROUTINES 
CALLED ; 

None 


ERROR 
RETURNS : 

None 
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SUBROUTINE SORVEL 


PURPOSE : 


METHOD : 


USE: 


To calculate the three components of velocity 
induced at a given control point by constant 
and linearly varying source distributions on 
a swept quadrilateral panel. The subroutine 
calculates the velocity components induced by 
one corner of the panel. 

Formulas for the three components of velocity 
UC, VC, WC induced by a constant source distri- 
bution, and UL, VL, WL induced by a source 
distribution having a linear chordwise variation 
are given in Part I of this report. The sub- 
routine obtains the information necessary to 
evaluate these formulas through COMMON block 
COMPS, and returns the velocity components 
through the parameter list. It is assumed that 
the Gothert Rule compressibility transformation 
has been applied to the geometrical data prior 
to calling this subroutine. Both sets of veloc- 
ity components are expressed in terms of the 
influence functions FI, Gl, G2 , and G3 which 
depend only on the geometrical relationship of 
the control point to the corner point, the 
sweep angle, and Mach number. 

The subroutine contains three separate branches 
for evaluating the velocity components , corre- 
sponding to the general case, a special case 
for supersonic leading edges, and a special 
case if the control point lies in the plane of 
the panel. 

CALL SORVEL (UC, VC, WC , UL, VL, WL) 

Input : 

DELTAY , y and z coordinates of control point 
DELTAZ in reference coordinate system 

X, Coordinates of control point with 

Y, reference to corner point 

Z 

B Leading edge slope 

BROS Leading edge slope sign parameter 

(logical) 
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SUB 


Subsonic Mach number parameter (logical) 



COST, 

SINT 

Sine and cosine of panel dihedral 
angle 9 


Output ; 



SUP 

Supersonic Mach number parameter 
( logical) 


SUPLE 

Supersonic leading edge parameter 
(logical) 


BNEG 

Leading edge slope sign parameter 
( logical) 


BTERM 

Leading edge slope parameter 


D 

Distance from control point to 
corner point 


R 

Distance from control point to 
side edge 


TZ 

Distance from control point to 
leading edge 


Cl, 
T2 , 
T3, 
AT 3 

Geometrical parameters 


FI, 
Gl, 
G2 , 
G3 

Influence functions 


UC, 

VC, 

WC 

Velocity components at control point 
induced by constant source distribution 


UL, 

VL, 

WL 

Velocity components at control point 
induced by source distribution with 
linear chordwise variation 

SUBROUTINES 
CALLED : 

None 


ERROR 
RETURNS : 

None 
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PROGRAM TALPAN 


PURPOSE : 
METHOD : 


USE: 


To revise the spanwise panel spacing on fin, 
canard, or tail and compute the panel geometry. 

The program first tests to determine if the 
component under consideration is a fin (vertical 
tail) , a canard, or a horizontal tail. The 
program initializes the appropriate constants, 
then rewinds TAPE 7 , reads the wing geometry 
arrays from that file, and stores them in COMMON 
block POINT. Each horizontal or vertical tail 
component is then treated as an additional wing 
segment, following the steps described under 
subroutine WNGPAN. 

At the completion of the calculation, the com- 
bined wing and tail geometry arrays are stored 
in COMMON block POINT, and the entire sequence 
of arrays are written as a single record back 
on TAPE 7 . The augmented CHORD and SLOPE arrays 
are also written on TAPE 7 at this point. The 
remaining wing and tail geometry parameters are 
stored in COMMON blocks PARAM and SEG. Finally, 
if the print option is positive, the fin, canard 
or tail panel corner point coordinates, control 
point coordinates, inclination angles, areas, 
and chords are written on the output file. 

CALL OVERLAY (LWB, 1, 7) 


Input : 



LBC 

Boundary condition option 

(logical) 

PRINT 

Print option 


K4 

Fin definition parameter 


K5 

Tail or canard definition 

parameter 

NF 

Number of fins 


NK 

Number of tails or canards 

NCPT 

Number of control points ' 

on wing 

NWING , 

Number of panels in wing 

(initial value) 

NINIT 
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KF 


KAN 

YK 

KFINOR 

KANOR 

TALORG 

XAFK 

WAFORK 

TZORK 
DZTDXK 
DZCDXK 
Output ; 
Note : 


FIN 

NTAL 

NT 

KK 

KL 

KWAF 


Number of streamwise panel edges on 
fin 

Number of streamwise panel edges on 
tail or canard 

Array of spanwise locations of fin or 
tail panel streamwise edges 

Number of ordinates defining fin panel 
leading and trailing edges 

Number of ordinates defining tail or 
canard panel leading and trailing 
edges 

Array of origin and chord length of 
each fin, tail or canard airfoil 
(x, y, z, c) 

Array of percent chords for panel 
leading and trailing edges 

Array of airfoil half-thickness 
ordinates 

Array of airfoil camber ordinates 
Array of airfoil half— thickness slopes 
Array of airfoil camber slopes 


In the following descriptions , the 
words tail or tail segment may refer 
to any fin, canard, or tail component 

Fin identification variable (logical) 

Number of fins , tails , or canards 

Tail segment niimber 

Fin identification integer 

Leading edge identification integer 

Number of panel streamwise edges on 
tail segments 



KWAFOR 

WAFORG 

NWING 

NCPT 

NP 

NC 

NSEG 

NROW 

NCOL 

KOL 

BL, 

BLE 

BT, 

BTE 

TH 

SINS 

COSS 

NWT 

XK 

ZK 

CK 


Number of ordinates defining panel 
leading and trailing edges on tail 
segment 

Array of origin and chord length of 
each tail segment airfoil (x, y, z, c) 

Total number of wing and tail panels 

Total number of wing and tail control 
points 

Panel number 

Control point number 

Number of wing and tail segments 

Number of rows of panels in segment 

Number of columns of panels in segment 

Number of columns of panels on wing 
and tail 

Arrays of segment leading edge slopes 


Arrays of segment trailing edge slopes 


Array of segment dihedral angles 

Array of sine of segment dihedral angle 

Array of cosine of segment dihedral 
angle 

Array of wing and tail indicator 
parameters 

Array of x-coordinates of origins of 
tail panel streamwise edges 

Array of z-coordinates of origins of 
tail panel streamwise edges 

Array of chord lengths of tail panel 
streamwise edges 


CL 

L 

SJ 

IP, 

IQ 

XC, 

YC 

ZC 

ZU 

CR 

CT 

RI, 

RO 

XLE, 

XLEW 

XTE 

CHORD 

SPN, 

SPNW 

AREA 


Chord length of tail panel streamwise 
edge divided by one hundred 

Tail surface indicator 

L = 1 indicates upper (inner) surface 
L = 2 indicates lower (outer) surface 
Tail surface sign parameter 
Panel identification constants 


Arrays of tail panel corner point x 
and y coordinates 

Array of tail panel corner point z 
coordinates or lower (outer) surface 
coordinates for the non planar boundary 
condition option 

Array of upper (inner) surface co- 
ordinates for the non planar boundary 
condition option 

Panel root chord 

Panel tip chord 

Centroid ratios 


X coordinate of intersection of panel 
leading edge with streamwise line 
through centroid 

X coordinate of intersection of panel 
trailing edge with streamwise line 
through centroid 

Array of panel chord lengths passing 
through centroids 

Array of panel spans 


Array of panel areas 



SUBROUTINES 
CALLED : 

ERROR 
RETURNS : 


XPT, Array of panel control point coordinates 

YPT, 

ZPT 

THET Array of panel dihedral angles 

DZCDX Array of tail camber slopes at panel 

edges (planar boundary condition option) 

DELTA Array of tail camber slopes at panel 
control points (planar boundary 
condition option) or panel incidence 
angle (non planar boundary condition 
option) 

DZTDX, Array of tail half-thickness slopes at 

SLOPE panel edges (planar boundary condition 
option) 

SLE Leading edge slope for round leading 

edge airfoils 

XE Array of x coordinates of panel control 

points 

XS , Array of point source origins (non planar 

YS, boundary condition option) 

ZS 


PANEL 

The program calls EXIT if NWING > 600 
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SUBROUTINE TRANS 


PURPOSE : 


METHOD : 


USE: 


To transform the three components of velocity 
from the panel coordinate system to the refer- 
ence coordinate system, to combine the contri- 
butions of the left and right wing panels, and 
to calculate the normal velocity at the control 
point . 

The axial and vertical velocity components are 
transformed by a rotation of the coordinate 
system about the horizontal axis by the angle o . 
The axial velocity components induced by the 
left and right wing panels are added directly 
to determine the resultant axial velocity u. 

Two additional coordinate rotations about the x 
axis are required before the v and w components 
induced by the left and right wing panels can 
be combined. The first rotation transforms the 
V and w components from the influencing panel 
coordinate system to the control point panel 
coordinate system, and the second transforms 
the combined normal and tangential velocity at 
the control point to v and w velocity components 
in the reference coordinate system. 

CALL TRANS (UR, VR, WR, UL, VL, WL, U, V, W, A) 


Input : 

UR, Three components of velocity at control 

VR, point in right wing panel coordinate 

WR system 


UL , Three components of velocity at control 

VL, point in left wing panel coordinate 

WL system 


Output : 

U, Three components of velocity at control 

point in reference coordinate system 


Normal velocity at control point 
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SUBROUTINES 


CALLED : 

None 

ERROR 


RETURNS ; 

None 
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SUBROUTINE TRAP 


PURPOSE : 
METHOD : 

USE: 


SUBROUTINES 
CALLED : 

ERROR 
RETURNS : 


To evaluate an integral by the trapezoidal 
rule . 

The X and y coordinates of a curve are read 
and the integral obtained by sununing the areas 
within each interval for i = 1, NT 

NT 

SUM = j ^ (^i " ^i-1^ ^^i ^i-l^ 

i:=l 

CALL TRAP (XT, YT, SUM, NT) 

Array of x coordinates (abscissa) 

Array of y coordinates (ordinates) 
Number of coordinates 

Output : 

SUM Integral J y dx by trapezoidal rule 

None 

None 


Input : 
XT 
YT 
NT 



PROGRAM US S AERO 


PURPOSE : 
METHOD : 

USE: 

PROGRAMS 
CALLED : 


This program controls the sequence of calculations 
required to determine the aerodynamic character- 
istics of wing-body-tail configurations in subsonic 
or supersonic potential flow. 

The input card deck is read and listed on the out- 
put file. The three primary overlay programs GEOM, 
VELCMP , and SOLVE are then called in sequence to 
perform the remaining calculations. The lengths 
of the principal COMMON blocks are also specified 
in this program. 


OVERLAY 

(LWB, 

0, 0) 


LWB is 1 

overlay 

■ file 

name . 

OVERLAY 

(LWB, 

1, 0) 

(GEOM) 

OVERLAY 

(LWB, 

2, 0) 

(VELCMP) 

OVERLAY 

(LWB, 

3, 0) 

(SOLVE) 
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PROGRAM VELCMP 


PURPOSE : 


METHOD : 


To compute the velocity components u, v, and w 
at panel control points, and form the aero- 
dynamic influence coefficient matrices. 

The program reads the Mach number and angle of 
attack from the input file. If the Mach number 
is negative, or the same as the previous case, 
a return is executed. Otherwise, the program 
proceeds to compute the velocity components . 

For wing alone, or wing-body configurations, a 
preliminary calculation is made to determine if 
the wing control points require relocation, and 
to compute the number and size of the wing di- 
agonal blocks for later use in the matrix calcu- 
lations. For wing-body configurations, the body 
geometry is first placed in temporary storage 
on TAPE 10. The program then proceeds to recal- 
culate the chordwise locations of the wing control 
points for wings having supersonic edges , provided 
the planar boundary condition option has been 
selected. (An edge is defined to be supersonic 
if the component of the Mach number normal to 
the edge is greater than one.) Considering one 
wing segment at a time, the program tests to 
determine if either the leading or trailing edge 
is supersonic. If all edges are subsonic, the 
control points retain their original locations 
at the panel centroids. If the leading edge is 
subsonic and the trailing edge is supersonic, the 
control points in a given column of panels are 
adjusted uniformly between the centroid of the 
leading edge panel and the trailing edge of the 
last panel in the column. If both edges are 
supersonic, the control points are relocated at 
the panel leading edges, and the trailing edge of 
the last panel in the column. A wing supersonic 
trailing edge indicator array is also computed at 
this point in the program. The revised control 
points are stored in COMMON block POINT, and the 
entire wing geometry array written on TAPE 7. 

The body geometry temporarily stored on TAPE 10 
is then rewritten on TAPE 7 following the wing 
geometry arrays. 
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USE: 


The velocity component calculations are sub- 
divided into four steps. For wing-alone or 
body-alone configurations, only the first step 
is necessary, otherwise all four steps are 
included. Each step involves the calculation 
of the influence coefficients of one partition 
of the complete aerodynamic matrices. The first 
partition gives the influence of the body panels 
at the body control points (or the influence of 
the wing panels at the wing control points for 
wing-alone configurations) . The second partition 
gives the influence of the wing panels at the 
body control points , the third gives the influ- 
ence of the body panels at the wing control 
points , and the fourth partition gives the 
influence of the wing panels at the wing control 
points for wing-body configurations . The pro- 
gram calculates the partition number, reads 
the appropriate geometry arrays from TAPE 7, 
and calls the wing or body panel velocity com- 
ponent program to obtain the influence coeffi- 
cients . 

On completion of the influence coefficient 
calculations, and if the order of any partition 
is greater than 60, the program writes the di- 
agonal blocks of the aerodynamic matrix on TAPE 7 
(following the geometry arrays) in preparation 
for the iterative solution. The number and size 
of the body diagonal blocks is calculated at 
this time, and stored with the wing diagonal 
block data and other matrix constants in COMMON 
block VELCOM. 


CALL OVERLAY (LWB, 2, 0) 


Input : 


LBC 

Linear boundary condition option para- 
meter (logical) 

PRINT 

Print option parameter 

MACH 

Mach number (real) 

ALPHA 

Angle of attack 

NMAX 

Maximum order of diagonal block matrices 
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NBODY 

KFUS 

KRADX 

KFORX 

NWING 

NCPT 

NSEG 

NROW 

NCOL 

BL 

BT 

ARRAY 

CHORD 

DZTDX 

XLE 

XPT, 

YPT, 

ZPT 

DELTA 

D 

Output : 
MATIN 
SUB 
SUBLE 


Number of body panels 

Number of body segments 

Array of body panel meridian lines in 
segment 

^jrray of body panel axial stations in 
segment 

Number of wing and tail panels 

Number of wing and tail control points 

Number of wing and tail segments 

Number of rows of panels in segment 

Number of columns of panels in segment 

Leading edge sweep of wing segment 

Trailing edge sweep of wing segment 

Wing or body geometry arrays on TAPE 7 

Array of wing panel chords 

Array of wing thickness slopes 

Array of chordwise locations of wing 
control points 

Arrays of panel control point coordinates 

Array of panel incidence angles 
Array of diagonal block matrix 


Matrix inversion indicator 
Subsonic indicator (logical) 

Subsonic leading edge indicator (logical) 
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SUBTE 

EM 

BETA 

NPOINT 

NPANEL 

NWBLOK 

NBBLOK 

NWROW, 

NK 

NBROW 

NC 

NR 

BLE 

BTE 

IT 

XPT 

XCPT 

LOCPT 

NPART 

XBT, 

YBT, 

ZBT 

TIME 


Subsonic trailing ©dge indicator 
(logical) 

Mach number 

Mach number parameter 

Number of control points 

Total number of panels 

Number of diagonal block matrices in 
wing partition 

Number of diagonal block matrices in 
body partition 

Number of rows in wing diagonal block 
matrix 

Number of rows in body diagonal block 
matrix 

Number of columns of panels in segment 

Number of rows of panels in segment 

Leading edge sweep of wing segment 

Trailing edge sweep of wing segment 

Array of wing supersonic trailing edge 
indicators 

Array of x coordinates of wing control 
points 

Array of chord fractions for control 
point location 

Array of control point location 
indicators 

Matrix partition number 

Temporary array of panel control point 
coordinates 

Elapsed time in seconds 



PROGRAMS 
CALLED : 


OVERLAY 

(LWB, 

2, 

1) 

(BODVEL) 

OVERLAY 

(LWB, 

2, 

2) 

(LINVEL) 

OVERLAY 

SECOND 

(LWB, 

2, 

3) 

(WNGVEL) 


ERROR 

RETURNS ; None 



SUBROUTINE VORPAN 


PURPOSE : 


METHOD : 


To calculate the three components of velocity 
induced at a given control point by constant 
and linearly varying vortex distributions on 
a swept quadrilateral panel. In addition, the 
subroutine calculates the three components of 
velocity induced by the concentrated vortex 
lying along the downstream extension of the 
inboard edge, and the vortex sheet located down- 
stream of the trailing edge between the inboard 
edge and the intersection of the leading and 
trailing edges of the panel. 

Formulas for the velocity components UC, VC, WC 
induced by a constant vortex distribution, UL, 

VL, WL , and ULT , VLT, WLT induced by the leading 
and trailing edge corners respectively of a 
vortex distribution having a linear chordwise 
variation are given in Part I of this report. 

The subroutine obtains the information necessary 
to evaluate these formulas through COMMON block 
COMPS, and returns the velocity components through 
the parameter list. It is assumed that the 
Gothert Rule compressibility transformation has 
been applied to the geometrical data prior to 
calling this subroutine. 

The subroutine first performs the coordinate 
transformations and calculates the geometrical 
parameters required by the velocity component 
formulas . It then evaluates the downwash veloc- 
ity induced by the trailing vortex sheet by 
numerical integration. Eleven chordwise stations 
are used in the trapezoidal rule integration. 

Three separate branches are provided for evalu- 
ating the velocity coefficients. The first 
branch is a special case for supersonic leading 
edges, the second contains the formulas for the 
general case, and the third contains special 
formulas used if the control point lies in the 
plane of the panel. In the latter two branches, 
the velocity components are expressed in terms 
of the six influence functions Fl, Gl, G2 , G3, 

HI , and H2 which depend on the geometrical 
relationship of the control point to the corner 
point, the leading edge sweep angle, and the 
Mach number. 
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USE: 


The V and w velocity components induced by 
the vortex sheet in the wake are expressed 
in terms of the influence functions Fl and 
H2, while those induced by the concentrated 
vortex in the wake are expressed in terms of 
the parameter C6 . The wake vortices induce 
no axial component of velocity. 

CALL VORPAN (UC, VC, WC , UL, VL, WL, ULT , 
VLT, WLT, VE, WE, VA, WA) 


Input : 

DELTAY, y and z coordinates of control points 
DELTAZ in reference coordinate system 

X, Coordinates of control point with 

Y, reference to leading edge corner 

Z point 

A Difference between leading and trailing 

edge slopes 

B Leading edge slope 

C Panel chord length along inboard edge 

EPOS Leading edge slope sign parameter 

(logical) 

SUB Subsonic Mach number parameter (logical) 

LBC Planar boundary condition option 

parameter (logical) 

COST, Sine and cosine of panel dihedral 

SINT angle 6 

ML Panel leading or trailing edge indicator 

MAX Number of arguments in numerical evalu- 

ation of downwash integral 


Output : 

SUP Supersonic Mach number parameter (logical) 

SUPLE Supersonic leading edge parameter 

( logical) 
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AB, Geometrical parameters 

AZ , 

CC 

R Distance from control point to side 

edge 

D Distance from control point to corner 

point 

E Distance from control point to inter- 

section of leading and trailing edges 
of panel 

T8 

TZ Distance from control point to 

leading edge squared 

B1 , Mach number parameters 

SBl 

B2 Geometrical parameters 

•I- 

B4 

Cl Geometrical parameters 

+ 

C6 

Tl Geometrical parameters 

+ 

T9 

XI Array of x coordinates used in evalu- 

ation of downwash integral 

Q, 

QX Array of arguments used in evaluation 

of downwash integral 

WQ Value of downwash integral at leading 

edge 

WQT Value of downwash integral at trailing 

edge 

SL, Intermediate values of velocity functions 

TL at leading edge 
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ULS, 

VLS, 

WLS, 

TT, 

TLT 

Intermediate values of velocity functions 
at trailing edge 


FI, 

Gl, 

G2, 

G3, 

HI, 

H2 

Influence functions 


VS, 

WS 

Intermediate values of velocity 
components 


UC, 

VC, 

wc 

Velocity components at control point 
induced by constant vortex 
distribution 


UL, 

VL, 

WL 

Velocity components at control point 
induced by the leading edge corner of 
a vortex distribution with linear 
chordwise variation 


ULT, 

VLT, 

WLT 

Velocity components at control point 
induced by the trailing edge corner of 
a vortex distribution with linear 
chordwise variation 


VA, 

VB, 
WA, 
WB 

Velocity components at control point 
induced by the trailing vortex sheet 


VE, 

VD, 

WE, 

WD 

Velocity components at control point 
induced by concentrated vortex in 
wake 

SUBROUTINES 
CALLED ; 

TRAP 


ERROR 
RETURNS : 

None 
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SUBROUTINE VORVEL 


PURPOSE : 


METHOD : 


To calculate the three components of velocity 
induced at a given control point by constant 
and linearly varying vortex distributions on 
a swept quadrilateral panel. The subroutxne 
calculates the velocity components induced by 
the leading and trailing corners of one edge 
of the panel. 

Formulas for the velocity components UC, VC, 

WC induced by a constant vortex distribution, 
UL, VL, WL, and ULT, VLT , WLT induced by the 
leading and trailing edge corners respectively 
of a vortex distribution having a linear chord 
wise variation are given in Part I of this 
report. The subroutine obtains the information 
necessary to evaluate these formulas through 
COMMON block COMPS, and returns the velocity 
components through the parameter list. It is 
assumed that Gothert Rule compressibility 
transformation has been applied to the geo- 
metrical data prior to calling this subroutine. 


The subroutine first performs the coordinate 
transformations and calculates the geometrical 
parameters required by the velocity component 
formulas. It then evaluates the downwash veloc- 
ity induced by the trailing vortex sheet by 
numerical integration. Eleven chordwise stations 
are used in the trapezoidal rule integration. 

Three separate branches are provided for evalu- 
ating the velocity coefficients. The first 
branch is a special case for supersonic leading 
edges , the second contains the formulas for the 
general case, and the third contains special 
formulas used if the control point lies in the 
plane of the panel. In the latter two branches , 
the velocity components are expressed in terms 
of the six influence functions Fl, Gl, G2, G3 , 

HI, and H2 which depend on the geometrical 
relationship of the control point to the corner 
point, the leading edge sweep angle, and the 
Mach number. 
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USE: 


CALL VORVEL (UC, VC, WC, UL, VL, WL, 
ULT, VLT, WLT) 


Input : 

DELTAY, y and z coordinates of control points 
DELTAZ in reference coordinate system 

X, Coordinates of control point with 

Y, reference to leading edge corner 

Z point 

A Difference between leading and trailing 

edge slopes 

B Leading edge slope 

C Panel chord length along inboard edge 

EPOS Leading edge slope sign parameter 

(logical) 

SUB Subsonic Mach number parameter (logical) 

LBC Planar boundary condition option 

parameter (logical) 

COST, Sine and cosine of panel dihedral 

SINT angle 0 

ML Panel leading or trailing edge indicator 

MAX Number of arguments in numerical evalu- 

ation of downwash integral 


Output : 

SUP Supersonic Mach number parameter (logical) 

SUPLE Supersonic leading edge parameter 

( logical) 

AB, Geometrical parameters 

AZ, 

CC 

R Distance from control point to side edge 

D Distance from control point to corner 

point 
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E 


T8 

TZ 


Bl, 

SBl 

B2 

B4 

Cl 

+ 

C6 

T1 

+ 

T9 

XI 


Q, 

QX 

WQ 


WQT 


SL, 

TL 

ULS, 

VLS, 

WLS, 

TT, 

TLT 

FI, 
Gl, 
G2 , 
G3, 
HI, 
H2 


Distance from control point to inter- 
section of leading and trailing edges 
of panel 

Distance from control point to 
leading edge squared 

Mach number pareuneters 


Geometrical parameters 


Geometrical parameters 


Geometrical parameters 


Array of x coordinates used in evalu- 
ation of downwash integral 

Array of arguments used in evaluation 
of downwash integral 

Value of downwash integral at leading 
edge 

Value of downwash integral at trailing 
edge 

Intermediate values of velocity functions 
at leading edge 

Intermediate values of velocity functions 
at trailing edge 


Influence functions 
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vs, Intermediate values of velocity 

WS components 

UC, Velocity components at control point 

VC, induced by constant vortex 

WC distribution 

UL, Velocity components at control point 

VL, induced by the leading edge corner of 

WL a vortex distribution with linear 

chordwise variation 

ULT, Velocity components at control point 

VLT, induced by the trailing edge corner of 
WLT a vortex distribution with linear 

chordwise variation 


SUBROUTINES 


CALLED: 

TRAP 

ERROR 


RETURNS : 

None 
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PROGRAM WNGPAN 


PURPOSE : 
METHdO : 


To revise the spanwise panel spacing on the wing 
and compute the panel geometry. 

The program first checks if the spanwise panel 
spacing is to be revised. If so, an array of 
revised panel edge locations is read in; otherwise, 
the panel edges are used as originally defined. 

The wing panel geometry is established by consid- 
ering regions defined by sequential pairs of the 
originally defined airfoil sections. The leading 
and trailing edge slopes and dihedral angle of the 
region are calculated, and the origins and chord 
lengths of any intermediate panel edges obtained 
by linear interpolation in the spanwise direction. 

The individual panel geometry is then calculated. 
For the planar boundary condition option, the 
corner points and control points are calculated 
in the plane of the wing, while the wing camber 
and thickness slopes at the panel edges are ob- 
tained by a linear interpolation of the slopes 
determined in the program NEWORD. For the non 
planar boundary condition, the corner points and 
control points are calculated on the upper and 
lower surfaces of the wing, and the panel inclina- 
tion angles determined by subroutine PANEL. In 
addition, both options calculate the panel area, 
chord, span, and leading edge x coordinate. 

The same procedure is followed for each of the 
regions between the remaining airfoil sections. 
Prior to each step, the leading and trailing edge 
slopes and dihedral angles of the region are 
compared with those calculated for the previous 
region. If all these quantities are the same, 
the calculation proceeds normally. Otherwise, a 
new wing segment is defined, and the leading and 
trailing edge slopes, sine and cosine of the di- 
hedral angle , and a wing indicator parameter for 
the segment are stored in a special array before 
continuing the calculations. The progrcim also 
computes the number of rows and columns of panels 
in each wing segment, the total number of panels, 
and the total number of segments on the wing. 
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USE; 


The three coordinates of the control points , the 
panel dihedral angles 0, the panel inclination 
angles 6 , the three coordinates of the panel 
corner points, the panel areas, and x coordinates 
of the panel leading edges are stored in the 
COMMON block POINT, and the entire sequence of 
arrays written as a single record on TAPE 7. The 
CHORD and SLOPE arrays are also written on TAPE 7 
at this point. The remaining wing geometry para- 
meters are stored in COMMON blocks PARAM, and SEG. 
Finally, if the print option is positive, the 
corner point coordinates, control point coordinates, 
inclination angles, areas, and chords are written 
in the output file for reference. 

CALL OVERLAY (LWB, 1, 3) 


Input ; 
LBC 
PRINT 
KL 

NWAF 

KWAF 

KWAFOR 

WAFORG 

XAFK 

YK 

WAFORK 

TZORK 

DZTDXK 

DZCDXK 


Boundary condition option (logical) 
Print option 

Leading edge radius parameter 

Number of wing airfoil sections 

Number of wing panel streamwise edges 

Number of ordinates defining wing 
panel leading and trailing edges 

Array of origin and chord length of 
each wing airfoil (x, y, x, c) 

Array of percent chord locations of 
panel leading and trailing edges 

Array of spanwise locations of wing 
panel streamwise edges 

Array of airfoil half-thickness 
ordinates 

Array of airfoil camber ordinates 
Array of airfoil half-thickness slopes 
Array of airfoil camber slopes 
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Output ; 


NWING 

Total nuinber of wing panels 

NCPT 

Total number of control points 

NP 

Panel number 

NC 

Control point number 

NSEG 

Number of wing segments 

NROW 

Number of rows of panels in segment 

NCOL 

Number of columns of panels in 
segment 

KOL 

Number of wing panel strearawise edges 

BL, BLE 

Leading edge slope of wing segment 
array 

BT , BTE 

Trailing edge slope of wing segment 
array 

TH 

Dihedral angle of wing segment array 

SINS 

Sine of segment dihedral angle array 

COSS 

Cosine of segment dihedral angle 
array 

NWT 

Wing indicator parameter array 

XK 

Array of x-coordinates of origins of 
wing panel streamwise edges 

ZK 

Array of z-coordinates of origins of 
wing panel streamwise edges 

CK 

Array of chord lengths of wing panel 
streamwise edges 

CL 

Chord length of wing panel streaunwise 
edge divided by one hundred 

L 

Wing surface indicator 
L = 1 indicates upper surface 
L = 2 indicates lower surface 


SJ 

IP, IQ 
XC, YC 

ZC 

ZU 

CR 

CT 

RI, RO 

XLE, 

XLEW 

XTE 

CHORD 

SPN, 

SPNW 

AREA 

XPT, 

YPT, 

ZPT 

THET 

DZCDX 


Wing surface sign parameter 

Panel identification constants 

Arrays of wing panel corner point x 
and y coordinates 

Array of wing panel corner point z 
coordinates or lower surface z co- 
ordinates for the non planar boundary 
condition option 

Array of upper surface z coordinates 
for the non planar boundary condition 
option 

Panel root chord 
Panel tip chord 
Centroid ratios 

X coordinate of intersection of panel 
leading edge with streamwise line 
through centroid 

X coordinate of intersection of panel 
trailing edge with streamwise line 
through centroid 

Array of panel chord lengths passing 
through centroids 

Array of panel spans 


Array of panel areas 

Arrays of panel control point 
coordinates 


Array of panel dihedral angles 

Array of wing camber slopes at panel 
edges (planar boundary condition 
option) 



DELTA Array of wing camber slopes at panel 

control point (planar boundary 
condition option) or panel incidence 
angle (non planar boundary condition 
option) 

DZTDX, Array of wing half-thickness slopes 

SLOPE at panel edges (planar boundary 

condition option) 

SLE Leading edge slope for round leading 

edge airfoils 

XE Array of x coordinates of panel control 

points 

XS, Arrays of point source origins (non 

YS , planar boundary condition option) 

ZS 


SUBROUTINES 
CALLED : PANEL 

ERROR 

RETURNS: The program calls EXIT if NWING > 600 
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PROGRAM WNGVEL 


PURPOSE : 


METHOD : 


To calculate the three components of velocity 
induced at specified control points by vortex 
panels located on wing or tail surfaces. 

The program first applies the Gothert rule 
compressibility transformation to the tangent 
of the panel inclination angles, and computes 
trigonometric functions of the revised angles. 

If the product 3 tan 6 is greater than one in 
supersonic flow, the panel lies outside the 
Mach cone from its apex, an error message is 
written and the program terminated. 

The three coordinates of the first control 
point, and the corresponding panel inclination 
angles 0 and 6 are read from COMMON block POINT. 
If the control point is on the body, the incli- 
nation angle 9 is obtained from COMMON block 
BTHET . 

The program then computes the influence of 
each panel at the control point. The panels 
on the upper surface of each chordwise column 
are considered first, followed by those on 
the lower surface. This process is repeated 
for each column of panels on a wing segment, 
starting with the inboard panel, and continued 
until all wing and tail segments have been 
included . 

The coordinates of the four corner points of 
the influencing panel are obtained from COMMON 
block POINT in the reference coordinate system. 
They are indexed according to the panel row 
and column numbers . They are first used to 
calculate the leading and trailing edge slopes 
and the chord lengths of the inboard and out- 
board edges of the panel in a panel coordinate 
system lying in the plane of the panel and 
originating at the inboard leading edge corner. 
The control point is also transformed to the 
panel coordinate system, and the velocity com- 
ponents induced at the control point by each 
of the four corners computed by subroutine 
VORPAN. The subroutine is called twice for 
each corner point to obtain the contributions 
of both left and right wing panels. 
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The contribution of a wake consisting of two 
concentrated edge vortices with a constant 
strength vortex sheet between them is calcu- 
lated following the last panel in each column. 
The wake vortices are all oriented in a stream- 
wise direction, and are assumed to lie in a 
plane parallel to the reference axis and con- 
taining the trailing edge of the last panel 
in the column. The velocity components at the 
control point induced by the upstream corners 
of the wake are obtained by four additional 
calls to VORPAN. 

The velocity components induced by the four 
corners of the panel and the wedce are now 
combined to obtain the resultant velocities 
at the control point. The velocity components 
induced by the right and left wing panels are 
combined and the results transformed back to 
the reference coordinate system by subroutine 
TRANS. This subroutine calculates the u, v, 
and w velocity components and the normal veloc- 
ity at the control point. A similar procedure 
is applied to calculate the transformed velocity 
components induced by the three components of 
the wake. The wake velocity components are then 
multiplied by the appropriate strength factors 
and added to obtain the net contribution of the 
wake. The wake velocities are then added to 
the panel velocities to obtain the final values 
of the velocity components at the control point. 

Special rules are applied to obtain the velocity 
components of the leading and trailing edge 
panels in each coliamn. These rules are designed 
to provide a continuous vortex distribution 
around the nose of the airfoil, and to enforce 
the Kutta condition at the trailing edge. 

The procedure is repeated for each column of 
panels of each wing segment. When all panel 
influences have been computed, the u, v, and w 
components of velocity are written as a single 
record on TAPE 8, and the normal velocities 
written in one array on TAPE 9. If the control 
point is in the same column of panels on the 
wing as the influencing panel, and the wing has 
more than 60 panels, the normal velocity at the 
control point is written on TAPE 10, and its 


USE: 


value set equal to zero in the array written 
on TAPE 9. This procedure sets up the diagonal 
blocks of the aerodynamic matrix for later use 
in the iterative solution procedure. Finally, 
if the print option is selected, the axial and 
normal velocity component arrays are written 
on the output file. 

This process is repeated for each control point. 
CALL OVERLAY (LWB, 2, 3) 


Input : 

Note: The word wing includes any tail, fin, 

or canard in the following descriptions . 


MACH 

PRINT 

NPART 

NMAX 

NWING 

NPOINT 

NSEG 

NROW 

NCOL 

NWT 

XPT, 

YPT, 

ZPT 

THET, 

THETI 

DELTA , 
DELTI 

XC, 

YC 


Mach number 

Print option parameter 

Matrix partition number 

Maximum order of diagonal block matrices 

Number of wing panels 

Number of control points 

Number of wing segments 

Number of rows of panels in segment 

Number of columns of panels in segment 

Tail segment identification parameter 

Arrays of control point coordinates 

Array of panel inclination angles 

Array of panel incidence angles 

Arrays of x and y coordinates of wing 
panel corner points 
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ZC Array of z coordinates of lower surface 

wing panel corner points 

ZU Array of z coordinates of upper surface 

wing panel corner points 

XS , Arrays of coordinates of point source 

YS, origins 

ZS 

Output ; 

I Control point index 

J, Wing panel index 

JJ 

L Panel row index 

N Panel column index 

NSIDE Column upper and lower surface index 

NS Wing segment number 

BETA Mach number parameter 

SUB Subsonic flow parameter (logical) 

SGN Supersonic flow sign parameter 

CON, Vortex panel constants 

BCON 

NR Number of rows of panels in segment 

NRl NR + 1 

NR2 , 2NR 

NRS 

NC Number of columns of panels in segment 

NCI NC + 1 

NT Tail segment identification parameter 

NI Number of first column in segment 

N2 Number of last column in segment 
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JL, 

Jl, 

JSl 

JT 

J2 , 

JS2 

11 

12 

JK, 

JM 

M 

K 

BD, 

TANBD , 
TAND 

SINBD, 

SIND, 

COSBD, 

COSD 

SINTI 

COSTI 

THETA 

SINT 

COST 

COSTD 

CONTD 

COSTD 

CONTDD 


Number of first vortex distribution on 
upper surface of column 

Number of last vortex distribution on 
upper surface of column 

Number of last vortex distribution on 
lower surface of column 

Number of last panel on upper surface 
of column 

Number of last panel on lower surface 
of column 

Temporary panel indices 

Panel leading or trailing edge index 
Panel side edge index 

Tangent of transformed panel incidence 
angle , 6 tan 6 

Trigonometric functions of transformed 
panel incidence angle 

sin 9 (I) 
cos 9 (I) 

Inclination angle of panel J 
sin 9 (J) 
cos 9 (J) 

cos 9(J)/(1. + (3 tan 
[(3 tan 6)^ + cos^9(J)]^ 
1./(C0SD*C0NTD) 

1 . /CONTD 



Coordinates of control point I 


XI, 

YI, 

ZI 

DXC , Differences between panel corner points 

DYC , in reference coordinate system 

DZC 

DXL, Differences between panel corner points 

DYL, in panel coordinate system 

DZL 

BL Panel edge sweep 

BLE Panel leading edge sweep 

BTE Panel trailing edge sweep 

AL, Difference between panel leading and 

A trailing edge sweeps 

CL Panel edge chord 

Cl Chord of inboard edge 

CO Chord of outboard edge 

DX, Control point coordinates with reference 

DY, to panel corner point 
DZ 

XJ, Control point coordinates in panel 

YJ, coordinate system 
ZJ 

X Dummy variable 

UCIR, Velocity components induced by inboard 
VCIR, leading edge corner of right wing 
WCIR panels containing constant vortex 

distribution 

ULIR, Velocity components induced by inboard 
VLIR, leading edge corner of right wing 
WLIR panels containing linearly varying 

vortex distribution 

RCIR, Same as UCIR, VCIR, WCIR, for inboard 

SCIR, trailing edge corner of right wing 

TCIR panel 
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RLIR, 

SLIR, 

TLIR 

UCIL, 

VCIL, 

WCIL 

ULIL, 

VLIL, 

WLIL 

RCIL, 

SCIL, 

TCIL 

RLIL, 

SLIL, 

TLIL 

VEIR, 

WEIR 


SEIR, 

TEIR 

VEIL, 

WEIL 

SEIL, 

TEIL 

VEOR, 
WE OR 


SEOR, 

TEOR 


Same as ULIR, VLIR, WLIR for inboard 
trailing edge corner of right wing 
panel 

Same as UCIR, VCIR, WCIR for left 
wing panels 

Same as ULIR, VLIR, WLIR for left 
wing panels 

Same as RCIR, SCIR, TCIR for left 
wing panels 

Same as RLIR, SLIR, TLIR for left 
wing panels 

Velocity components induced by concen- 
trated vortex from leading edge along 
inboard edge of right wing panel 

Same as VEIR, WEIR for vortex from 
trailing edge 

Same as VEIR, WEIR for left wing panel 
Same as SEIR, TEIR for left wing panel 


Velocity components induced by concen- 
trated vortex from leading edge along 
outboard edge of right wing panel 

Same as VEOR, WEOR for vortex from 
trailing edge 


VEOL, Same as VEOR, WEOR for left wing panel 
WEOL 


SEOL, Same as SEOR, TEOR for left wing panel 
TEOL 


VAIR, Velocity components induced by vortex 
WAIR sheet from inboard leading edge corner 
of right wing panel 


127 



SAIR, Same as VAIR, WAIR for vortex sheet 
TAIR from trailing edge 

VAIL, Same as VAIR, WAIR for left wing panel 
WAIL 

SAIL, Same as SAIR, TAIR for left wing panel 
TAIL 

VAOR, Same as VAIR, WAIR for outboard corner 
WAOR of right wing panel 

SAOR, Same as SAIR, TAIR for outboard corner 
TAOR of right wing panel 

VAOL, Same as VAOR, WAOR for left wing panel 
WAOL 

SAOL, Same as SAOR, TAOR for left wing panel 
TAOL 

UAR, Velocity components induced by vortex 

VAR, sheet behind right wing panels 

WAR 

UAL, Same as above for left wing panels 
VAL, 

WAL 

UIR, Velocity components induced by inboard 
VIR, concentrated vortex behind right wing 
WIR panels 

UIL, Same as above for left wing panels 

VIL, 

WIL 

UOR, Velocity components induced by outboard 
VOR, concentrated vortex behind right wing 
WOR panels 

UOL, Same as above for left wing panels 

VOL, 

WOL 

ULR, Velocity components induced by linearly 

VLR, varying vortex distribution having zero 

WLR strength along leading edge on right 

wing panels 


SUBROUTINES 
CALLED : 

ERROR 
RETURNS : 


ULL, Same as above for left wing panels 
VLL, 

WLL 

UCR, Velocity components induced by linearly 

VCR, varying vortex distribution having zero 

WCR strength along trailing edge on right 

wing panels 

UCL, Same as above for left wing panels 
VCL, 

WCL 

UC , Arrays of velocity components induced 

VC, by vortex panels at control point I 

WC 

US, Velocity components induced by point 

VS, sources at control point I 

WS 

AC Array of normal velocities induced by 

vortex panels at control point I 

AS Normal velocity induced by point sources 

at control point I 

DC Array of normal velocities induced by 

vortex panels in diagonal block matrices 


VORPAN, TRANS 


Program calls EXIT if 3 tan 6 > 1. in supersonic 
flow 
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APPENDIX II 


PROGRAM LISTING 
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160 FORMAT tIOFT.O) M1290 

170 FORMAT CIH »39HERR0R - BODY HAS MORE THAN 20 MERIDIANS) M1300 

180 FORMAT (IH tAAHERROR - BODY HAS MORE THAN 60 AXIAL STATIONS) M1310 

END M1320 
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REVISE CHOROWISE PANEL SPACING ON FIN (VERTICAL TAIL) OR CANARD 
(HORIZONTAL TAiL) AND COMPUTE NEW AIRFOIL ORDINATES 
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PROGRAM VELCMP 

COMPUTE THE VELOCITY CCMPONENTS (UtV*VIJ AT THE PANEL CONTROL 
POINTS AND FORM THE AERODYNAMIC INFLUENCE COEFFICIENT MATRICES 
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COMPOTE THE THREE COMPONENTS OF VELOCITY INDUCED AT SPECIFIED U 40 
CONTROL POINTS BY SOURCE AND VORTEX DISTRIBUTIONS ON PANELS U SO 
LOCATED IN THE PLANE OF THE WINGt FIN (VERTICAL TAIL) t OR CANARD U 60 
(HORIZONTAL TAiL) SURFACE U 70 
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COMPUTE THE THREE COMPONENTS OF VELOCITY INDUCED AT A SPECIFIED 
CONTROL POINT BY CONSTANT AND LINEARLY VARYING VORTEX 
DISTRIBUTIONS ON A SWEPT QUAOR-ILATERAL PANEL. VORVEL CALCULATES 
THE VELOCITY INCUCED BY THE LEADING AND TRAILING CORNERS OF ONE 
EDGE OF THE PANEL. 
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IF (C.EQ.O.I GO TO 330 Y3760 
RETURN Y3770 
END Y3780 
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APPENDIX III 


SAMPLE CASE 



UNIFIED SUBSONIC-SUPERSONIC AEROOYNANICS PROGRAM 


VERSICN AGO 


LIST OF INPUT CAROS 


0000000001 ill 11 ilXl2222222Z2233333J3333AA4A4AA44A55555S5555666666666677777 7777 78 
1234S678901234567890123456789012345678901234S676901234367890123A5678901234567890 


OGIVE CVLIKOER BODY M|TH 45 DEGREE SmEEP NACA 65A004 HlO-kiING 
0 -I -1 -1 2 26 1 5 28 
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100. 
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10. 
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27.65 12. 
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2. 
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0. .3075 
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1.6505 
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1.657 

1. 728 
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1.3815 
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1.174 .949 

.715 
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0. .08605 

.32573 

.69221 

1. 160951. 709202.315422. 960603.627134.29759 

FLSARDl 

4.95834 5. 595676. 197 1 76. 75269 7. 25261 7.689558.055558. 3458 78.556108.68411 

FUSAR02 

8.72 7 8.727 

8.727 

8.727 

8.727 

8. 72 7 

i.727 

8. 727 



FUSAR03 

SINGULARITY PANELING 1 

FOR SAMPLE CASE 






1 1 -2 










1 3 1 
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UGIVE CYLINDER BODY UITH 45 DEGREE SrEEP NACA 65A004 MIO-RlNG 



MING PANEL COHNER POINT COORDINATES 
I ANO 3 INDICATE MiNG PANEL LEADING-EDGE POINTSt 


2 ANO h INDICATE TRAlLlNGr-EOGE POINTS 


PAN EL 
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1.66700 
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1.66700 
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7 
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8 
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10 
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11 
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12 
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13 
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14 
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16 
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1 7 
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19 
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20 
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2 L 
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22 
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2 3 
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24 
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25 
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27 
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35 

24.60700 

7. 73000 
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43 
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2.97000 
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20.92803 
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0.00000 
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0.00000 
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0.00000 
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0.00000 
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0.00000 
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0. 00000 
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0. ooooo 
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0.00000 
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0. ooooo 
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25.13500 
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0. ooooo 

25.69300 
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0.000 00 

20.55700 
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0.00000 

22.66833 

7.73000 

0.00000 
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0.00000 
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7.73000 

0.00000 
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0.00000 
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0. 00000 
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0.00000 

24.12233 
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0.00000 
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0. ooooo 

24.60700 
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0. 00000 
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0. ooooo 
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7.73000 

0.00000 
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5.37000 

0. ooooo 
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7.73000 

0.00000 

25.05100 
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0.00000 

26.06100 

7.73000 

0.00000 

25.69300 
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0.00000 
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7.73000 

0.00000 
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0.00000 
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0.00000 
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0.00000 
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0.00000 

25.09167 

7.73000 

0.00000 
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10.10000 

0.00000 

25.57633 

7.73C00 

c. ooooo 

27.06667 

10.10000 

0.00000 

26.06100 

7.73000 

0.00000 

27.3933^ 

10.10000 

0.00000 

26.54567 

7. 73000 

0.00000 

27.72000 

10.10000 

0.00000 

27.03033 

7. 73GOO 

0.00000 

28.04667 

10.10000 

0.00000 

27.51500 

7. 73000 

0.00000 

28.37333 

10.10000 

0.00000 

25.76000 

1 0. 10000 

0. ooooo 

27.65000 

12.00000 


0,00000 26.0db67 lO. 10000 0.00000 27.«5000 12.00000 
0.00000 26.41333 10.10000 0.00000 28.0S000 12.000Q0 
0.00000 26,74000 10.10000 0.00000 28.25000 12.00000 
0.00000 27.06667 10.10000 0.00000 28.45000 12.00000 
0.00000 27.39333 10.10000 0.00000 28.65000 12.00000 
0-00000 27.72000 10.10000 0.00000 28.85000 12.00000 
0.00000 28.04667 10.10000 0.00000 29.05000 12.00000 
0.00000 28.37333 10.10000 0.00000 29.25000 12.00000 
0.00000 28.70000 10.10000 0.00000 29.45000 12.00000 


z 

3 


0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0,00000 
0.00000 
o.poooo 
o.oopoo 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0,00000 
0,00000 
0.00000 
0.00030 
0.00000 
0. 00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0,00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 


X 

4 


1^.91700 
18.71900 
19.52100 
20.32300 
21.12500 
21.92700 
22.72900 
23.53100 
24.33300 
25,13500 
20.55700 
21.19900 
21.84100 
22.48300 
23^12500 
23f 76700 
24.40900 
25.05100 
25.69300 
26.33500 
23.15300 
23.63767 
24.12233 
24.«60 700 
25.09167 
25.57633 
26.06100 
26.54567 
27.03033 
27.51^00 
25.76000 
26.08667 
26.41333 
26.74000 
27.06667 
27.39333 
27.72000 
28.64667 
28.37333 
26.70000 

27.85000 

28.05000 

28.25000 

28.45000 

28.65000 

28.85000 

29.05000 

29.25000 

29.45000 

29.65000 


Y 

4 

2.97000 
2.97000 
2.97000 
2.97000 
2.97000 
2.97000 
2.97000 
2.97000 
2.97000 
2.97000 
5.37000 
5.37000 
5.37000 
5.37000 
5.37000 
5.37000 
5.37000 
5.37000 
5.37000 
5.37000 
7.73000 
7.73000 
7.73000 
7. 73000 
7.73000 
7.73000 
7.73000 
7.73000 
7.73000 
7.73000 
10.10000 
10.10000 
10.10000 
10.10000 
10.10000 
10.10000 
10.10000 
10.10000 
10.10000 
10.10000 
12.00000 
12.00000 
12.00000 
12.00000 
12.00000 
12.00000 
12.00000 
12.00000 
12.00000 
12.00000 


I 

4 


0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
o.ooobo 
o.ooo'oo 
0.00000 
0.00000 
0. 00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0 . 0000 0 
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WING 

PANEL CONTkOL POINTS ANO 

POINT 

X y 


CP CP 


1 

16.34190 

2.30734 

2 

17.18808 

2.30734 

3 

18.03425 

2.30734 

4 

18.88043 

2.30734 

5 

19.72661 

2.30734 

6 

20.57279 

2.30734 

> 

21.41896 

2.30734 

8 

22.26514 

2.30734 

9 

23.11132 

2.3C734 

10 

23.95749 

2.30734 

11 

24.80367 

2.30734 

12 

18.46329 

4.12568 

U 

19. 18825 

4.12568 

14 

19.91320 

4.12566 

15 

20.63816 

4.12568 

16 

21.36311 

4.12568 

17 

22.08807 

4.12568 

18 

22.81302 

4.12568 

19 

23.53798 

4.12568 

20 

24.26293 

4.12568 

21 

24.98788 

4.12568 

22 

25.71284 

4.12568 

23 

21.22759 

6.49507 

24 

21.79458 

6.49507 

25 

22.36158 

6.49507 

26 

22. 92857 

6.49507 

27 

23.49557 

6.49507 

28 

24.06256 

6.49507 

29 

24. 62956 

6.49507 

30 

25.19655 

6 . 4950 7 

31 

25.76355 

6.49507 

32 

26.33054 

6.49507 

33 

26.89754 

6.49507 

34 

23.96109 

8.83808 

35 

24.37188 

8.83808 

36 

24.78268 

6.63808 

37 

25.19347 

8.63608 

38 

25.60427 

8.63808 

39 

26.01506 

8.83808 

40 

26.42586 

8.83808 

41 

26.83665 

8.83808 

42 

27.24745 

8.83808 

43 

27.65824 

8.83808 

44 

28.06904 

8.83808 

45 

26,45281 

10.97384 

46 

26. 72122 

10.97384 

47 

26.98963 

10.97384 

48 

27.25805 

10.97384 

49 

27.52646 

10.97384 

50 

27.79487 

10.97384 

51 

28. 06328 

10.97384 

52 

28.33169 

10.97384 

53 

28.60010 

10.97384 

54 

28.86851 

10.97384 

55 

29.13692 

10.97384 


INCLINATION ANGLES 


1 

THETA 

CAH8EK 

CP 


SLOPE 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

o.coooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

o.coooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

o.coooo 

0. 00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

o.coooo 

0.00000 

0.00000 

0. ooooo 

0.00000 

0.00000 

0.00000 

0. 00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0. 00000 

0.00000 

0.00000 

0,00000 

0.00000 

0.00000 

0.00000 

0.00000 

0. ooooo 

0.00000 

0.00000 

o.coooo 

0.00000 

0.00000 

0.00000 

0. 00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

o.coooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0,00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0. 00000 

0,00000 

0.00000 

0.00000 

0.00000 

0. ooooo 

0. 00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0. ooooo 

0.00000 

0.00000 

0.00000 

0.00000 

o.coooo 

0. ooooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

o.coooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

o.coooo 

0.00000 

0.00000 

o.coooo 

0.00000 

0.00000 

o.coooo 

o.coooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

o.coooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0. coooo 

0.00000 

o.coooo 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

o.coooo 

0.00000 

0.00000 

0.00000 


THICKNESS 

SLOPE 


• 178 79 
•05696 
.03265 
•01709 
•00315 
-.01379 
-.02918 
-•03946 
-.04623 
-.04704 
-.04710 
.17879 
.05696 
•03265 
•01709 
•00315 
-.01379 
-.02918 
-.03946 
-.04623 
-.04704 
-.04710 
• 17879 
.05696 
.03265 
.01709 
.00315 
-.01379 
-.02918 
-.03946 
-.04623 
-.04704 
-.04710 
.17879 
.05696 
•03265 
.01709 
•00315 
-.01379 
-.02918 
-.03946 
-.04623 
-.04704 
-.04710 
.17879 
.05696 
.03265 
•01709 
.00315 
-.01379 
-.02918 
-•03946 
-.04623 
-.04704 
-.04710 



iilHlf PANEL AREAS AND CriQHOS 


PANEL 

AREA 

CHORD 

1 

1.10160 

.64618 

2 

1. 10160 

.84618 

3 

1.10160 

. £4618 

4 

1. 10160 

.84618 

5 

1.10160 

•84618 

6 

1*10160 

.84618 

7 

1.10160 

.84618 

8 

1.10160 

.£4618 

9 

1. 10160 

.84618 

10 

I. 10160 

•64613 

U 

1.73280 

.72495 

12 

1. 73280 

.72495 

13 

1. 73280 

.72495 

19 

1.73280 

.72495 

15 

1. 73280 

.72495 

16 

1. 73280 

.72495 

17 

1.73280 

.72495 

18 

1. 73280 

. 7 249 5 

19 

1.73280 

.72495 

20 

1. 73280 

.72495 

21 

1. 32947 

•56700 

22 

1. 32947 

.56700 

23 

1.32947 

.56700 

24 

1.32947 

.56700 

25 

1. 32947 

.56700 

26 

1.32947 

•56700 

27 

1.32947 

.56700 

28 

1.32947 

.56700 

29 

1.32947 

.56700 

30 

1.32947 

.56700 

31 

.96143 

.41079 

32 

.96143 

.41079 

33 

.96143 

.41079 

34 

.96143 

.41079 

35 

,96143 

.41079 

36 

.96143 

.41079 

37 

•96143 

.41079 

38 

.96143 

•41C79 

39 

.96143 

.41079 

40 

.96143 

.41079 

41 

.50033 

.26841 

42 

.50033 

.26841 

43 

.50033 

•2 6841 

44 

.50033 

.26841 

45 

.50033 

.26641 

46 

. 50033 

•26841 

47 

.50033 

.26841 

48 

.50033 

.26841 

49 

.50033 

.26841 

50 

.50033 

.26341 
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PANEL 


BODY PANEL CORNER POINT COORDINATES 

I *NO 3 .NDICAIE 800r PANEL LEADING-EDGE POINTS. E AND A INDICATE TRAILING-EDGE POINTS 


1 

2 

3 

4 

5 

6 
7 

a 

9 

10 

11 

12 

13 

14 

15 

16 
17 

la 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 
56 

59 

60 


0.00000 

> 0.00000 

0.00000 

0.00000 

0. OOOOO 

0. OOOOO 

0.00000 

0.00000 

1. 50000 

.COOOO 

1.50000 

.28724 

1.50000 

.40622 

1.50000 

.26724 

4. 50000 

.00000 

4.50000 

. 73883 

4.50000 

1.04487 

4.50000 

• 73883 

7.50000 

•OOOOO 

7. 50000 

1.03047 

7.50000 

1.45730 

7.50000 

1.03047 

10.50000 

•OOOOO 

10. 50000 

1. 16694 

10.50000 

1.65030 

10.50000 

1.16694 

11.66700 

•OOOOO 

11.66700 

1.17854 - 

11.66700 

1.66670 

11.66700 

1.17854 

15. 59480 

•COOOO • 

15.59480 

1.17854 - 

15. 59480 

1 . 66 6 70 

15.59460 

1. 17854 

17.3 726 0 

.00000 - 

17.37260 

1.17854 - 

17.37260 

1.66670 

17.37260 

1.17854 

19. 15030 

.00000 - 

19.15030 

1.17854 - 

19.15030 

1.66670 

19.15030 

1.17854 

20.92800 

.00000 - 

20.92800 

1.17854 - 

20.92800 

1.66670 

20.9280 0 

1. 17854 

22. 70580 

•OOOOO - 

22.70580 

1.1/854 - 

22,70580 

1.66670 

22. 70580 

1.17854 

24. 48350 

•COOOO - 

24.48350 

1.17854 - 

24.48350 

1.66670 

24.48350 

1.17854 

26.^8000 

•OOOOO 

26.28000 

1. 17854 - 

26.28000 

1.66670 

26.28000 

1.17634 J 

29.40000 

.00000 - 

29.40000 

1.17854 - 

29.40000 

1.66670 

29.40000 

1. 17634 

33.00000 

•OOOOO - 

33.00000 

1.1/854 - 

33.00000 

L.66670 

33.00000 

1.17854 1 


0.00000 
0.00000 
0.00000 
0.00000 
-.40622 
-.28724 
•00000 
.28724 
-1.04487 
-.73683 
.00000 
• 73883 
-1.45730 
-1.03047 
.00000 
1.03047 
-1.65030 
-1.16694 
.00000 
1. 16694 
-1.66670 
-1.17854 
. 00000 
1.17354 
-1.66670 
-1. 17854 
. OOOOO 
1.17854 
-1.66670 
-1. 17854 
.00000 
1. 17854 
-1.66670 
-1.17854 
.00000 
1. 17654 
-1.66670 
-1. 17854 
.00000 
I. 17654 
1.66670 
1.17854 
.00000 
1.17SS4 
1.66670 
1.17854 
.00000 
1. 17854 
1.66670 
1. 17854 
.00000 


1.66670 
1.17854 
.00000 
. 17854 
1.66670 
.17854 
.00000 
. 17854 


1.50000 

1.50000 

1.50000 

1.50000 

4.50000 
4.50000 
4.50000 

4.50000 

7.50000 
7.50000 
7.50000 
7.50000 

10.50000 
10.50000 
10.50000 
10.50000 
11.66700 
11.66700 
11.66700 
11.66700 
15.5 9480 
15.59480 
15.59480 
15.59480 
17.37260 
17.37260 
17.37260 
17.37260 
19.15030 

19. 15030 

15. 15030 

19.15030 
20.92800 
20.92800 
20.92800 
20.92 800 
22. 70580 
22.70580 
22. 70580 
22.70580 
24.48350 
24.48350 
24.48350 
24.48350 
26.28000 
26.28000 
26.28000 
26.28000 
29.40000 
29.40000 
29.40000 
29.40000 
33.00000 
33.00000 
33.00000 
33.00000 
36.50000 
36.50000 
36.50000 
36.50000 


• COOOO 
*26 724 
.40622 
. 28 724 
•OOOOO 
.73863 
1.04487 
. 73883 
. OOOOO 
1.03047 
1.45730 
1.03047 
.00000 
1. 16694 
1.65030 
1* 16694 
.OOGOO 
1. 17854 
1.66670 
1.17854 
•OOOOO 
1.17854 
1.66670 
1. 17854 
.00000 
1. 17854 
1.66670 
1.17854 
.00000 - 
1. 17854 - 
1.66670 
1.17854 
•COOOO - 
1.17854 < 
1.66670 
1.17854 
.COOOO - 
1.17854 - 
1. 66670 
1. 17854 
.00000 - 
1.17854 - 
1.66670 
1. 17854 
.COOOO 
1. 17854 
1.66670 
1. 17854 
.00000 
1.17854 
1.66670 
1. 17854 
.00000 
1.17854 
1.66670 
1.17854 
.00000 
1. 17854 
1. 66670 
1.17854 


-.40622 
-.28724 
.00000 
.28724 
-1.04487 
-. 73883 
•OOOOO 
• 73883 
-1.45730 
•1.03047 
•OOOOO 
1.03047 
1.650 30 
1. 16694 
.00000 
1. 16694 
1.66670 
1. 17854 
.00000 
1. 17854 
-1.66670 
-1. 17854 
•OOOOO 
1. 17854 
-1.66670 
-1. 17854 
•OOOOO 
1. 17854 
-1.66670 
-1.17854 
.00000 
1. 17854 
-1.66670 
-1. 17854 
•OOOOO 
1. 17854 
-1.66670 
-1. 17854 
.00000 
1. 17854 
-1.66670 
-1. 17854 
.00000 
1. 17854 
-1.66670 
-1. 17854 
•OOOOO 
1. 17854 
-1.66670 
-1. 17854 
•OOOOO 
1.17854 
-1.66670 
-1. 17854 
.00000 
1. 17854 
-1.66670 
1. 17854 
•OOOOO 
1. 17854 


0.00000 
0.00000 
0.00000 
0.00000 
1.50000 
1.5COOO 
1.50000 

1.50000 

4. 50000 
4.50000 
4.50000 

4.50000 

7.50000 
7.50000 
7. 50000 
7.50000 

10. 50000 
10.50000 
10.50000 
10.50000 
11.66700 
11.66700 
11.66700 
11.66700 
15.59480 
15. 59480 
15.59480 
15. 59480 
17.37260 
17.37260 
17.37260 
17.37260 
19.15030 
19.15030 
19. 15030 
19. 15030 
20. 92800 
20.92600 
20.92800 
20. 92800 
22. 70580 
22.70580 
22. 70580 
22.70580 
24.48350 
24.483 50 
24.48350 
24. 48350 
26.28000 
26.28000 
26.28000 
26.28000 
29.40000 
29.40000 
29.40000 
29.40000 
33.00000 
33.00000 
33.00000 
33.00000 


0.00000 
0.00000 
0.00000 
O.OQOOO 
.28724 
.40622 
.28724 
-•OOOOO 
.73883 
1.04487 
. 73883 
-.00000 
1.03047 
1.45730 
1.03047 
-.00000 
1.16694 
1.65030 
1.16694 
-.00000 
1.17854 
1.66670 
1.17854 
-.00000 
1.17854 
1.66670 
1.17854 
-.00000 
1.17854 - 
1.66670 
1.17854 
-.00000 
1.17854 - 
1.66670 
1.17854 
-•OOOOO 
1.17654 - 
1.66670 
1.17854 
-.00000 
1.17854 - 
1.666T0 
1.17854 
-.00000 
1.17654 
1.66670 
1.17854 
-•OOOOO 
1.17854 
1.66670 
1.17854 
-.00000 
1.17854 
1.66670 
1.17854 
-.00000 
1.17854 
1.66670 
1.17854 
-.00000 


0.00000 
0.00000 
0.00000 
0.00000 
-.28724 
.00000 
.28724 
•40622 
-.73883 
.00000 
.73883 
1.04487 
1.03047 
.00000 
1.03047 
1.45730 
-1.16694 
•OOOOO 
1.16694 
1.65030 
-1.17854 
.00000 
1.17854 
1.66670 
-1.17854 
•OOOOO 
1.17854 
1.66670 
-1. 17854 
•OOOOO 
1. 17854 
1.66670 
-1. 17854 
•OOOOO 
1.17854 
1.66670 
-1. 17854 
•OOOOO 
1.17854 
1.66670 
-1.17854 
.00000 
1. 17854 
1.66670 
-1.17854 
.00000 
1. 17854 
1.66670 
-1.17854 
•OOOOO 
1.17854 
1.66670 
-1. 17854 
• OOOOO 
1.17854 
1.66670 
1. 17854 
• OOOOO 
1.17854 
1.66670 


1.50000 

1.50000 

1.50000 

1.50000 

4.50000 
4.50000 
4.50000 

4.50000 

7.50000 
7.50000 
7.50000 
7.50000 

10.50000 
10.50000 
10.50000 
10. 50000 
11.66700 
11.66700 
11.66700 
11.66700 
15.59480 
15.59480 
15.59460 
15.59480 
17.37260 
17.37260 
17.37260 
17.37260 
19. 15030 
19.15030 
19.15030 
19.15030 
20.92000 
20.92800 
20.92800 
20. 92800 
22.70580 
22.70580 
22.70580 
22.70500 
24.48350 
24.48350 
24.46350 
24.48350 
26.28000 
26.28000 
26.28000 
26.28000 
29.40000 
29.40000 
29.40000 
29.40000 
33.00000 
33.00000 
33.00000 
33.00000 
36.50000 
36.50000 
36.50000 
36.50000 


.28724 
•40622 
.28724 
-.00000 
.73683 
1.04487 
.73863 
-.00000 
1.03047 
1.45730 
1.03047 
-.00000 
1.16694 • 
1.65030 
1.16694 
-.00000 
1.17854 - 
1.666 70 
1.17854 
-.00000 
1.17854 - 
1.66670 
1.17854 
-•OOOOO 
1.17854 - 
1.66670 
1. 17854 
-.00000 
1.17854 - 
1.66670 
1.17854 
-.00000 
1.17854 - 
1.66670 
1.17854 
-.00000 
1.17854 
1.66670 
1.17854 ] 

-.00000 1 
1.17854 -j 
1.66670 
1-17854 ] 

-.00000 ] 
1.17854 
1 . 66 670 
1. 17654 1 

-.00000 1 
1.17854 -1 
1.66670 
1.17854 1 

-.00000 1 
1.17854 >1 
1.66670 
1.17854 1 

-.00000 1 
1.17854 -1 
1.66670 
1.17854 1 

-.00000 1 


I 

4 


-.28724 
.00000 
.28724 
.40622 
-.7388 3 
•OOOOO 
. 73883 
1.0448 7 
-1.0304 7 
.00000 
1.03047 
1.45730 
-1. 16 69 4 
.00000 
1.16694 
1.65030 
-1.17854 
•OOOOO 
1.17854 
1.66670 
1.17854 
.00000 
1.17854 
1.66670 
1. 17854 
.OOOOO 
1. 17854 
1.66670 
1.17854 
.00000 
1. 17854 
1.66670 
1.17854 
.00000 
1.17854 
1.66670 
L.17854 
.00000 
L.17854 
L. 66670 
L. 17854 
-OOOOO 
1.17854 
. .66670 
.17854 
.00000 
.17854 
•66670 
. 17854 
.00000 
.17854 
.6667 0 
.17854 
.00000 
.17854 
.66670 
.17854 
•OOOOO 
.17854 
.66670 
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bOCy PANEL CONTROL POINT COORDINATES 


POINT 

X 

CP 

i 

1.00000 

2 

1.00000 

3 

1.00000 

A 

1.00000 

b 

3.22006 

b 

3.22006 

7 

3.22006 

8 

3.22006 

9 

6.06242 

iO 

6.08242 

u 

6.08242 

12 

6.08242 

13 

9.03105 

14 

9.03105 

15 

9.03105 

lb 

9.03105 

17 

1 L. 06446 

Itt 

11.06446 

19 

11.08446 

20 

11.06446 

21 

13.63090 

22 

13.63090 

23 

13.63090 

24 

13.63090 

25 

16.46370 

2b 

16.48370 

27 

16.46370 

26 

16.48370 

29 

18.26145 

30 

18.26145 

31 

18. 26145 

32 

18.26145 

33 

20.03915 

34 

20.03915 

35 

20.03915 

3b 

20.03915 

37 

21.81690 

36 

21.81690 

39 

21.81690 

40 

21. 81690 

41 

23. 59465 

42 

23.59465 

43 

23.59465 

44 

23. 59465 

45 

25.38175 

46 

25.38175 

47 

25.36175 

46 

25.38175 

49 

27.64000 

50 

27.84000 

51 

27.84000 

52 

27. 84000 

53 

31.20000 

54 

31. 20000 

55 

31.20000 

56 

31.20000 

57 

34.75000 

58 

34.75000 

59 

34.75000 

60 

34. 75000 


CP 

CP 

.09575 

-.23116 

• 23116 

-.09575 

.23116 

.09575 

.09575 

.23116 

.27308 

-.65928 

.65928 

-.27308 

.65928 

.27308 

.27308 

.65928 

•44633 

-1.07754 

1.C7754 

-.44633 

1.07754 

,44633 

.44633 

1.07754 

•55006 

-1.32796 

1.32796 

-.55006 

1.32796 

•55006 

.55006 

1.32796 

.58637 

-1.41563 

1.41563 

-.58637 

1.41563 

.58637 

•58637 

1.41563 

.58927 

-1.42262 

1.42262 

-.58927 

1.42262 

.58927 

.58927 

1.42262 

.58927 

-1.42262 

1.42262 

-.58927 

1.42262 

.58927 

.58927 

1.42262 

• 58927 

-1.42262 

1.42262 

-.58927 

1.42262 

.58927 

.58927 

1.42262 

.58927 

-1.42262 

1.42262 

-.58927 

1.42262 

.58927 

.58927 

1.42262 

•58927 

-1.42262 

1.42262 

-.58927 

1.42262 

.58927 

.58927 

1.42262 

.58927 

-1.42262 

1.42262 

-.50927 

1.42262 

.58927 

.58927 

1.42262 

•58927 

-1.42262 

1.42262 

-.58927 

1.42262 

•58927 

.50921 

1.42262 

.58927 

-1.42262 

1.42262 

-.58927 

1.42262 

.58927 

.58927 

1.42262 

.58927 

-1.42262 

1.42262 

-.58927 

1.42262 

.58927 

•58927 

1.42262 

.50927 

-1.42262 

1.42262 

-.58927 

1.42262 

.58927 

.50927 

1.42262 
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BODY PANEL AREAS AND INCLINATION ANGLES 


PANEL AREA DELTA THEIA 


1 

•24037 

2 

.24037 

3 

•24037 

4 

.24037 

5 

1.69784 

b 

1.69784 

7 

1.69784 

B 

I. 6S784 

9 

2.89570 

iO 

2.89570 

U 

2.89570 

12 

2.89570 

13 

3.57398 

14 

3. 57398 

15 

3. 57398 

16 

3.57398 

17 

1.4B147 

IB 

1.48147 

19 

1. 48147 

20 

1. 48147 

21 

5.01045 

22 

5.01045 

23 

5.01045 

24 

5.01045 

25 

2.26783 

26 

2.26783 

27 

2.26783 

26 

2. 26783 

29 

2*26770 

30 

2.26770 

31 

2.26770 

32 

2,26770 

33 

2.26770 

34 

2.26770 

35 

2.26770 

36 

2.26770 

37 

2. 26783 

36 

2.26763 

39 

2.26783 

40 

2.26783 

41 

2.26770 

42 

2.26770 

43 

2.26770 

44 

2.26770 

45 

2.29168 

46 

2.29168 

47 

2.29168 

4B 

2.29168 

49 

3. 97999 

50 

3.97999 

51 

3. 97999 

52 

3.97999 

53 

4. 59229 

54 

4.59229 

55 

4.59229 

56 

4.59229 

57 

4.46473 

58 

4. 464 73 

59 

4.464 73 

60 

4.46473 


.24517 

-2. 74889 

.2451? 

-I.S6350 

.24517 

-1.17810 

-24517 

-.39270 

.19420 

-2. 74889 

.19420 

-I.S6350 

.19420 

-1.17810 

• 19420 

-.39270 

.12634 

-2.74889 

.12634 

-1.96350 

• 12634 

-1.17810 

.12634 

-.39270 

.05937 

-2.74869 

.05937 

-1.96350 

.05937 

-1.17810 

.05937 

-.39270 

.01298 

-2. 74889 

.01298 

-1.96350 

.01298 

-1.17810 

.01298 

-.392 70 

0.00000 

-2.74889 

o.coooo 

-1.96350 

0.00000 

-1.17810 

o.coooo 

-.39270 

0.00000 

-2.74889 

0.00000 

-1.96350 

o.coooo 

-1. 17810 

o.coooo 

-.392 70 

0.00000 

-2.74689 

0.00000 

-1.96350 

0.00000 

-1.17810 

0.00000 

-.392 70 

0.00000 

-2. 74889 

0.00000 

-1.96350 

0,00000 

-1.17610 

o.coooo 

-.39270 

0.00000 

-2 . 74889 

0.00000 

-1.96350 

0.00000 

-1.17810 

o.coooo 

-.39270 

0.00000 

-2.74689 

0.00000 

-1.96350 

o.coooo 

-1. 17810 

O.QOOOO 

-.392 70 

0.00000 

-2.74889 

0.00000 

-1.96350 

0.00000 

-1.17810 

0.00000 

-.39270 

O.COOOO 

-2. 74889 

O.QOOOO 

-1 .96350 

O.OGOOO 

-1.17810 

0.00000 

-.39270 

o.coooo 

-2.74669 

O.COOOO 

-1.96350 

0.00000 

-1.17810 

0.00000 

-.39270 

0.00000 

-2.74889 

0.00000 

-1.96350 

0.00000 

-1. 17810 

0.00000 

-.39270 


PARTITION « 1 TINE « 86.6A100 

INFLUENCE OF BODY ON BODY 


PARTITION • 2 TIME « U6.B6700 

INFLUENCE OF WING ON BODY 


PARTITION « 3 TIME « 146.2B100 

INFLUENCE OF BODY ON WING 


PARTITION « A TIME • 173.92100 
INFLUENCE OF WING ON WING 

TIME * 204. 25900 

TIME • 206.69900 


TIME - 224.44500 


VEtOClTlES CN aOOYf HACH « 2 » 0i0 ALPHA * 


0.000 


PANEL 

SOURCE 

AAIAL 

LATERAL 

VERTICAL 

NORMAL 

NO . 

STRENGTH 

velocity 

VELOCITY 

VELOCITY 

VELOCITY 

I 

.196 42 

-.09678 

.08648 

-.20879 

•25020 

2 

. 196 42 

-.09678 

• 208 79 

-.08648 

•25020 

3 

.196 42 

-.09678 

• 208 T9 

•08648 

.25020 

4 

.19642 

-.09678 

•08648 

.20879 

.25020 

5 

.17261 

-.07016 

•06998 

-.16896 

• L9668 

6 

.17261 

-.07016 

. 16896 

-.06998 

•19668 

7 

.17261 

-• Q70I6 

. 16896 

•06998 

•19668 

8 

•17261 

-.07016 

. 06998 

.16896 

•19666 

9 

.11296 

-.03539 

.04689 

-.11319 

.12701 

iO 

.11296 

-.03539 

.11319 

-.04689 

•12701 

U 

.11296 

-.03539 

• 11319 

.04689 

•12701 

12 

.11296 

-.03539 

.04689 

•11319 

.12701 

13 

•05198 

-.00181 

.02270 

-.05481 

.05944 

14 

.05198 

-.00181 

.05481 

-.02270 

•05944 

15 

.05198 

-.00161 

.05481 

•02270 

.05944 

16 

.05198 

-.00181 

.022 70 

.05481 

.05944 

17 

-.02286 

.02050 

.00507 

-.01224 

•01298 

18 

-.02286 

.02050 

.01224 

-.00507 

.01298 

19 

-.02286 

.02050 

.01224 

.00507 

.01298 

20 

-.02286 

.02050 

.00507 

•01224 

.01298 

21 

-.00277 

.01807 

• 00000 

.00000 

-.00000 

22 

-.00277 

.01807 

.00000 

• 00000 

-.00000 

23 

-.00277 

.01807 

-• 00000 

.00000 

-.00000 

24 

-.00277 

.01807 

-. 00000 

•00000 

-.00000 

25 

.02891 

•00714 

. 00000 

•00000 

-.00000 

26 

.02891 

.00714 

-. 00000 

-.00000 

.00000 

27 

.02891 

.00714 

. 00000 

-.00000 

-.00000 

28 

.02891 

.00714 

-.00000 

•ooooo 

-.00000 

29 

-.00208 

.00751 

-.00000 

-.00000 

-.00000 

30 

.023 71 

-.00814 

-.00807 

-.01949 

-.00000 

31 

.02371 

-.00814 

-.00807 

•01949 

-.00000 

32 

-.00208 

.00751 

. 00000 

-.00000 

•ooooo 

33 

-.02323 

.00037 

-. 01946 

-.00806 

•OOOOO 

34 

.01089 

.00055 

-.00635 

-.01534 

-.00000 

35 

.01089 

.00055 

-• 006 35 

.01534 

-.00000 

36 

-.02323 

•00037 

-.01946 

.00806 

•ooooo 

37 

-.04014 

-.00101 

-.01001 

-.00415 

•ooooo 

38 

-.01999 

.01434 

.00025 

•00060 

•OOOOO 

39 

-.01999 

.01434 

.00025 

-.00060 

.00000 

40 

-.04014 

-.00101 

-.01001 

•00415 

•ooooo 

41 

.00854 

-.00006 

.00244 

•00101 

-•ooooo 

42 

-.00627 

.01918 

,00711 

.01716 

.00000 

43 

-.00627 

.01918 

.00711 

-.01716 

.00000 

44 

.00854 

-.00006 

. 00244 

-.00101 

.00000 

45 

.03622 

•00587 

•01490 

•00617 

-.00000 

46 

-.01324 

•01811 

.01223 

.02952 

.00000 

47 

-.01324 

.01811 

.01223 

-.02952 

.00000 

48 

.036 22 

•00587 

.01490 

-.00617 

-.00000 

49 

•01694 

.01450 

.02036 

•00843 

-.00000 

50 

-.06066 

.00777 

. 00965 

•02329 

•ooooo 

51 

-.06066 

.00777 

.00965 

-.02329 

.00000 

52 

•01694 

.01450 

.02036 

-.00843 

-.00000 

53 

.03094 

.01317 

.02397 

.00993 

•ooooo 

54 

-.05393 

.01141 

.01118 

.02698 

.00000 

55 

-.05393 

•01141 

.01118 

-.02690 

.00000 

56 

.03094 

.01317 

•02397 

-.00993 

•ooooo 

57 

•01268 

.00842 

.020 74 

•00859 

•ooooo 

58 

-.07801 

•00925 

•01240 

.02993 

-ooooo 

59 

-.07801 

• 0092 5 

.01240 

-.02993 

.00000 

60 

.01268 

.00842 

•02074 

-.00859 

.00000 
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OGIVE CYLINDER 8QOY WITH 45 DEGREE SWEEP NACA ASA004 NlD>WlN& 
SINGULARITY PANELING POR SAMPLE CASE 


INTEGRATION OF THE PRESSURE DISTRIBUTION 
ON THE BODY 



MACH* 2. 

0100 



ALPHA- 0. 

0000 


POINT 

X 

Y 

2 

1 

1.00000 

.09575 

*.23116 

2 

1.00000 

.23116 

-.09575 

3 

1.00000 

.23116 

.09575 

4 

1.00000 

.09575 

.23116 

5 

3.22006 

.27308 

-.65928 

6 

3.22006 

•65928 

-.27308 

7 

3.22006 

.65928 

.27308 

B 

3.22006 

.27308 

•65928 

9 

6.06242 

.44633 

-1.07754 

10 

6.08242 

1.07754 

-.44633 

11 

6.06242 

1.07754 

•44633 

12 

6.08242 

•44633 

1.07754 

13 

9.03105 

.55006 

-1.32796 

14 

9.03105 

1.32796 

-.55006 

15 

9.03105 

1.32796 

.55006 

16 

9.03105 

.55006 

1.32796 

17 

11.08 446 

•58637 

-1.41563 

16 

11.08446 

1.4 1563 

-.58637 

19 

11.08446 

1.41563 

.58637 

20 

11.08446 

.56637 

1.4 1563 

21 

13.63090 

.58927 

-1.42262 

22 

13.63090 

1.42262 

-.58927 

23 

13.63090 

1.42262 

•58927 

24 

13.63090 

.58927 

1.42262 

25 

16.48370 

.58927 

-1.42262 

26 

16.48370 

1.42262 

-.58927 

27 

16.48370 

1.42262 

• 58927 

28 

16.48370 

.58927 

1.42262 

29 

18.26145 

•58927 

-1.42262 

30 

18.26145 

1.42262 

-.58927 

31 

18.26145 

1.42262 

.58927 

32 

18.26145 

.58927 

1.42262 

33 

20.03915 

.58927 

-1.42262 

34 

20.03915 

1.42262 

-.58927 

35 

20.03915 

1.42262 

•58927 

36 

20.03915 

.58927 

1.42262 

37 

21.81690 

.58927 

-1.42262 

38 

21.81690 

1.42262 

-.58927 

39 

21.81690 

1.42262 

.58927 

40 

21.81690 

.58927 

1.42262 

4 1 

23.59465 

.58927 

-1.42262 

42 

23.59465 

1.42262 

-.58927 

43 

23.59465 

1.42262 

.58927 

44 

23.59465 

.58927 

1.42262 

45 

25.38175 

•58927 

-1.42262 

46 

25.38175 

1.42262 

-.58927 

47 

25.38175 

1.42262 

•58927 

48 

25.38175 

. 5 892 7 

1.42262 

49 

27.84000 

.58927 

-1.42262 

50 

27.84000 

1.42262 

-.58927 

51 

27.84000 

1.42262 

.58927 

52 

27.84000 

•58927 

1.42262 

53 

31.20000 

•58927 

-1.42262 

54 

31.20000 

1.42262 

-.58927 

55 

31.20000 

1.42262 

.58927 

56 

31.20000 

.58927 

1.42262 

57 

34.75000 

•58927 

-1.42262 

58 

34.75000 

1.42262 

-.58927 

59 

34.75000 

1.42262 

.58927 

60 

34.75000 

.58927 

1.42262 


TOTAL COEFFICIENTS 


ON THE 

BODY 


REFA- 

144.0000 

REFU- 

REFX« 

20.8130 

REFi* 

CN* 

• 0000 


CT* 

.0037 


CM- 

-.0000 


CL» 

• 0000 


CD* 

.0037 


XCP* 

-1.5870 



X/C 

2Y/B 

2/C 

CP 

.02740 

.02875 

-.06942 

•15199 

•02740 

.06942 

-.02875 

.15199 

.02740 

•06942 

•02875 

.15199 

•02740 

•02675 

.06942 

•15199 

.088 22 

.08201 

-.19798 

.11288 

.08822 

•19798 

-.08201 

•11286 

•08822 

.19798 

.08201 

.11288 

.08822 

.08201 

•19798 

.11286 

.16664 

.13403 

-.32359 

.05759 

.16664 

•32359 

-.13403 

.05759 

« 16664 

.32359 

.13403 

.05759 

.16664 

.13403 

•32359 

.05759 

.24743 

•16518 

-.39879 

.00010 

.24743 

•39879 

-.16518 

•00010 

.24743 

.39879 

.16518 

•00010 

•24743 

•16518 

.39879 

.00010 

.30368 

.17609 

-.42511 

-.03988 

.30368 

.42511 

-.17609 

-.03968 

.30368 

.42511 

.17609 

-.03988 

.30368 

.17609 

.42511 

-.03968 

.37345 

.17696 

-.42721 

-.03515 

.37345 

• 42721 

-.17696 

-.03515 

.37345 

.42721 

.17696 

-.03515 

.37345 

.17696 

.42721 

-.03515 

•45161 

•17696 

-.42721 

-.01413 

.45161 

.42721 

-.17696 

-.01413 

.45161 

• 42721 

.17696 

-.01413 

.45161 

•17696 

•42 721 

-.01413 

.50031 

•17696 

-.42721 

-.01486 

.50031 

.42721 

-.17696 

•01602 

.50031 

•42721 

.17696 

.01602 

.50031 

.17696 

•42721 

-.01486 

.54902 

•17696 

-.42721 

-.00118 

•54902 

.42721 

-.17696 

-.00137 

.54902 

•42721 

•17696 

-.00137 

.54902 

•17696 

.42721 

-.00118 

.59772 

.17696 

-.42721 

.00191 

•59772 

.42721 

-.17696 

-.02805 

.59772 

•42721 

.17696 

-.02805 

.59772 

.17696 

•42721 

.00191 

.64643 

•17696 

-.42721 

•GOOlO 

.64643 

.42721 

-.17696 

-.03756 

.64643 

.42721 

•17696 

-.03756 

•64643 

. 17696 

.42721 

.00010 

.69539 

.17696 

-.42721 

-.01188 

.69539 

.42721 

-.17696 

-.03616 

.69539 

.42721 

•17696 

-.03616 

.69539 

•17696 

.42721 

-.01188 

.76274 

•17696 

-.42721 

-.02881 

.76274 

.42721 

-.17696 

-.01596 

. 762 74 

.42721 

.17696 

-.01596 

.762 74 

.17696 

.42721 

-.02881 

.85479 

• 17696 

-.42721 

-.02645 

.854 79 

.42721 

-.17696 

-.02323 

. 854 79 

• 42721 

•17696 

-.02323 

.85479 

.17696 

•42721 

-.02645 

• 952 05 

.17696 

-.42721 

-.01712 

.95205 

•42721 

-.17696 

-.01924 

.95205 

.42721 

.17696 

-.01924 

.95205 

.17696 

.42721 

-.01712 

3.3300 

0.0000 

REFL* 

36. 5000 



CN 

CT 

CM 

POINT 

•03274 

.00687 

.64669 

1 

•01356 

•00887 

•26787 

2 

-.01356 

•00887 

-.26787 

3 

-.03274 

•00887 

-.64669 

4 

.17374 

.03699 

3.03221 

5 

.07197 

.03699 

1.25598 

6 

-.07197 

.03699 

-1.25598 

7 

-.17374 

.03699 

-3.03222 

8 

.15283 

•02101 

2.22865 

9 

•06330 

.02101 

.92314 

10 

-.06330 

.02101 

-.92314 

11 

-.15283 

•02101 

-2.22865 

12 

•00034 

.00002 

.00392 

13 

•00014 

•00002 

•00163 

14 

-.00014 

•00002 

-.00163 

15 

-.00034 

•00002 

-.00392 

16 

-.05458 

-.00077 

-.52986 

17 

-.02261 

-.00077 

-.21948 

18 

.02261 

-.00077 

.21948 

19 

•05458 

-.00077 

.52986 

20 

-.16270 

0.00000 

-1.16854 

21 

-.06739 

0.00000 

-.48402 

22 

.06739 

0.00000 

.48402 

23 

.162 70 

0.00000 

1.16854 

24 

-.02960 

0.00000 

-.12817 

25 

-.01226 

0.00000 

-.05309 

26 

.01226 

0.00000 

.05309 

27 

.02960 

0.00000 

.12817 

28 

-.03112 

0.00000 

-.07942 

29 

•01390 

0.00000 

.03547 

30 

-.01390 

0.00000 

-.03547 

31 

.03112 

0.00000 

.07942 

32 

-.00248 

0.00000 

-.00192 

33 

-.00119 

0.00000 

-.00092 

34 

.00119 

0.00000 

•00092 

35 

•00248 

0.00000 

.00192 

36 

.00400 

0.00000 

-.00401 

37 

-.02434 

0.00000 

.02443 

38 

.02434 

0.00000 

-.02443 

39 

-.00400 

0.00000 

.00401 

40 

.00022 

0.00000 

-.00061 

41 

-.03259 

0.00000 

.09066 

42 

.03259 

0.00000 

-.09066 

43 

-.00022 

0.00000 

.00061 

44 

-.02515 

0.00000 

.11491 

45 

-.03171 

0.00000 

.14490 

46 

•03171 

0.00000 

-.14490 

47 

.02515 

0.00000 

-.11491 

48 

-.10593 

0.00000 

. 74438 

49 

-.02431 

0.00000 

.17085 

50 

.02431 

0.00000 

-.17085 

51 

.10593 

0.00000 

-.74438 

52 

11222 

O.OOQOO 

1.16565 

53 

-.04083 

0.00000 

.42408 

54 

•04083 

0. 00000 

-.42408 

55 

.11222 

0.00000 

-1.16565 

56 

-.07061 

0.00000 

.98404 

57 

-.03288 

0.00000 

.45826 

58 

.03288 

0.00000 

45826 

59 

.07061 

0.00000 

-.98404 

60 
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VELOCITIES ON MING UPPER SURFACEf MACH -2.010 ALPHAS 0.000 


PANEL 

NO . 


1 

2 

J 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 


vortex axial lateral VERTICAL 

STRENGTH VELOCITY VELOCITY VELOCITY 


-.00000 

-.12191 

-.00000 

.02325 

.00000 

•00695 

.00000 

.00518 

.00000 

.00868 

-.00000 

.01821 

-.00000 

.02206 

•00000 

.02414 

•00000 

.02864 

•00000 

.02582 

-.00000 

•02247 

•00000 

-.12687 

-.00000 

-.03264 

-.00000 

-.00618 

-.00000 

.02367 

-.00000 

.01651 

.00000 

.02284 

.00000 

.02992 

-.00000 

.03571 

-.00000 

.03733 

-.00000 

.03439 

-.00000 

•03269 

-.00000 

-.12591 

-.00000 

-.03502 

.00000 

-.01856 

-.00000 

-.00934 

-.00000 

.00056 

-.00000 

.02291 

-.00000 

.04714 

-.00000 

•04349 

-.00000 

.04270 

-.00000 

•04277 

-.00000 

•04111 

-.00000 

-.12906 

.00000 

-.03 764 

-.00000 

-.01850 

-.00000 

-.00791 

- . 00000 

-.00023 

-.00000 

•01061 

•00000 

.01939 

-.00000 

•02434 

-.00000 

.02815 

-.00000 

.02809 

.00000 

.03065 

-.00000 

-.14159 

-.00000 

-.04752 

•00000 

-.01946 

. 00000 

-.01057 

-.00000 

-.00287 

.00000 

.00763 

.00000 

.01648 

-.00000 

•02364 

-.00000 

.02826 

-.00000 

.02754 

-.00000 

.02755 


.14872 

•17879 

-.06309 

.05696 

-.03468 

•03265 

-.02358 

•01709 

-.01760 

•0031 5 

-.02061 

-.01379 

-.01522 

-.02918 

-.00917 

-.03946 

-.01025 

-.04623 

-.00434 

-.04704 

• 00144 

-.04710 

. 15178 

. 17879 

• 04060 

•05696 

•00164 

•03265 

-.05002 

.01709 

-.03580 

.00315 

-* 04163 

-.01379 

-.04226 

-.02918 

-.04406 

-.03946 

-.04194 

-.04623 

-.03595 

-.04704 

03250 

-.04710 

. 14680 

.17879 

•04206 

•05696 

•02123 

.03265 

•00966 

.01709 

-.00363 

•00315 

-.03670 

-.01379 

-.07500 

-.02918 

06282 

-.03946 

-.05353 

-.04623 

-.05200 

-.04704 

-• 04856 

-.04710 

.15068 

.17879 

.04531 

.05696 

.01967 

•03265 

.00540 

•01709 

-.00397 

•00315 

-.01688 

-.01379 

-.02695 

-.02918 

-.03331 

-.03946 

-.03980 

-.04623 

-.04334 

-.04704 

-.04243 

-.04710 

. 17134 

•17879 

.06159 

.05696 

.02050 

.03265 

.00896 

.01709 

-.00036 

.00315 

-.01249 

-.01379 

-.02250 

-.02918 

-.03263 

-.03946 

-.04052 

-.04623 

-.04298 

-.04704 

-.04707 

-.04710 
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UGIVE CYLINDER dOUY WITH DEGREE SWEEP NACA 65A00A MIO'WING 
SINGULARITY PANELING FOR SAMPLE CASE 


INTEGRATiON UF THE PRESSURE OISTRUUTION 
ON THE WING UPPER SURFACE 


HACH» ^.0100 
ALPHA* 0.0000 


POINT 

X 

V 

1 

x/c 

2Y/8 

2/C 

CP 

CN 

CT 

CM 

POINT 

I 

16. 76499 

2.30734 

0.00000 

•05000 

•19228 

0.00000 

.07631 

-.08627 

.01017 

-.34921 

1 

2 

17.61117 

2.30734 

0.00000 

.15000 

• 19228 

0.00000 

-.03365 

.03707 

-.00166 

.11869 

2 

3 

18.45734 

2.30734 

0.00000 

•25000 

•19228 

0.00000 

-.01353 

.01490 

-.00037 

•03511 

i 


19.30352 

2.30734 

0.00000 

.35000 

.19228 

0.00000 

-.01427 

.01572 

-.00016 

.02373 

4 

5 

20.14970 

2.30734 

0.00000 

.45000 

.19228 

0.00000 

-.02670 

.02941 

.00016 

.01951 

5 

6 

20.99587 

2.30734 

0.00000 

.55000 

.19228 

0.00000 

-.03980 

.04384 

.00094 

-.00602 

6 

7 

21.84205 

2.30734 

0.00000 

.65000 

.19228 

Q. CO 000 

-.04580 

.05046 

.00173 

-.05192 

7 

8 

22.68823 

2.30734 

0.00000 

.75000 

•19228 

0.00000 

-.05237 

.05769 

.00247 

10819 

8 

9 

23.53441 

2.30734 

0.00000 

.65000 

.19228 

0.00000 

-.05418 

.05969 

•00278 

-.16243 

9 

10 

24.38058 

2. 30734 

0.00000 

.95000 

.19228 

0.00000 

-.04851 

.05343 

.00252 

19063 

10 

11 

1 a . 82 5 7 7 

4.12568 

0.00000 

.05000 

•34381 

0.00000 

. 14091 

-.24418 

•02878 

-.48523 

11 

12 

19.55072 

4.12568 

0.00000 

•15000 

•34381 

0.00000 

.03717 

-.06441 

•00289 

-.08131 

12 

13 

20.27568 

4.12568 

0.00000 

.25000 

.34381 

0.00000 

-.01839 

.03187 

-.00079 

•01712 

13 

14 

21.00063 

4. 12568 

0.00000 

.35000 

.34381 

0.00000 

-.04077 

.07065 

-.00071 

-.01326 

14 

IS 

21 . 72559 

4. 12568 

0.00000 

•45000 

• 3438 1 

0.00000 

-.03961 

.06864 

.00037 

-.06264 

15 

16 

22.45054 

4.12568 

0.00000 

.55000 

•34381 

0.00000 

-.05263 

.09119 

.00196 

-.14933 

16 

17 

23.17550 

4.12568 

0.00000 

.65000 

•34381 

0.00000 

-.06497 

.11258 

.003 86 

-.26596 

17 

18 

23.90045 

4. 12568 

0.00000 

. 75000 

•34381 

0.00000 

-.07210 

. 12494 

.00535 

-.38575 

13 

19 

24.62541 

4. 12568 

0.00000 

.85000 

•34381 

0.00000 

-.07093 

.12292 

.00573 

-.46861 

19 

20 

25.35036 

4.12568 

0.00000 

.95000 

.34381 

0.00000 

-.06657 

.11534 

•00543 

-.52336 

20 

21 

21.51108 

6.4950/ 

0.00000 

.05000 

•54126 

0.00000 

.14331 

-.19052 

.02246 

. 13300 

21 

22 

22.07808 

6.49507 

0.00000 

.15000 

.54126 

0.00000 

.05227 

-.06950 

.00311 

.08792 

22 

23 

22.64507 

6.49507 

0.00000 

.25000 

.54126 

0.00000 

.02754 

-.03662 

.00091 

.06708 

23 

24 

23.21207 

6.49507 

0.00000 

•35000 

.54126 

0.00000 

.00870 

-.01157 

.00012 

.02 7 76 

24 

2S 

23.77906 

6.49507 

0.00000 

.45000 

•54126 

0.00000 

-.02337 

.03107 

•0001 7 

-.09215 

25 

26 

24.34606 

6.49507 

0.00000 

.55000 

.54126 

0.00000 

-.06915 

-09194 

.00198 

-.32482 

26 

27 

24.91305 

6.49507 

0.00000 

.65000 

•54126 

0.00000 

-.08921 

•11860 

.00407 

-.48628 

27 

28 

25.48005 

6-49507 

0.00000 

.75000 

•54126 

0.00000 

-.08484 

. 11279 

•00483 

-.52638 

28 

29 

26.04704 

6.49507 

0.00000 

.85000 

.54126 

0.00000 

-.08398 

.11164 

•00521 

-.58435 

29 

30 

26.61404 

6.49507 

0.00000 

.95000 

.54126 

O.OQOQO 

-.08242 

.10958 

•00516 

-.63568 

30 

31 

24.16649 

8.83808 

0.00000 

.05000 

.73651 

0.00000 

.14910 

-.14335 

.01690 

.48071 

31 

32 

24.57728 

8.63808 

0.00000 

.15000 

.73651 

0.00000 

.05497 

-.05285 

.00237 

.19693 

32 

33 

24.98808 

8.83808 

0.00000 

.25000 

.73651 

0.00000 

.02608 

-.02507 

.00062 

.10468 

33 

34 

25.39887 

8.63808 

0.00000 

.35000 

• 73651 

0.00000 

•00806 

-.00775 

.00008 

•03552 

34 

35 

25.80967 

8.83808 

0.00000 

.45000 

.73651 

0.00000 

-.01044 

.01004 

.00005 

-.05018 

35 

36 

26.22046 

8.83808 

0.00000 

.55000 

•73651 

0.00000 

-.03020 

.02904 

•00062 

-.15701 

36 

37 

26.63126 

8.83808 

0.00000 

•65000 

.73651 

0.00000 

-.04417 

.04247 

.00146 

-.24711 

37 

38 

27.04205 

6.63808 

0.00000 

.75000 

•73651 

Q. 00000 

-.05322 

.05117 

.00219 

-.31073 

38 

39 

27.45285 

8.83808 

0.00000 

.85000 

.13651 

0.00000 

-.05728 

.05507 

•00257 

-.36568 

39 

40 

27.86364 

8.83808 

0.00000 

.95000 

.73651 

0.00000 

-.05967 

-05737 

.00270 

-.40449 

40 

41 

26.58702 

10.97384 

0.00000 

.05000 

.91449 

0.00000 

.16990 

-.08500 

.01002 

.49082 

41 

42 

26.85543 

10.97364 

0.00000 

•15000 

.91449 

0.00000 

.06594 

-.03299 

.00148 

. 19936 

42 

43 

27.12384 

10.97384 

0.00000 

.25000 

.91449 

0.00000 

.02978 

-.01490 

.00037 

.09402 

43 

44 

27.39225 

10.97384 

0.00000 

.35000 

♦91449 

0.00000 

.01342 

-.00672 

.00007 

•04419 

44 

45 

27.66066 

10.97384 

0.00000 

.45000 

•91449 

0.00000 

-.00483 

•00242 

•00001 

-.01656 

45 

46 

27.92907 

10.97384 

0.00000 

.55000 

•91449 

0.00000 

-.02441 

.01221 

.00026 

-.08692 

46 

47 

28.19748 

10.97384 

0.00000 

.65000 

•91449 

0.00000 

-.04067 

.02035 

.00070 

15026 

47 

48 

28.46589 

10.97364 

0.00000 

•75000 

•91449 

0.00000 

-.05268 

.02636 

.00113 

-.20172 

48 

49 

28.73430 

10.97384 

0.00000 

.85000 

•91449 

0.00000 

-.05689 

.02847 

.00133 

-.22549 

49 

50 

29.00271 

10.97384 

0.00000 

.95000 

.91449 

0.00000 

-.05653 

.02828 

•00133 

-.23164 

50 
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VELOCITIES ON HINC LOWER SURFACE* HACH-2.010 ALPHA* 0.000 


PANEL 

VORTEX 

AXIAL 

LATERAL 

VERTICAL 

NO. 

STRENGTH 

VELOCITY 

VELOCITY 

VELOCITY 


1 

-.00000 

-.12191 

. 14872 

-.17879 

2 

-.00000 

•02325 

-. 06309 

-.05696 

3 

.00000 

• C 0695 

-.03468 

-.03265 

A 

.00000 

.00518 

-. 02358 

-.01709 

5 

.00000 

•00868 

-.01760 

-.00315 

6 

-.00000 

.01821 

-.02061 

•01379 

7 

-.00000 

.02206 

-.01522 

•02918 

8 

•00000 

•02414 

-.00917 

•03946 

9 

. 00000 

•02864 

-.01025 

•04623 

10 

•00000 

.02582 

-.00434 

.04704 

11 

-.00000 

.02247 

.00144 

.04710 

12 

.00000 

-.12687 

•15178 

-.17879 

13 

-.00000 

-.03264 

•04060 

-.05696 

14 

-.00000 

-.00616 

.00164 

-.03265 

15 

-.00000 

.02367 

-.05002 

-.01709 

16 

-.00000 

.01651 

-.03580 

-.00315 

17 

•00000 

.02284 

-.04163 

•01379 

18 

.00000 

.02992 

-.04226 

•02918 

19 

-.00000 

.03571 

-.04406 

•03946 

20 

-.00000 

.03733 

-.04194 

.04623 

21 

-.00000 

.03439 

-.03595 

.04704 

22 

-.00000 

•03269 

-.03250 

•04710 

23 

-•coooo 

-.12591 

• 14680 

-.17879 

24 

-.00000 

-.03502 

•04206 

-.05696 

25 

•00000 

-.01856 

•02123 

-.03265 

26 

-.00000 

-.00934 

•00966 

-.01709 

27 

-.00000 

•00056 

-.00363 

-.00315 

28 

-.00000 

.02291 

-.03670 

.01379 

29 

-.00000 

.04714 

-.07500 

.02918 

30 

-.00000 

•04349 

-.06282 

•03946 

31 

-.00000 

•04270 

-.05353 

•04623 

32 

-.00000 

.04277 

-. 05200 

.04704 

33 

-.00000 

.04111 

-.04856 

•04710 

34 

-.00000 

-.12906 

.15068 

-.17879 

35 

.00000 

-.03764 

.04531 

-.05696 

36 

-.00000 

-.01850 

.01967 

-.03265 

37 . 

-.00000 

-.00791 

.00540 

-.01709 

38 

-.00000 

-.00023 

-.00397 

-.00315 

39 

-.00000 

.01061 

-.01688 

•01379 

40 

•00000 

•01939 

-.02695 

•02918 

41 

-.00000 

•02434 

-.03331 

•03946 

42 

-.00000 

.02815 

-.03980 

•04623 

43 

-.00000 

•02809 

-.04334 

.04704 

44 

.00000 

.03065 

-.04243 

•04710 

45 

-.00000 

-.14159 

. 17134 

-.17879 

46 

-.00000 

-.04752 

.06159 

-.05696 

47 

.00000 

-.01946 

.02050 

-.03265 

48 

.00000 

-.01057 

.00896 

-.01709 

49 

-.00000 

-.00287 

-.00036 

-.00315 

50 

•ooooo 

.00763 

-.01249 

•01379 

51 

.00000 

.01648 

-. 02250 

.02918 

52 

-.00000 

.02364 

-.03263 

.03946 

53 

-.00000 

.02826 

-.04052 

.04623 

54 

-.00000 

.02754 

-.04298 

•04 T 04 

55 

-.00000 

.02755 

-.04707 

.04710 
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aoiVE CYLINDER BCDV I^ITH DEGREE SmEEP NACA 6SA004 NIO-WING 
SINGULARITY PANELING FOR SAHPLE CASE 

INTEGRATION OF THE PRESSURE DISTRIBUTION 
ON THE MING LOmER SURFACE 


MACH>i 2.0100 
ALPHA- O.COOO 


POINT 

X 

Y 

1 

x/c 

2Y/B 

1 

16. 76499 

2.30734 

0.00000 

.05000 

.19226 

2 

17.61117 

2.30734 

0.00000 

.15000 

•19228 

3 

10.45734 

2.30734 

0.00000 

.25000 

.19220 

4 

19.30352 

2.30734 

0.00000 

.35000 

• 1922 0 

5 

20.14970 

2. 30734 

0.00000 

.45000 

•19220 

6 

20.99507 

2.30734 

0.00000 

.55000 

. 19226 

7 

21.04205 

2.30734 

0.00000 

.65000 

•19220 

b 

22 .60023 

2.30734 

0.00000 

.75000 

.19220 

N 

2J.53441 

2.30734 

0.00000 

.85000 

•19220 

10 

24.30050 

2.30734 

0.00000 

*95000 

.19226 

11 

10.02577 

4.12568 

0.00000 

.05000 

.34361 

12 

19.55072 

4. 12560 

0.00000 

.15000 

.34381 

13 

20.27560 

4.12566 

0.00000 

.25000 

.34301 

14 

21.00063 

4. 12566 

0.00000 

.35000 

•34301 

15 

21.72559 

4.12568 

0.00000 

.45000 

.34301 

16 

22.45054 

4.12560 

0.00000 

.55000 

.34301 

17 

23.17550 

4.12560 

0.00000 

.65000 

•34301 

IB 

23. 90045 

4.12566 

0.00000 

.75000 

.34381 

19 

24.62541 

4.12560 

0,00000 

.05000 

.34301 

20 

25.35036 

4.12560 

0.00000 

.95000 

.34301 

21 

21.51100 

6.49507 

0.00000 

.05000 

•54126 

22 

22.O70C0 

6.4950/ 

0.00000 

.15000 

.54126 

23 

22.64507 

6.49507 

0.00000 

.25000 

.54126 

24 

23.21207 

6.49507 

0.00000 

.35000 

.54126 

25 

23.77906 

6.49507 

0.00000 

.45000 

.54126 

26 

24.34606 

6.49507 

0.00000 

.55000 

.54126 

27 

24.91305 

6.49507 

0.00000 

.65000 

• 54126 

20 

25.40005 

6.49507 

0.00000 

.75000 

•54126 

29 

26.04704 

6.49507 

0.00000 

•05000 

•54126 

30 

26.61404 

6.49507 

0.00000 

.95000 

•54126 

31 

24. 16649 

0.03000 

0.00000 

.05000 

.73651 

32 

24.57720 

6.03000 

0.00000 

.15000 

. 7365 1 

33 

24.90000 

0.03808 

0.00000 

.25000 

•73651 

34 

25.39007 

0.03600 

0.00000 

•35000 

•73651 

35 

25.00967 

0.83800 

0.00000 

.45000 

•73651 

36 

26.22046 

0.83008 

0.00000 

•55000 

.73651 

37 

26.63126 

8.83800 

0.00000 

•65000 

•73651 

3S 

27.04205 

8.83008 

0.00000 

.75000 

•73651 

39 

27.45205 

8.03608 

0.00000 

.05000 

• 73651 

40 

27.06364 

0.03600 

0.00000 

•95000 

•73651 

41 

26.50702 

10.97384 

0.00000 

•05000 

•91449 

42 

26.05543 

10.97364 

0.00000 

.15000 

.91449 

43 

27.12304 

10.97384 

0.00000 

•25000 

•91449 

44 

27.39225 

10.97304 

0.00000 

.35000 

•91449 

45 

27.66066 

10.973 84 

0.00000 

.45000 

• 91449 

46 

27.92907 

10.97384 

0.00000 

.55000 

•91449 

47 

20.19740 

10.97384 

0.00000 

•65000 

•91449 

40 

20.46589 

10.97304 

0.00000 

.75000 

•91449 

49 

20.73430 

10.97304 

0.00000 

.85000 

.91449 

50 

29.00271 

10.97304 

0.00000 

•95000 

.91449 


2/C 

CP 

CN 

CT 

CH 

POINT 

0.00000 

.07031 

•08627 

• 01017 

.34921 

1 

0.00000 

-.03365 

-.03707 

-.00166 

-.11869 

2 

0.00000 

-.01353 

-.01490 

-.00037 

-.03511 

3 

0.00000 

-.01427 

-.01572 

-.00016 

-.02373 

4 

0.00000 

-.02670 

-.02941 

•00016 

-.01951 

5 

0. 00000 

-.03900 

-.04304 

•00094 

.00802 

6 

0.00000 

-.04500 

-.05046 

.00173 

.05192 

7 

0.00000 

-.05237 

-.05769 

•00247 

.10819 

8 

0.00000 

-.05410 

-.05969 

.00278 

.16243 

9 

0.00000 

-.04051 

-.05343 

.00252 

.19063 

10 

0.00000 

.14091 

.24410 

•02878 

•46523 

11 

0.00000 

.03717 

•06441 

•00209 

.08131 

12 

0.00000 

-.01639 

-.03107 

-.00079 

-.01712 

13 

0.00000 

-.04077 

-.07065 

-.00071 

.01326 

14 

0.00000 

-.03961 

-.06864 

.00037 

.06264 

15 

0.00000 

-.05263 

-.09119 

•00196 

.14933 

16 

Q.OOOQO 

-.06497 

-.11258 

.00386 

.26596 

17 

0.00000 

-.07210 

-.12494 

•00535 

.38575 

10 

0.00000 

-.07093 

12292 

.00573 

.46661 

19 

0.00000 

-.06657 

-.11534 

.00543 

.52336 

20 

0.00000 

.14331 

•19052 

•02246 

13300 

21 

0.00000 

.05227 

.06950 

•00311 

-.08792 

22 

0.00000 

•02754 

•03662 

.00091 

-.06708 

23 

0.00000 

•00070 

•01157 

.00012 

-.02776 

24 

0.00000 

-.02337 

-.03107 

.00017 

.09215 

25 

0.00000 

-.06915 

-.09194 

•00196 

.32482 

26 

0.00000 

-.08921 

-.11860 

.00407 

•48620 

27 

0.00000 

-.06484 

-.112 79 

.00483 

•52638 

20 

0.00000 

-.06396 

-.11164 

.00521 

.50435 

29 

0,00000 

-.08242 

-.10958 

.00516 

•63568 

30 

0.00000 

.14910 

.14335 

•01690 

-.40071 

31 

0.00000 

.05497 

.05285 

•00237 

-.19893 

32 

0.00000 

.02606 

•02507 

•00062 

10468 

33 

0.00000 

•00806 

.00775 

•00008 

-.03552 

34 

0.00000 

-.01044 

-.01004 

.00005 

.05018 

35 

0.00000 

-.08020 

-.02904 

.00062 

.15701 

36 

0.00000 

-.04417 

-.04247 

•00146 

.24711 

37 

0.00000 

-.05322 

-.05117 

•00219 

.31873 

30 

0.00000 

-.05728 

-.05507 

.00257 

.36566 

39 

0.00000 

-.05967 

-.05737 

.00270 

.40449 

40 

0.00000 

.16990 

.06500 

•0 1002 

-.49082 

41 

0.00000 

•06594 

.03299 

•00148 

-. 19936 

42 

0.00000 

•02978 

•01490 

.00037 

-.09402 

43 

0.00000 

.01342 

.00672 

.00007 

-.04419 

44 

0.00000 

-.00483 

-.00242 

.00001 

•01656 

45 

0.00000 

-.02441 

-.01221 

•00026 

.06692 

46 

0.00000 

-.04067 

-.02035 

.00070 

.15026 

47 

0.00000 

-.05268 

-.02636 

.00113 

.20172 

40 

0.00000 

-.05689 

-.02847 

.00133 

.22549 

49 

0.00000 

-.05653 

-.02828 

.00133 

.23164 

50 
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TOTAL 

COEFFICIENTS 


ON THE 

MING 


REFA- 

LA4.0000 

REFB> 

REFX- 

20.8130 

REFZ' 

CN» 

CT- 

CM« 

CL» 

CO- 

XCP- 

-.0000 

• 0046 

• 0000 
-.0000 

.0046 
-, 1056 



U.OOOO REFC* 

0.0000 


TOTAL 

COEFFICIENTS 



ON THE 

COMPLETE CONFIGURATION 


REFA- 

144.0000 

REF8- 

12.0000 

RE FX- 

20.8130 

PEFZ- 

0.0000 

CN- 

CT- 

CM- 

CL- 

Cl>- 

XCP- 

.0000 

.0083 

-.0000 

.0000 

.0083 

-3.5710 




6.8900 


6.8900 
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SECTIUN COEfflCICNTS 
ON THE MING 


UELV- 

1.3030 

8EFL 

CN- 

-.0000 


CT- 

.0034 


CM« 

-.0000 


CL- 

-.0000 


co» 

.0034 


XCP- 

• 3336 



UELY« 

2.4000 

REFL 

CN- 

-.0000 


CT- 

.0061 


CM- 

• 0000 


CL- 

-.0000 


co« 

• 0061 


XCP- 

-.1531 



OELY- 

2.3600 

«£FL 

CN- 

-.0000 


CT- 

.0072 


CM- 

• 0000 


CL- 

-.0000 


co- 

.00 72 


XCP- 

-.4477 



6«8900 XLE* 


6.8900 XLE« 


6.8900 XLE* 


SECTION COEFI^ICIENTS 
ON THE MING 


□ ELY- 

2.3700 

CN- 

-.0000 

CT- 

.0061 

CM- 

.0000 

CL- 

-.0000 

CD- 

• 0061 

XCP- 

-.6329 


OELY- 

1.9000 

CN- 

.0000 

CT- 

.0065 

CM« 

• 0000 

CL- 

.0000 

CO- 

.0065 

XCF- 

3.0337 


CFSTAG • 2.656S0 CFCAIT - 

time - 253.82300 
time « 256.24700 
TIME • 274.05300 


MEft- 




1.13092 


6.8900 XLE- 


6.8900 XLE» 


CF9AC ■ -.33360 


16.3419 


16.4633 


21.2276 


23.9611 


26.4528 
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VtLaCITieS ON BOUVi MACH»2*010 alpha* 


5.000 


PANEL 

SOURCE 

AXIAL 

NO. 

STRENGTH 

VELOCITY 

i 

.29821 

-.12592 

2 

.23814 

-.10863 

3 

.15320 

-.08418 

4 

.a9313 

-.06689 

5 

•29647 

-.09585 

6 

.22353 

-.08065 

7 

. L2036 

-.05914 


•04744 

-.04393 

9 

.25703 

-.05884 

10 

. 17238 

-.04502 

11 

.05267 

-.02549 

12 

-.03198 

-.01167 

13 

.228 72 

-.01979 

14 

.12507 

-.00925 

15 

-.02150 

.00564 

16 

-.12515 

.01618 

17 

. 16034 

.00840 

18 

.05307 

.01544 

19 

-.09862 

•02540 

£0 

-.20589 

.03245 

21 

.19295 

.01316 

£2 

.07831 

.01601 

23 

-.08383 

.02000 

24 

-.19848 

•02282 

25 

.23512 

•00731 

26 

•11426 

• 00 720 

27 

-.05666 

.00703 

28 

17752 

.00692 

29 

-.15629 

•Q0058 

30 

-.07665 

-.03603 

31 

•12388 

.01975 

32 

. 15214 

.01440 

33 

.00748 

-.03345 

34 

.00491 

-.03146 

35 

.01679 

.03255 

36 

-.053 76 

•03414 

37 

•05941 

-.04515 

38 

-.00215 

-.00961 

39 

-.03768 

.03816 

40 

-.13940 

.04306 

41 

10269 

-.01962 

42 

-.09683 

-.01116 

43 

•08433 

.04942 

44 

.11970 

• 01950 

45 

.09772 

-.03097 

46 

- .03546 

-.03552 

47 

•00907 

•07167 

48 

-.02556 

.04271 

49 

.04061 

-.02359 

50 

-.06008 

-.00391 

51 

-.06077 

.01942 

52 

-.00687 

.05255 

53 

.16144 

-.00624 

54 

-.03317 

-.00266 

55 

-.07429 

.02541 

56 

-.09979 

•03250 

57 

. 16721 

.00256 

58 

-.01301 

•01136 

59 

-.14241 

.00707 

60 

-.14195 

.01420 


LATERAL 

VERTICAL 

NORMAL 

VELOCITY 

VELOCITY 

VELOCITY 

. 12742 

-.27006 

.32977 

.24926 

-.06568 

.28260 

. 166 72 

.10663 

•21590 

.04488 

•14592 

• 168 7 3 

• 11639 

-.23061 

.27645 

.21498 

-.03868 

.22928 

.12165 

.10076 

•16258 

.02305 

.10602 

.11541 

.09712 

-.17579 

.20705 

.16 316 

-.00891 

.15988 

.06234 

.08451 

.09316 

-.00371 

.04973 

•04601 

.07956 

-.11702 

•13973 

. 11155 

.02886 

.09256 

-.00234 

•07409 

.02586 

-.03433 

-.00781 

-.02131 

•06631 

-.07381 

.09346 

.07345 

•05586 

.04629 

-.04907 

•06596 

-.02042 

-.05621 

-.04942 

-.06759 

.06479 

-.06032 

.08052 

.06479 

•06925 

.03335 

-.06479 

•06925 

-.03335 

-.06479 

-.06032 

-.06052 

.06546 

-.05996 

•08052 

• 06566 

.07137 

.03335 

-.06566 

.07137 

-.03335 

-.06566 

-.05996 

-.06052 

-.00204 

-.08800 

.06052 

.00504 

-.07486 

.03330 

-.02128 

-.03592 

-.03340 

.00204 

-.08800 

-.08052 

-.02645 

-.09813 

.06054 

-.02400 

-.14499 

.03331 

•01123 

-.11439 

-.03340 

-.01240 

-.08200 

-.08050 

-.01938 

-.09522 

.08056 

-.00666 

-.10331 

.03338 

.00723 

-.10454 

-.03333 

-.00056 

-.06689 

-.08049 

•03603 

-.07224 

.08053 

.01946 

-.04030 

.03340 

-.00515 

-.07459 

-.03330 

-.03113 

-.07425 

-.08051 

.01107 

-.08255 

.08050 

-. 01010 

-.11170 

.03342 

•03465 

-.17065 

-.03329 

•01666 

-.09491 

-.08034 

.04121 

-.07005 

.08049 

•03051 

-.01369 

.03343 

-.01112 

-.06012 

-.03328 

-.00061 

-.08694 

-.08055 

•06840 

-.05875 

.08053 

.04641 

.02451 

.03349 

02383 

-.02926 

-.03321 

-.02073 

-.07855 

-.06051 

•06929 

-.05847 

.08053 

.06502 

.06939 

.03351 

-.03996 

•00974 

-.03319 

-.02790 

-.07559 

-.06051 
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OGIVE CYLINDER 800V KITH 45 DEGREE SWEEP NACA 65A004 HIO-WlNG 
SINGULARITY PANELING FOR SAMPLE CASE 


INTEGRATION OF THE PRESSURE DISTRIBUTION 
ON THE BODY 


MACH- 2.0100 
ALPi^A» 5.0000 


POINT 

X 

Y 

2 

X/C 

2V/8 

1 

1.00000 

.09575 

-.23116 

.02740 

.02875 

2 

1.00000 

.23116 

-.09575 

•02740 

.06942 

3 

1.00000 

.23116 

•09575 

.02740 

•06942 

4 

1.00000 

.09575 

.23116 

.02740 

•02875 

5 

3.22006 

.27308 

-.65928 

.08822 

.08201 

6 

3.22006 

.65928 

-.27306 

•08822 

.19798 

7 

3.22006 

•65928 

.2 7308 

•08822 

•19798 

8 

3.22006 

•27300 

•65928 

.08822 

.06201 

9 

6.08242 

•44633 

-1.07754 

. 16664 

.13403 

10 

6.06242 

1.07754 

-.44633 

.16664 

•32359 

11 

6.08242 

1.07754 

.44633 

.16664 

.32359 

12 

6.08242 

.44633 

1.07754 

•16664 

•13403 

13 

9.03105 

.55006 

-1.32796 

.24743 

.16518 

14 

9.03105 

1.32796 

-.55006 

.24743 

.39879 

15 

9.03105 

1.32796 

•55006 

« 24743 

.39879 

16 

9.03105 

.55006 

1.32796 

•24743 

.16518 

17 

11.08446 

.36637 

-1.41563 

•30368 

•17609 

18 

11.08446 

1.41563 

-.58637 

.30368 

.42511 

19 

11.08446 

1.41563 

.56637 

.30368 

.42511 

20 

11.08446 

.58637 

1.41563 

.30368 

• 17609 

21 

13.63090 

.58927 

-1.42262 

.37345 

.17696 

22 

13.63090 

1.42262 

-.58927 

.37345 

• 42721 

23 

13.63090 

1.42262 

.58927 

.37345 

.42721 

24 

13.63090 

.58927 

1.42262 

.37345 

•17696 

25 

16.48370 

.56927 

-1.42262 

.45161 

.17696 

26 

16.48370 

1.42262 

-.58927 

•45161 

.42 721 

27 

16.48370 

1.42262 

.58927 

•45161 

.42721 

28 

16.48370 

.58927 

1.42262 

.45161 

.17696 

29 

18.26145 

.58927 

-1.42262 

.50031 

.17696 

30 

18.26145 

1.42262 

-.58927 

. 50031 

.42721 

31 

10.26145 

1.42262 

.56927 

.50031 

.42721 

32 

18.26145 

.58927 

1.42262 

. 50031 

•17696 

33 

20.03915 

•58927 

-1.42262 

.54902 

.17696 

34 

20.03915 

1.42262 

-.58927 

•54902 

.42721 

35 

20.03915 

1.42262 

.58927 

.54902 

• 42721 

36 

20.03915 

.58927 

1.42262 

.54902 

. 17696 

37 

21.81690 

.58927 

-1.42262 

.59772 

• 17696 

38 

21.81690 

1.42262 

-.58927 

.59772 

• 42721 

39 

21.81690 

1.42262 

.58927 

.59772 

•42721 

40 

21.81690 

.58927 

1.42262 

.59772 

.17696 

41 

23.59465 

.58927 

-1.42262 

.64643 

.17696 

42 

23.59465 

1.42262 

-.58927 

.64643 

.42721 

43 

23.59465 

1.42262 

.58927 

.64643 

.42721 

44 

23. 59465 

.58927 

1.42262 

•64643 

.17696 

45 

25.38175 

.58927 

-1.42262 

.69539 

. W696 

46 

25.38175 

1.4^262 

-.58927 

.69539 

.42721 

47 

25.30175 

1.42262 

•58927 

•69539 

•42721 

48 

25.38175 

.58927 

1.42262 

.69539 

.17696 

49 

27.84000 

.58927 

-1.42262 

.76274 

.17696 

50 

27.84000 

1.42262 

-.58927 

.76274 

.42721 

51 

27.84000 

1.422 62 

.58927 

• 762 74 

.42721 

52 

27.84000 

.58927 

1.42262 

. 762 7 9 

.17696 

53 

31.20000 

.58927 

-1.42262 

.85479 

.17696 

54 

31.20000 

1.42262 

-.58927 

.854 79 

.42721 

55 

31.20000 

1.42262 

.58927 

.85479 

.42721 

56 

31.20000 

.58927 

1.42262 

.854 79 

. 17696 

57 

34.75000 

.58927 

-1.42262 

• 95205 

•17696 

58 

34. 75000 

1.42262 

-.58927 

•95205 

.42 721 

59 

34.75000 

1.42262 

.58927 

.95205 

.42721 

60 

34.75000 

.58927 

1.42262 

-95205 

.17696 


Z/C 

CP 

CN 

CT 

CM 

POINT 

-.06942 

.23352 

.05031 

•01362 

.99362 

1 

-.02875 

.17365 

.01550 

.01013 

.30604 

2 

•02875 

.11402 

-.01017 

•00665 

-.20096 

3 

.06942 

.08674 

-.01869 

•00506 

-.36905 

4 

-.19798 

.16116 

•27883 

•05936 

4.86624 

5 

-.06201 

.12674 

.08080 

•04153 

1.41019 

6 

.08201 

•07717 

-.04920 

•02529 

-.85865 

7 

.19798 

.05850 

-.09004 

.01917 

-1.57145 

8 

-.32359 

.11548 

•30647 

.04213 

4.46908 

9 

-.13403 

•06657 

.07318 

•02429 

1.06719 

10 

.13403 

.02498 

-.02746 

.00911 

-.40042 

11 

.32359 

.01211 

-.03215 

.00442 

-.46678 

12 

-.39879 

•04100 

.13514 

•00869 

1.58062 

13 

-.16518 

.00005 

.00006 

.00001 

.00073 

14 

•16518 

-.02880 

.03933 

-.00611 

.45998 

15 

-39879 

-.03132 

.10324 

-.00664 

1.20759 

16 

-.42511 

-.01359 

-.01860 

-.00026 

-.18057 

17 

-.17609 

-.04665 

-.02656 

-.00090 

-.25787 

18 

•17609 

-.06477 

.03672 

-.00125 

.35647 

19 

.42511 

-.05862 

.08050 

-.00113 

.78153 

20 

-.42721 

-.02320 

-.10738 

0.00000 

-.77118 

21 

-.17656 

-.05041 

-.09667 

0.00000 

-.69426 

22 

•17696 

-.05761 

. 11046 

0. 00000 

•79336 

23 

.42721 

-.04146 

.19190 

0.00000 

1.37824 

24 

-.42721 

-.01193 

-.02500 

0. 00000 

10824 

25 

-.17696 

-.03493 

-.03031 

0.00000 

-.13124 

26 

.17696 

-.03462 

.03005 

0.00000 

.13009 

27 

•42721 

-.01116 

.02338 

0.00000 

.10121 

23 

-.42721 

•00649 

.01359 

0.00000 

.03468 

29 

-.17696 

.08422 

.07309 

0.00000 

.18649 

30 

.17696 

-.03398 

•02949 

0.00000 

.07524 

31 

.42721 

-.02084 

•04367 

0.00000 

.11142 

32 

-.42721 

•07774 

. L6288 

0.00000 

.12604 

33 

-.17696 

.06984 

.06061 

0.00000 

.04690 

34 

.17696 

-.05572 

.04836 

0.00000 

.03742 

35 

.42721 

-.05800 

• 12152 

0.00000 

.09404 

36 

-.42 721 

.10455 

.21906 

0.00000 

-.21991 

37 

-.17696 

.02705 

.02346 

0. 00000 

-.02357 

38 

.17696 

-.06541 

•05676 

0.00000 

-.05698 

39 

.42721 

-.07379 

.15461 

0.00000 

15522 

40 

-.42721 

.04685 

.09816 

0.00000 

-.27305 

41 

-.17696 

.02789 

.02420 

0.00000 

-.06733 

42 

.17696 

-.08500 

.07377 

0.00000 

-.20520 

43 

.42721 

-.03171 

•06643 

0.00000 

-.18478 

44 

-.42721 

.07302 

.15461 

0.00000 

70637 

45 

-.17696 

.08247 

.07233 

0.00000 

-.33046 

46 

•17696 

-.12755 

.11186 

0.00000 

-.51105 

47 

.42721 

-.07352 

.15565 

0.00000 

-.71114 

48 

-.42721 

.05484 

.20165 

0.00000 

-1.41697 

49 

-.17696 

.00912 

.01390 

0.00000 

-.09766 

50 

.17696 

-.03126 

.04764 

0.00000 

-.33476 

51 

.42721 

-.09019 

.33164 

0.00000 

-2.33045 

52 

-.42721 

•01468 

•06230 

0.00000 

-.64709 

53 

-.17696 

-.00172 

-.00303 

0.00000 

.03142 

54 

•17696 

-.04535 

.07970 

0.00000 

-.82789 

55 

.42 721 

-.05531 

.23467 

0.00000 

-2.43756 

56 

-.42721 

-.00317 

-.01308 

0.00000 

.18229 

57 

-.17696 

-.04199 

-.07174 

0.00000 

.99979 

58 

.17696 

-.01723 

.02943 

0.00000 

-.41022 

59 

.42721 

-.02134 

.08802 

0.00000 

- 1.22679 

60 


TOTAL COEFFICIENTS 
ON THE BODY 


REFA* 

144.0000 

REFD* 

3.3300 

REFX* 

20. 8130 

REFZ- 

0.0000 

CN* 

.0526 



CT- 

.0035 



CM* 

.0053 



CL* 

.0521 



CU- 

.0081 



XCP* 

.1021 
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VELOCLTIES ON NlNG UPPER SURFACE# MAC^^2•010 ALPHA* 5.000 


PANEL 

VORTEX 

AXIAL 

LATERAL 

VERTICAL 

NO. 

STRENGTH 

VELOCITY 

VELOCITY 

VELOCITY 

i 

•22608 

>.00993 

.01805 

•09163 

z 

.20022 

. 12328 

>•18020 

>•03020 

3 

•16902 

.09 166 

>• 13526 

>.05650 

A 

•12030 

.06533 

>.09911 

>.07007 

s 

•07135 

•06635 

>•07099 

-.08601 

6 

•03608 

.03523 

^.05783 

-.10095 

7 

•03636 

.06023 

>.05336 

-•11633 

s 

•06207 

•05516 

>.05672 

-.12661 

9 

.07180 

• 06669 

>•06097 

-.13338 

10 

•08669 

•06911 

>•05953 

-.13619 

1 1 

•11016 

.07769 

>•06001 

>•13626 

12 

. 16587 

>.06396 

.05506 

•09163 

13 

.16609 

.05036 

>•05625 

>•03020 


•17615 

.08086 

>.09951 

-.05450 

15 

• 16986 

.10855 

>• 16900 

-.07007 

16 

•15922 

.09612 

12989 

>•08601 

17 

.13567 

.09057 

>•12563 

>• 10095 

18 

.10363 

•08163 

11326 

-.11633 

19 

.06977 

.07059 

-.10269 

-.12661 

20 

.05319 

• 06391 

>•09506 

-.13338 

21 

•06661 

.05768 

>.08707 

-.13619 

22 

. 06233 

.05383 

>.00265 

-.13626 

23 

•16925 

-.05133 

•05979 

•09163 

26 

• 16901 

.03965 

-.06683 

-.03020 

25 

•15071 

.05676 

>.06657 

-.05650 

26 

.15332 

.06729 

-.07965 

-.07007 

27 

•15676 

.07891 

>.09635 

-.08601 

26 

.16070 

.10322 

129 10 

-.10095 

29 

• 16060 

.12730 

-. 16727 

-.11633 

30 

.15660 

.12167 

>. 15367 

-.12661 

31 

• 16868 

.11695 

-.16177 

-.13338 

32 

.13225 

.10889 

-.13537 

-.13619 

33 

. 11211 

.09716 

>.12656 

-.13626 

34 

•16331 

>.05763 

•06712 

.09163 

35 

.16362 

.03606 

>•03833 

-.03020 

36 

.16370 

.05333 

>•06612 

-.05650 

37 

.16626 

.06617 

-.07866 

-.07007 

38 

•16676 

.07211 

-.08825 

-•08601 

39 

.16566 

.08331 

>• 10167 

-.10095 

60 

. 16666 

.09269 

>. 11202 

-•11633 

61 

•16815 

.09839 

-.11893 

— • 12661 

62 

.16990 

. 10307 

-.12600 

-•13338 

43 

.15190 

•10606 

>• 13016 

-.13619 

44 

.15398 

.10761 

-.12978 

-.13626 

65 

.16022 

-.07167 

.08951 

•09163 

66 

•16038 

.02268 

-.02033 

-.03020 

47 

•16056 

.05080 

>, 06167 

-.05650 

60 

.16085 

.05986 

>.07316 

-.07007 

69 

• 16119 

.06769 

>•08261 

-.08601 

50 

•16151 

•07636 

-.09668 

-.10095 

51 

.16182 

•08737 

>• 10503 

-.11633 

52 

.16215 

.09669 

-.11528 

-.12661 

53 

.16254 

.09951 

>•12328 

-.13338 

56 

.16289 

.09895 

>. 125&6 

-.13619 

55 

.16325 

.09913 

>. 13002 

-.13626 
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□Give CVtINOEA BOOy 
SINGULARITY PANELING 


with 45 DEGREE SNEEP 
FOR sample case 


NACA 65A004 HlD-i4lNG 


INTEGRATION OF THE PRESSURE DISTRIBUTION 
ON THE NiNG UPPER SURFACE 


POINT 


MACH- 2.0100 
ALPHA. 5.0000 

X Y 


x/c 


1 

16.76499 

2.30734 

0.00000 

2 

17.61117 

2.30734 

0.00000 

3 

18.45734 

2.30734 

0.00000 

4 

19.30352 

2.30734 

0.00000 

5 

20.14970 

2.30734 

Q. 00000 

6 

20.99587 

2.30734 

0.00000 

7 

21.84205 

2. 30734 

0.00000 

8 

22 .68823 

2.30734 

0.00000 

9 

23.53441 

2.30134 

0.00000 

10 

24.38058 

2.30734 

0.00000 

11 

18.82577 

4.12568 

0.00000 

12 

19.55072 

9.12568 

0.00000 

li 

20.27568 

4. 12568 

0.00000 

14 

21.00063 

4.12568 

0.00000 

15 

21 .72559 

4.12568 

0.00000 

16 

22.45054 

9. 12568 

0.00000 

17 

23.17550 

4.12568 

0.00000 

18 

23.90045 

9.12568 

0.00000 

19 

24.62541 

4.12568 

0.00000 

20 

25.35036 

4. 12568 

0.00000 

21 

21,51106 

6.49507 

0.00000 

22 

22.07808 

6.49507 

0.00000 

23 

22.64507 

6.49507 

O.GOOOO 

24 

23.21207 

6.49507 

0.00000 

25 

23. 77906 

6.99507 

0.00000 

2 6 

24.34606 

6.99507 

0.00000 

27 

24.91305 

6.99507 

0.00000 

28 

25.48005 

6.49507 

0.00000 

29 

26.04704 

6.49507 

0.00000 

30 

26.61404 

6.49507 

0.00000 

3 1 

24. 16649 

8.83608 

0.00000 

32 

24.57728 

8.83808 

0.00000 

33 

24.98808 

8.83808 

0.00000 

34 

25.39887 

8.83808 

0.00000 

35 

25.80967 

8.83808 

0.00000 

36 

26.22046 

8.B3808 

0.00000 

37 

26.63126 

8.83608 

0.00000 

38 

27.04205 

8.83808 

0.00000 

39 

27.45285 

8.63808 

0.00000 

40 

27.86364 

6. 83806 

0.00000 

41 

26.58702 

10.97384 

0.00000 

42 

26.85543 

10.97384 

0.00000 

43 

27.12384 

10.97384 

0.00000 

44 

27.39225 

10.97384 

0.00000 

45 

2 7. 66 066 

10.97384 

0.00000 

46 

27,92907 

10.97384 

0.00000 

47 

28.19748 

10.97384 

0.00000 

48 

28.46589 

10.97364 

0.00000 

49 

28. 73430 

10.97384 

0.00000 

50 

29.00271 

10.97384 

0.00000 


.05000 

.15000 

•25000 

.35000 

.45000 

.55000 

.65000 

.75000 

•B5000 

-V5000 

.05000 

. 15000 

.25000 

.35000 

.45000 

.55000 

.65000 

.75000 

.85000 

.95000 

•05000 

. 15000 

.25000 

.35000 

.45000 

.55000 

.65000 

.75000 

.85000 

.95000 

.05000 

.15000 

.25000 

.35000 

.45000 

•55000 

.65000 

• 75000 

•85000 

.95000 

.05000 

.15000 

.25000 

.35000 

.45000 

.55000 

.65000 

.75000 

•85000 

.95000 


2Y/8 

Z/C 

. 19228 

0.00000 

.19228 

0.00000 

.19228 

0.00000 

•19228 

0.00000 

•19228 

0.00000 

.19228 

0.00000 

.19228 

0.00000 

.19228 

0.00000 

.19228 

0.00000 

.19228 

0.00000 

.3438 1 

0.00000 

• 34381 

0.00000 

.34381 

0.00000 

.34381 

0.00000 

• 3438 1 

0.00000 

.34361 

0.00000 

• 34381 

0.00000 

.34381 

0.00000 

•34381 

0.00000 

•34381 

0.00000 

•54126 

0.00000 

• 54126 

0.00000 

• 59126 

0.00000 

•54126 

0.00000 

•54126 

0.00000 

.54126 

0.00000 

.54126 

0.00000 

.54126 

0.00000 

.54126 

0.00000 

•54126 

0.00000 

• 73651 

0.00000 

. 73651 

0. 00000 

.73651 

0.00000 

. 73651 

0.00000 

• 73E51 

0.00000 

•73651 

0.00000 

.73651 

0.00000 

• 73651 

0.00000 

•73651 

0.00000 

.73651 

0.00000 

• 91449 

0.00000 

.91449 

0.00000 

.91449 

0.00000 

.91449 

0.00000 

.91449 

0.00000 

.91449 

0.00000 

.91449 

0.00000 

.91449 

0.00000 

.91949 

0.00000 

.91449 

0.00000 


CP 

CN 

-.10958 

.12071 

-.18922 

.20844 

-.14170 

.15610 

-.09971 

• 10964 

-.07159 

•07887 

-.06742 

.07427 

-.08517 

.09383 

-.10645 

•11727 

-.11814 

.13014 

-.12830 

.14133 

-.01494 

.02588 

-.11757 

.20373 

-.16626 

.28810 

-.17911 

.31036 

-.16501 

.28593 

15367 

•26627 

-.13760 

.23844 

-.12315 

•21339 

-.11204 

.19414 

-.10303 

.17853 

.00257 

-.00342 

-.08677 

.11535 

-.11058 

.14701 

-.12933 

•17194 

-.15894 

.21131 

-.19700 

•26191 

- .21101 

.28052 

-.20309 

.27001 

-.19389 

.25777 

-.17989 

•23916 

.01376 

-.01323 

-.07874 

.07570 

-.10507 

• 10102 

-.12121 

.11654 

-.13736 

•13206 

-.15448 

.14852 

— • 1 6686 

.16043 

-.17540 

.16863 

-.18025 

.17330 

-•18353 

.17646 

.03823 

-.01913 

-.06559 

.03282 

-.09898 

.04952 

-.11357 

•05682 

-. 12942 

.06475 

-.14607 

• 07308 

-.15978 

.07994 

-.17000 

.08505 

-.17392 

.08702 

-.17421 

.08716 


cr 

-.01423 

-.00934 

-.00388 

-.00111 

•00042 

.00160 

•00322 

.00502 

.00607 

•00665 

-.00305 

-.00913 

-.00716 

-.00314 

.00152 

.00572 

.00818 

.00914 

.00905 

•00640 

.00040 

-.00517 

-.00366 

-.00174 

.00112 

.00563 

•00963 

•01157 

.01202 

.01126 

.00156 

-.00339 

-.00251 

-.00118 

.00070 

•00319 

.00551 

.00722 

.00808 

.00831 

.00225 

-.00147 

-.00123 

-.00057 

.00034 

•00157 

.00274 

.00364 

.00406 

.00410 


CM 

. 48864 
.66740 
.36772 
. 16580 
.05231 
-.01358 
-. 09655 
-.21991 
-.3541 7 
-.50421 
.05143 
.25717 
• 15480 
-.05823 
-.26094 
-.43603 
-.56330 
-.65883 
-.74012 
-.81004 
.00239 
14593 
-.26934 
-.41251 
-.62676 
-.92533 
-i. 15017 
-1.26013 
-1. 34916 
-1.38737 
.04436 
-.28496 
-.42175 
-.53442 
-.65988 
-.80313 
93341 
-1.05041 
-1. 15065 
-1.24413 
.11044 
-. 19831 
-.3 1252 
-.37386 
-.44341 
-.52006 
-.59032 
-.65091 
-.68928 
-.71385 


POINT 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
2/ 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 
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VELOCITIES ON KING LOWER SURFACE, RACH-2.010 ALPHA- 5.000 


PANEL 

NO. 


VORTEX AXIAL 

STRENGTH VELOCI TV 


LATERAL 
VELOC ITY 


VERTICAL 

VELOCITY 


i 

•22408 

-.23401 

z 

.20022 

-.07694 

3 

• 16902 

-.07758 

A 

.12030 

-.05497 

5 

.07135 

-.02700 

6 

.03408 

•00115 

7 

.03636 

.00387 

8 

•06207 

-.00691 

9 

.07180 

-.00731 

10 

.08669 

-.01757 

1 1 

. 11014 

-.03265 

12 

.16587 

-.20984 

13 

.16609 

-.11573 

14 

. 17415 

-.09329 

15 

.16984 

-.06128 

16 

.15922 

-.06311 

17 

.13547 

-.04490 

18 

. 10343 

-.02179 

19 

.06977 

.00082 

20 

.05319 

.01072 

21 

•04661 

.01107 

22 

.04233 

.01131 

23 

.14925 

-.20038 

24 

. 14901 

-.10955 

25 

.13071 

-.09 395 

26 

.13332 

-.08603 

2 7 

.15678 

-.07786 

28 

.16070 

-.05748 

29 

• 16040 

-.03310 

30 

.156 40 

-.03473 

31 

•14848 

-.03154 

32 

.13225 

-.02336 

33 

.11211 

-.01493 

34 

•14331 

-.20075 

35 

.14342 

-.10936 

36 

.14370 

-.09038 

37 

.14424 

-.08007 

38 

.14474 

-.07263 

39 

.14546 

-.06215 

40 

• 14666 

-.05397 

41 

.14815 

-.04976 

42 

.14990 

*• 04683 

43 

.15198 

-.04792 

44 

.15398 

-.04636 

43 

. 14022 

-.21168 

46 

.14038 

-.11769 

47 

. 14054 

-.08974 

48 

.14085 

-.08102 

49 

. 14119 

-.07349 

50 

. 14151 

-.06315 

51 

• 14182 

-.03445 

32 

.14215 

-.04746 

53 

. 14254 

-.04303 

54 

.142 89 

-.04393 

55 

.14325 

-.04412 


.27S48 
.03416 
•06367 
.03238 
.03370 
.01637 
.02286 
.038 36 
.04059 
.03057 
.06299 
. 24637 
.13732 
.10283 
. 04903 
•03825 
,04213 
.02868 
.01455 
.01115 
.01517 
.01750 
.23392 
.12903 
. 10910 
.09883 
.08716 
•03580 
.01739 
.02806 
.03468 
.03133 
.02941 
.23432 
• 12899 
. 10330 
.08932 
.08038 
.06778 
.05819 
.05237 
.04647 
.04356 
•04500 
.25309 
. 14344 
. 10248 
.09111 
.08196 
•06996 
.06007 
•03006 
.04231 
.03996 
.03398 


-.26595 

-.14412 

-.11981 

-.10424 

-.09031 

-.07336 

-.05798 

-.04770 

-.04093 

-.04012 

-.04006 

-.26593 

-. 14412 

-.11981 

-.10424 

-.09031 

-.07336 

-.05798 

-.04770 

-.04093 

-.04012 

-.04006 

-.26595 

-.14412 

-.11981 

-. 10424 

-.09031 

-.07336 

-.05798 

-.04770 

-.04093 

-.04012 

-.04006 

-.26595 

-•14412 

-.11981 

-.10424 

-.09031 

-.07336 

-.05798 

-.04770 

-.04093 

-.04012 

-.04006 

-.26595 

-.14412 

-.11981 

-. 10424 

-.09031 

-.07336 

-.05798 

-.04770 

-.04093 

-.04012 

-.04006 


308 



OGIVE CYLiKOEft UOOY M I TH 4^ DEGREE SWEEP NACA 65A00<» MIO-NlNG 


SINGULARCTY PA^ELING FOR SAMPLE CASE 

INTEGRATION OF TFE PRESSURE DISTRIBUTION 
ON THE WING LOmEC SURFACE 


MACH* ^.0100 
ALPHA- 3.0000 


POINT 

X 

Y 

1 

x/c 

2Y/B 

1 

16.764 99 

2.30734 

0.00000 

.03000 

.19228 

Z 

17.61117 

2.3C734 

0.00000 

.13000 

.19228 

3 

18.43734 

2.30734 

0.00000 

.23000 

.19228 

4 

19.30332 

2.30734 

0.00000 

.33000 

.19228 

5 

20.14970 

2.30734 

0.00000 

.43000 

.19228 

6 

20.99387 

2.30734 

0.00000 

.35000 

. 19228 

1 

21.84203 

2.30734 

0.00000 

.63000 

. 19228 

8 

22.68823 

2.30734 

0.00000 

.73000 

.1922 8 

9 

23.33441 

2.30734 

o.coooo 

.83000 

. 19228 

10 

24.38038 

2.30734 

0.00000 

.93000 

. 1922 8 

11 

10.82377 

4.12368 

0.00000 

.03000 

.34381 

12 

19.33072 

4.12568 

0.00000 

.13000 

•34381 

13 

20.27368 

4. 12368 

0.00000 

.23000 

.34361 

14 

21.00063 

4.12368 

0.00000 

.33000 

.34381 

13 

21.72339 

4.12568 

0.00000 

.43000 

.34381 

16 

22.45054 

4.12368 

0.00000 

.55000 

.34381 

17 

23.17350 

4.12368 

0.00000 

.63000 

.34381 

18 

23.90043 

4.12368 

0.00000 

.75000 

.34381 

19 

24.62341 

4.12368 

0.00000 

.83000 

.34381 

20 

23.33036 

4.12368 

0.00000 

*95000 

.3438 1 

21 

21.31108 

6.49507 

0.00000 

.05000 

.54126 

22 

22.07808 

6.49507 

0,00000 

.13000 

.34126 

23 

22.64307 

6.49307 

0.00000 

.23000 

.54126 

24 

23,21207 

6.49307 

0.00000 

.33000 

•54126 

23 

23,77906 

6.49507 

0.00000 

.43000 

.54126 

26 

24.346C6 

6.49307 

0.00000 

.33000 

.54126 

27 

24.91303 

6.49507 

0.00000 

.63000 

.34126 

28 

23.48003 

6.49507 

0.00000 

.73000 

.54126 

29 

26.04704 

6.49507 

0.00000 

.83000 

.34126 

30 

26.61404 

6.49307 

0.00000 

.93000 

.34126 

31 

24.16649 

8.83808 

0.00000 

.03000 

.73651 

32 

24.37728 

8.83808 

0.00000 

.13000 

.73631 

33 

24.98808 

8.83808 

0.00000 

.23000 

.73651 

34 

25.39887 

8.83808 

0.00000 

.35000 

.73651 

33 

23.80967 

8.83808 

0.00000 

.43000 

.73631 

36 

26.22046 

8.83808 

0.00000 

.53000 

.73631 

37 

26.63126 

8.83808 

0.00000 

.65000 

.73651 

38 

27.04203 

8.83808 

0.00000 

.73000 

.73651 

39 

27.43285 

8.83808 

0.00000 

•63000 

.73631 

40 

27.86364 

8.83808 

0.00000 

.93000 

. 73651 


26.38702 

10.97384 

0.00000 

.03000 

•91449 

42 

26.85343 

10.97384 

0.00000 

.13000 

.91449 

43 

27. 12384 

10.97384 

0.00000 

.25000 

.91449 

44 

27.39223 

10*97384 

0.00000 

•35000 

.91449 

43 

27.66 066 

10.97384 

0.00000 

.45000 

.91449 

46 

27.92907 

10.97384 

0.00000 

.53000 

•91449 

47 

28.19748 

10.97384 

0.00000 

.63000 

.91449 

48 

28.46389 

10.97384 

0.00000 

.75000 

•91449 

49 

28. 73430 

10.97384 

0.00000 

.83000 

•91449 

30 

29.00271 

10.97384 

0.00000 

.95000 

.91449 


Z/C 

CP 

CN 

CT 

CM 

POINT 

0.00000 

.29517 

.32316 

.03833 

1.31626 

1 

0.00000 

.17381 

.19147 

•00838 

.61303 

2 

0.00000 

•14939 

•16478 

•00410 

.38817 

3 

0.00000 

.09334 

•10304 

•00104 

. 15554 

4 

0.00000 

.03394 

.03738 

-.00020 

•02480 

5 

0.00000 

.00168 

.00163 

-.00004 

-.00034 

6 

0.00000 

.00837 

.00944 

-.00032 

-.00971 

T 

0.00000 

.01863 

.02054 

-.00088 

-.03852 

8 

0.00000 

.02880 

•03172 

-.00148 

-.08633 

9 

0.00000 

•03438 

•05990 

-.00282 

-.21371 

10 

0.00000 

.31460 

.54313 

•06426 

1.08331 

11 

0.00000 

•22633 

.39257 

.01739 

.49553 

12 

0.00000 

.17245 

•29882 

•00743 

.16036 

13 

0.00000 

,14107 

•24443 

•00247 

-.04587 

14 

0.00000 

.12247 

.21221 

-.00113 

-.19366 

15 

0.00000 

.07674 

.1329 7 

-.00286 

-.21774 

16 

0.00000 

.02777 

.04812 

-.00165 

-.11367 

17 

0.00000 

-.00384 

-.0101 L 

•00043 

•03122 

18 

0.00000 

-.01631 

-.02825 

•00132 

.10772 

19 

0.00000 

-.01720 

-.02980 

.00140 

.13520 

20 

0.00000 

.30197 

.40146 

.04732 

-.26023 

21 

0.00000 

.22019 

.29274 

•01312 

-.37033 

22 

0.00000 

.19751 

•26238 

•00653 

-.48107 

23 

0.00000 

•18127 

.24099 

.00244 

-.57815 

24 

0.00000 

.15125 

.20109 

-.00107 

-.59643 

25 

0.00000 

.10311 

.13708 

-.00295 

-.48431 

26 

0.00000 

.07770 

.10330 

-.00334 

-.42354 

27 

0.00000 

•07472 

.09934 

-.00426 

-.46362 

28 

0.00000 

.06182 

.08219 

-.00383 

-.43018 

29 

0.00000 

.04414 

.03868 

-.00276 

-.34040 

30 

0.00000 

•30182 

.29018 

•03420 

-.97310 

31 

0.00000 

•21633 

.20799 

.00932 

-.78293 

32 

0.00000 

•18771 

•18047 

•00449 

-.75347 

33 

0.00000 

.16930 

•16296 

.00163 

-.74733 

34 

0.00000 

•L3022 

. 14442 

-.00077 

-.72164 

35 

0.00000 

.12932 

.12432 

-.00266 

-.67336 

36 

0.00000 

.11537 

.11092 

-.00381 

-.64536 

37 

0.00000 

.10711 

. 10298 

-.00441 

-.64146 

38 

0.00000 

.10512 

.10106 

-.00471 

-.67103 

39 

0.00000 

.10462 

.10039 

-.00473 

-.70920 

40 

0.00000 

.31643 

. 13832 

.01666 

-.91414 

41 

0.00000 

•22349 

.11182 

.00301 

-.67565 

42 

0.00000 

.18804 

•09408 

•00234 

-.59373 

43 

0,00000 

.17134 

.08373 

•00087 

-.56403 

44 

0.00000 

.13208 

•07609 

-.00041 

-. 52104 

45 

0.00000 

.13a93 

.06331 

-.00141 

-.46617 

46 

0.00000 

.11329 

.03668 

-.00195 

-.41859 

47 

0.00000 

•10032 

.03029 

-.00213 

-.38490 

48 

0.00000 

•09666 

•04836 

-.00226 

-.38310 

49 

0.00000 

•09814 

.04910 

-.00231 

-.40212 

50 
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TOTAL COEFfiCUNTS 


ON THE MING 


REFA« 

144.0000 

REFB* 

12.0000 

REFX« 

20.8130 

REFZ- 

0.0000 

CN- 

.1969 



CT» 

.0046 



CM- 

-.0705 



CL- 

.1957 



CO* 

• 0217 



XCP* 

-.3600 




TOTAL COEFFICIENTS 

ON THE CCHRLETE CONFI GOR AT ION 


REFA* 

144.0000 

REF8* 

12.0000 

REFC 

REFX* 

20.8130 

REFZ* 

0.0000 


CN- 

.2495 




CT* 

.0081 




CH- 

-.0651 




CL* 

-24 79 




CO* 

• 0298 




tCP- 

-.2628 





6*8900 


6*8900 
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SECTION COEFFICIENTS 


ON THE WING 


OELY- 

CN» 

CT- 

CM- 

CL- 

CO* 

XCP* 


1.3030 

.1974 
.0037 
.0336 
. 1963 
.0209 
• 1612 


REFL* 


6.8900 XLE- 


OELY- 

2.4000 

CN* 

.2305 

CT* 

• 0063 

CM* 

-.0135 

CL* 

.2291 

CD- 

• 0263 

XCP* 

-.0590 


REFL- 


6. 8900 XLE* 


DELY* 

2.3600 

CN- 

•2863 

CT* 

• 0069 

CM- 

-.1299 

CL- 

.2846 

CO- 

• 0318 

XCP* 

-.4563 


REFL* 


6.8900 XLE* 


SECTION CaeFFlCIENTS 
ON THE WING 


DELV- 

2.3700 

CN* 

.2841 

CT* 

.0058 

CM* 

-.2140 

CL* 

• 2825 

CD* 

• 0305 

XCP- 

-.7577 


OELY* 

1.9000 

CN- 

.2732 

CT- 

• 0062 

CM* 

-.2/62 

CL* 

.2716 

CD* 

• 0300 

XCP* 

-1.0171 


CPSTAG • 2.45650 CPCRIT * 

TIHE * 303.80100 


REFL* 6.8900 XLE- 


REFL* 6.8900 XLE* 


1.13092 CPWAC * -.35360 


16.3419 


18.4633 


21.2276 


23.9611 


26.4528 
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